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foreword 


This quarterly review of reactor development has been prepared at the request 
of the Division of Information Services of the U. S. Atomic Energy Commission. 
Its purpose is to assist interested organizations in the task of keeping abreast 
of new results in reactor technology for civilian application. 

The report is a concise discussion of selected phases of research and de- 
velopment for which there have been significant advances or a heightened inter- 
est in the past few months. It is not meant to be a comprehensive abstract of all 
material published during the quarter, nor is it meant to be a treatise on any 
part of the subject. The intention is to cover the various areas of reactor de- 
velopment from the general viewpoint of the reactor designer rather than from 
the more detailed points of view of specialists in the individual areas. However, 
papers which are thought to be of particular significance or particular useful- 
ness in specialized fields will be mentioned in short notes. In the over-all plan 
of the report, it is intended that various subjects will be treated from time to 
time and will be brought up to date at that time. 

Any interpretation of results which is given represents only the opinion of the 
editors of the report, who are General Nuclear Engineering Corporation person- 
nel. Readers are urged to consult the original references in order to obtain all 
the background of the work reported and to obtain the interpretation of the re- 
sults given by the original authors. 


W. H. ZINN 
General Nuclear Engineering Corporation 
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REACTOR PHYSICS 





Thermal-neutron Distributions 
Near Temperature Discontinuities 


For some time Hanford has been investigating, 
both theoretically and experimentally, the dis- 
tributions of thermal neutrons, in space and 
- energy, near a temperature discontinuity in a 
moderating medium. These investigations are 
interesting for the additional insight they give 
into the process of neutron thermalization, and 
they have practical implications for some 
power reactors, for example, those employing 
solid moderators in which there may be sizable 
variations of temperature with position. 

The theoretical problem has been treated by 
Kottwitz, in a number of Hanford quarterly 
reports, !-5 and, more recently, by Whitehead.° 
The results of an experimental investigation 
have been reported by Bennett). These re- 
sults will be summarized briefly. 

A column of graphite, about 2 ft square in 
cross section, was installed in a closely fitting, 
vacuum tank; the tank was inserted horizontally 
into a closely fitting square hole in the Physical 
Constants Testing Reactor (PCTR), where it 
was surrounded by the reactor fuel rods. The 
square graphite column was divided longitudi- 
nally into three rectangular columns which lay 
side by side, separated by gaps. The central 
column of the three could be either heated or 
cooled. The evacuated gaps insulated the cen- 
tral column from the two adjacent columns, 
and moderator temperature discontinuities were 
established at the gaps. Fission neutrons from 
surrounding fuel elements were thermalized in 
the graphite columns, and a current of thermal 
neutrons flowed from the columns to the walls 
of the vacuum tank and its surroundings. The 
spatial variation of thermal-neutron spectrum 
produced by the temperature discontinuity was 


inferred from the spatial distribution of activa- 
tion produced in 1/v detector foils (copper). 

Figure 1 shows traverses of the subcadmium 
activation of the copper foils, along a line from 
the center of the central (heated or cooled) 
graphite column, across the insulating gap, and 
to the outer edge of one of ** > outer columns. 
In the second curve from .ue top, all three 
graphite columns are at one temperature, and 
the activation curve follows the normal dis- 
tribution of thermal-neutron flux, as estab- 
lished by the production of thermal neutrons 
in the graphite and their leakage to the wall. 
When the temperature of the central graphite 
column is raised or lowered, respectively, the 
activation of the foils in the central region be- 
comes lower or higher relative to that of the 
foils in the outer graphite, because of the 1/v 
variation of the activation cross section. 

In analysis of the experiment, a representa- 
tion proposed by Selengut® was applied to fit 
the observed activation distributions. For the 
case of two regions at two different uniform 
temperatures T, and T, (as in the experiment), 
it is assumed that all thermal neutrons lie in 
one or the other of two energy groups and 
that the two energy groups are Maxwellian 
distributions of neutrons, with characteris- 
tic temperatures T; and T,. All thermal neu- 
trons born, by slowing down, in region 1 are 
considered to be born into the group char- 
acterized by 7;, whereas all neutrons born in 
region 2 are considered to be born into the 
group characterized by T,. Any neutrons cross- 
ing the boundary between the two regions are 
considered to be subject to transfer into the 
energy group characteristic of the new region, 
with a probability specified by a rethermaliza- 
tion cross section which is analogous to the 
intergroup slowing down cross section of the 
conventional multigroup formulation. 
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RADIAL POSITION IN UNITS OF 0.25 IN. 


Figure 1— Experimental thermal-activation traverses ingraphite with a temperature discontinuity.® 


Thus, mathematically, the distribution of 
neutron flux in space and energy can be ex- 
pressed as a sum of two fluxes, $,(x) and 
$,(x), having the Maxwellian energy distribu- 
tions M(E,T;) and M(E,T,) 


$(E,x) = b4(x) M(E,T;) + b,(*) M(E,7) 


The distribution of ¢,; and ¢, in space can be 
determined by solving the two-group equations 
for the appropriate source distribution and 
geometry. In region 1, for example, the two- 
group equations would be 


—D, Vv? $,+Za, 1 = 212 do +¢ 
—D,V? oy + (Zaz + X42) $2 = 0 


where Dis the diffusion constant, Zq_ is the 
(macroscopic) absorption cross section, Dy is 
the rethermalization cross section for transfer 
of neutrons from group 2 to group 1, andgq is 
the source density of thermal neutrons. 

For a foil with a 1/v activation cross sec- 
tion, the effective activation cross section ina 
neutron flux having a Maxwellian energy dis- 


tribution characterized by the temperature T 
will be proportional to T~%. The spatial dis- 
tribution of foil activation, A(x), arbitrarily 
normalized, will then be given in terms of 4, 
and ¢, as 


A(x) = VT,/T, $4(x) + $(x) 


In analyzing the experimental results, the 
thermal source distribution, g, was first de- 
termined by finding a source distribution which 
would fit the activation distribution measured 
when all the graphite was at the same tem- 
perature. Using this source distribution, each 
of the cases involving a temperature discon- 
tinuity was analyzed, and a value of the re- 
thermalization cross section was sought that 
would give flux distributions, $,(x) and ¢,(x), 
consistent with the observed activation dis- 
tribution. In each case it was assumed that the 
rethermalization cross section 24) for transfer 
from group 2 to group 1 was equal to the cross 
section &,,;, for the reverse transfer, although 
the reference points out that this equality is 
not necessarily true. 
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It was found that the apparent value of the 
rethermalization cross section varied from 
experiment to experiment. Theoretically, the 
rethermalization cross section can be expressed 
in terms of the scattering cross section, Zs, 
and the ratio, u, of the neutron mass to the ef- 
fective mass of the moderator atom:? 


Zret = 2uzs 


In Table I-1 are listed, for four cases meas- 
ured, the values of the rethermalization cross 


temperature was nearly linear over this range, 
Ap/AT varying from about 1.4 x 10“/°C at 
20°C to about -—0.5x 1074/°C at 90°C. The 
temperature at which the coefficient reached 
zero and changed sign from positive to nega- 
tive was about 72°C. 

The major difficulty in the theoretical cal- 
culation was considered to be the geometrical 
complexity of the assembly. Calculations based 
on three different geometrical models were 
compared with the measured values. One of 
these models could be handled by one-dimen- 











Table I-1 RETHERMALIZATION PROPERTIES OF GRAPHITE® 
Temp., °K 

Experiment Region 1 Region 2 Macroscopic Relaxation Effective 

No. (T;) (T) cross section length mass 
(Z ret), 1079 em (L ret), Cm (Meq), amu 

1A 108 + 2 285 +2 1.9 +0.5 21.0 + 2.8 420 + 110 
2A 491+2 320 + 2 6.4+1.5 11.4+1.4 125 + 30 
2B 491 +2 320 + 2 3.8 + 0.9 14.8 + 1.7 208 + 47 
3A 666 + 2 350 + 2 14.5 + 2.6 7.5 + 0.7 55 +10 





Yq= 2.822 x 1074 em (macroscopic absorption cross section of graphite). 


D = 0.8086 cm (graphite diffusion coefficient). 


p = 1.66 g/cm? (graphite density). 


section, yet; the values of the rethermaliza- 
tion relaxation length, Lre, = VD/Zret ; and the 
corresponding values of the effective mass of 
the graphite “atom’’ as a scatterer. 


Temperature Coefficients 
of Reactivity 


Although the temperature coefficient of re- 
activity in the temperature range near room 
temperature is not of the greatest importance 
for power reactors, a comparison of theoretical 
and experimental values for these coefficients 
constitutes a test of theory and is, therefore, 
of considerable interest. Reference 7 reports 
comparisons of this kind for a highly enriched 
assembly made up of plate type uranium- 
Zircaloy elements in a slab geometry, 7.4 by 
29.7 by 54 in. high. The slab was controlled by 
a centrally located row of four cruciform con- 
trol rods operated as a bank. The critical 
position of the bank of rods was approximately 
32 in. withdrawn. 

The temperature coefficient was measured 
over the temperature range from about 20 to 
90°C. The variation of the coefficient with 


sional calculations, whereas the other two re- 
quired two-dimensional calculations. Reason- 
able agreement was obtained with all three, but 
a two-dimensional calculation in the x-z direc- 
tions (« being the direction of the thin dimen- 
sion of the slab and z being the direction of the 
control-rod axes) was considered best. The 
two-group approximation was used in all cases, 
with Maxwell-averaged thermal constants. 
Although the reference describes the calcula- 
tion methods in some detail, it does not at- 
tempt to account for the results in terms of 
the various physical processes involved. Such 
an analysis would be interesting, since it seems 
rather surprising that an assembly of such 
high neutron leakage would demonstrate a posi- 
tive temperature coefficient of reactivity at 
temperatures up to 72°C. The assembly did 
contain a rather large fraction of its water in 
the control-rod channels, and it is possible 
that a high thermal-flux peaking in the chan- 
nels may account for much of the positive 
component of the temperature coefficient. 
Measurements of the temperature coefficients 
of reactivity in several highly enriched water- 
moderated assemblies at elevated temperatures 
(up to about 540°F) are reported in reference 
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8. The measurements were made in the Pres- 
surized Test Reactor (PTR) at Knolls Atomic 
Power Laboratory (KAPL). The assemblies 
were relatively uncomplicated, and a simplified 
analysis was applied in which it was assumed 
that the temperature dependence of reactivity 
arises only from changes in L’ and 7, the ther- 
mal-diffusion area and the Fermi age. It was 
found that this approximation was good enough 
to support predictions of the reactivity varia- 
tion for a given assembly over a range of 
temperatures, based on measurements of the 
temperature coefficient of reactivity at two 
different temperatures. However, neglect of 
the change of geometrical buckling of reflected 
reactors when the temperature is changed ap- 
pears to be questionable. 


Measurements Related 
to Thermal Breeding 


In the last issue of Power Reactor Technol- 
ogy, 3(2): 21 (March 1960), a summary was 
given of the Conference on the Physics of 
Breeding, at Argonne National Laboratory, 
Oct. 19-21, 1959. Two recent reports from 
Oak Ridge National Laboratory (ORNL) give 
more detailed accounts of two of the investiga- 
tions discussed at the Conference. 


Reference 9 reports the measurements of 7 
(neutrons produced per neutron absorbed) for 
thermal-neutron irradiation of U?** and U?** by 
the manganese-bath method. In concept, the 
method is simple and direct. A beam of ther- 
mal neutrons is directed through a tube to the 
center of a large spherical tank filled with an 
aqueous solution of MnSO,. After a sufficient 
time the neutron beam is stopped, the MnSO, 
solution is thoroughly mixed, and the manga- 
nese activity is determined by counting a sample 
of the solution. A target of U’*’ or U* is then 
installed at the center of the sphere where the 
neutron beam will impinge upon it. The sample 
is sufficiently thick to absorb very nearly all 
the incident neutrons, and the spherical tank 
is so large that very nearly all the fission neu- 
trons produced in the uranium are slowed down 
and absorbed in the MnSO, solution. For equal 
neutron beam strengths (as determined by 
monitors) and equal irradiation times, 7 is ob- 
tained (except for small corrections) from the 
ratio: (manganese activity produced when the 


uranium is present)/(manganese activity pro- 
duced when the uranium is absent). 

The experimental technique and the correc- 
tions are discussed in the reference in con- 
siderable detail. The results, corrected to 
2200 m/sec values, are given in Table I-2. 


Table I-2 THERMAL 7 VALUES DERIVED FROM 
MANGANESE BATH MEASUREMENTS® 
(2200 M/SEC VALUES) 








Material n 
7233 2.296 + 0.010 
y235 2.077 + 0.010 


Natural uranium 1.329 + 0.009 





Also given is 7 for natural uranium, based on 
the measured value of 7 for U* and on recent 
determinations of the absorption cross sec- 
tions of U** and vU*** (681+ 7 barns’® and 
2.71 + 0.02 barns,'! respectively). 

The status of critical experiments with aque- 
ous solutions of fissionable isotopes, for the 
determination of 7 of U***, is reported in refer- 
ence 12. Analyses of the data were not finished 
at the time of publication of the reference. The 
measurement judged to give the best measure- 
ment in spherical geometry gave n = 2.295+ 
0.036 for U**?, whereas measurements in cy- 
lindrical geometry gave values from 2.281 + 
0.016 to 2.290+0.025. These are evidently 
the effective 7 values for the neutron spectra 
in the critical assemblies. 

Reference 13 gives the details of the calcula- 
tion of the “fast effect’’ in beryllium oxide, 
which was reported at the Physics of Breeding 
Conference. The net fast effect, «€, analogous, 
in the neutron balance, to the fast fission effect 
in U*8, is calculated to be between 1.027 and 
1.046. This net effect includes the neutron pro- 
duction by the (”,2 7) reaction and the neutron 
loss by the (n,a) reaction but does not include 
the neutron loss to Li® manufactured in the 
BeO by the (”,q@) reaction. The latter effect is 
discussed also in the reference. 


Power Flattening by Variation 
of Moderator Properties 


The results of a calculational study of power 
flattening in a Sodium Graphite Reactor (SGR) 
by means of spatial variation of the moderator 
properties is reported in reference 14. It is 
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proposed that the properties of the moderator 
might be varied, as a function of position, 
either by varying the density of the graphite or 
by replacing some of the graphite by BeO. It 
is considered that the latter alternative might 
be accomplished either by mixing the BeO with 
the graphite at some stage in the manufacture 
or by placing solid BeO inserts in holes in the 
graphite. ; 

To flatten the normal power distribution in 
the sodium graphite reactor, the effective mod- 
erator-to-fuel ratio should be increased pro- 
gressively from the center of the core to the 
edges. This could be done either by decreasing 
the graphite density progressively toward the 
center or by increasing the BeO addition pro- 
gressively toward the edges, or by a combina- 
tion of the two variations. 

Table I-3 gives the characteristics of the 
sodium graphite reactor investigated. The. ref- 
erence gives calculated results for a number 
of different moderator distributions. An im- 


Table I-3 CHARACTERISTICS OF SGR CORE FOR POWER 
FLATTENING STUDY"! 





Dimensions 
Core height, ft 12.5 
Core diameter, ft 12.5 
Reflector thickness, in.: 
Axial 24 
Radial 17 
Moderator can size, in. 16.75 


Fuel Element 


Fuel rods per element 37 

Fuel rod diameter, in. 0.434 (10 g UO,/cem’) 
Fuel enrichment 3.5 at.% 235 
Process tube inside diameter, in. 4.00 


Unit Cell Volume Fractions 





Uranium dioxide 0.0466 
Sodium 0.0835 
Graphite 0.8510 
Stainless steel 0.0080 
Zirconium 0.0106 
Void 0.0003 
Two-group Constants 
Reflector 
Core Axial Radial 
Lay, em 0.000617 
ZR4,> em! 0.002368 0.00308 0.00322 
Dy, em 0.986 1.044 0.958 
Za,, cm! 0.005728 0.002764 0.001277 
€vZf.,cm  0,008521 
Dy, cm 0.887 0.949 0.879 





provement of 22 per cent in average specific 
power was shown for the case in which the 
power was flattened radially by incorporating 
40 vol.% of BeO in the peripheral annulus 
described by the outer 20 per cent of the core 
radius. An increase in average specific power 
of about 27 per cent was shown for the case in 
which the graphite density was varied radially 
from 80 per cent at the center to 100 per cent 
of normal density in the outer annulus. In- 
creases of up to 44 per cent in average specific 
power were shown for combinations of radial 
variations of graphite density and BeO content. 
Substantial improvements were also shown for 
axial variations of the moderator properties. 
The changes in moderator properties resulted 
in changes of reactivity of magnitude up to 
about 3 or 4 per centers. As would be ex- 
pected, the addition of BeO increased reactivity, 
whereas the reductior of graphite density de- 
creased reactivity. 


Analogue for Three- 
dimensional Reactor 


It has been reported!® that an analogue method 
was used for some of the three-dimensional 
flux-distribution determinations in connection 
with the design of the British Trawfynydd re- 
actor. The analogue is a thermal one which 
simulates the one-group reactor equation. The 
instantaneous temperature distribution in a 
cooling solid obeys a diffusion equation identical 
to the one-group critical equation for the neu- 
tron flux. In the analogue the reactor core is 
represented by a hot solid of similar shape 
which loses heat to the boundaries, and the 
temperature distribution in the solid is ob- 
served. The model is made ofa suitable wax. 
Spatial variations in the material buckling 
B’[=(k-1)/M?] are simulated by the use of 
waxes of different thermal diffusivities, where- 
as control rods and other lumped absorbers 
are simulated by cooled holes drilled in the 
wax. 
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Hot-channel Factors 


The concept of hot-channel factors has been 
discussed in the September 1959 and March 
1960 issues of Power Reactor Technology. In 
reference 1 the application of hot-channel hot- 
spot factors to the Engineering Test Reactor 
(ETR) was demonstrated, and the hot-channel 
calculations indicated a wall temperature at the 
hot spot of 395° F, some 30° above saturation. 
Reference 2 presents another method of taking 
into account the hot-channel factors and in- 
vestigates the calculations presented in refer- 
ence 1. 

The statistical method presented in reference 
2 is similar to that shown by Bonilla,’ in that 
the hot-spot wall temperature is written as the 
sum of two terms 


T wHC = Tw + No (1) 


where T,,7¢ = hot-spot wall temperature as- 
sociated with N number of 
sigmas 

» =the “nominal’’ wall temperature 

o = variance (standard deviation) 


ti 


U 


If N = 1, the probability that the actual hot-spot 
wall temperature exceeds Tc is 0.16; ifN = 
2, the probability is 0.023, etc. In reference 3, 
Ty is the wall temperature calculated without 
the application of any hot-channel factors. 
Sigma is calculated as the square root of the 
sum of the squares of the product of the indi- 
vidual fractional root-mean-square variation 
and the corresponding temperature rise. Al- 
though it is assumed that all hot-channel fac- 
tors correspond to two fractional root-mean- 
square variations, this would, of course, be 
modified by actual experience. For example, 
a recent reactor design report lists hot-channel 
factors that are three standard deviations from 
the design values.‘ 

The methods used in calculating T,, andoy in 
reference 2 are somewhat different. The usual 
calculation for Ty is defined by 


Ty = Ty + uA (2) 


where q/A is the heat flux and h is the heat- 
transfer coefficient at film temperature as cal- 
culated by some appropriate relation, e.g., the 
modified Colburn equation (h = bB of that equa- 
tion). This form is neither explicit in T,, (since 
h is a function of Ty and T; is a function of Ty) 
nor does it lend itself to a statistical error 
analysis. To circumvent this difficulty, utilizing 
the modified Colburn equation, triples of num- 
bers T,, Tp, and K*=(q/A)/B which solve 
Eq. 2 were determined, and by statistical re- 
gression analysis the equation 


T, = 46.68 + 0.877, + 16.99 x 10“K* (3) 


was derived which approximates to Eq. 2 within 
10°F in the range 250° F <= Tw = 550°F. Theap- 
plication of this equation to hot-spot temperature 
and associated uncertainty estimation is accom- 
plished by defining 7, and K* in terms of the 
basic components entering into their evaluation 
and through factors applied to the components 
where appropriate to translate them from the 
normal condition to the hot-channel or hot-spot 
condition. Gaussian distributions with known 
means and variances are postulated for each of 
the components and factors entering into the 
computation, and, by means of Eqs. 2 and 3 and 
the propagation-of-errors formula, the estimate 
and the associated standard deviation in Ty, 
(o7,,) are determined. Inferences can then be 
made relative to the true hot-spot wall tempera- 
ture exceeding any given value. 

A typical ETR calculation serves as an exam- 
ple of the procedure. Using the data of refer- 
ence 1, risks of exceeding various temperatures 
are determined. For 395° F, the maximum value 
determined in reference 1, the risk was found 
to be 0.02. 

The procedure has application to any calcula- 
tion of this same general type. 

Another paper on the statistical analysis 
technique is presented in reference 5. The 
section on the analysis of individual factors 
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points out that the engineering hot-channel fac- 
tors arise from distributions which are often 
normal distributions of mechanical data; and, 
if such is the case, the hot-channel factors can 
be calculated by the use of tolerance limits for 
a specified confidence level, the number of 
observations or measurements of the particular 
parameter of a fuel element being considered. 
The fuel-density hot-channel factor is illustrated 
and defined as 


p+p 





Fdensity = 


where pis nominal density and pis variation 
from nominal. 

Table IIl-1 shows the probability for F gensity 
being exceeded, with a confidence level of 
99 per cent, assuming 500 observations are 
recorded. The figures in Table II-1 illustrate 
that, as pincreases, the probability of ex- 
ceeding p + p decreases. 


Table II-l1 PROBABILITY FOR FUEL-DENSITY 
HOT-CHANNEL FACTOR BEING EXCEEDED 








F density Probability of being exceeded 
1.0022 0.10 

1.0027 0.05 

1.0038 0.01 

1.0051 0.001 





Other hot-channel factors can be calculated 
in much the same way, and the problem be- 
comes one of how to combine the various fac- 
tors in a meaningful manner. Consideration 
must be given as to whether any of the factors 
are related, e.g., eccentricity and fuel meat 
thickness. Reference 5 illustrates both a graphi- 
cal and a numerical technique for combining 
two independent factors and briefly discusses 
how to combine related factors, and the com- 
bination of more than two factors, some of 
which are related. The results of the multiplica- 
tion of several independent factors yield a 
table similar to Table II-1. The “probability 
of being exceeded’’ column may then be multi- 
plied by the number of coolant channels in the 
reactor to obtain the expected number of hot 
channels, that is, the number of channels whose 
temperature exceeds that specified by the over- 
all hot-channel factor. Presumably, for a con- 
servative design, one would make this number 
considerably less than unity. It is to be re- 


membered, of course, that hot channels and 
hot-spot factors are used only to take account 
of those deviations whose occurrence and loca- 
tion in the reactor cannot be predicted. Any 
predictable deviation from the average thermal 
situation, such as the deviations due to pre- 
dictable neutron-flux distributions or to pref- 
erential orificing of coolant flow, must be ac- 
counted for in an appropriate manner. The 
approximate result will be, for most reactors, 
that the probability of being exceeded will not 
be multiplied by the total number of coolant 
channels but by some considerably smaller 
number made up of those channels computed 
to be the “hottest’’ on the basis of the pre- 
dictable quantities. An alternative, andperhaps 
preferable, approach would be to decide on 
some acceptable probability of exceeding the 
design condition somewhere in the reactor. The 
ETR, for example, apparently has been suc- 
cessfully operated with about a 2 per cent 
probability of the wall temperature exceeding 
395°F. The corresponding probability for the 
wall temperature exceeding saturation is about 
10 per cent. 


Geometrical Complications 


Reference 6 contains recent information on 
the analysis of turbulent flow and heat transfer 
in noncircular passages. The introduction con- 
tains several references to papers dealing with 
heat transfer and flow wherein turbulent axial 
flow between rods and in eccentric annuli was 
studied. The particular geometries studied were 
a triangular duct and square duct, for Reynolds 
number 10‘ to 10°. The analytical investigation 
proceeded by integrating the differential equa- 
tions for shear stress and heat transfer for 
fully developed flow, using appropriate functions 
for the momentum eddy diffusivity and assum- 
ing the heat-transfer eddy diffusivity to be 
equal to the momentum diffusivity. The results 
were compared with a limited number of experi- 
mental data and were found tobe “in reasonably 
good agreement.’’ The following results are 
quoted from the reference: 

1. The velocities and shear stresses in the re- 
gion near the corner were lower than the average 
values and went to zero at the corner. 

2. When the passages were heated, the heat 
transfer to the fluid in the region near the corner 
was lower than the average value and went to zero 
at the corner. 


_——-— —__ - 


er Ww VS —~ 


oo -/ 


~ we ‘YP We —_—_ 


— . 





HEAT TRANSFER AND FLUID FLOW 9 


8. The friction factors for the noncircular pas- 
sages were somewhat lower than those for a circu- 
lar tube. 

4. When uniform heat sources in the passage 
wall and uniform heat transfer at the surface were 
assumed to occur, the difference between the maxi- 
mum and average wall temperatures was directly 
proportional to the heat flux. 

5. The average Nusselt numbers for the noncir- 
cular passages were somewhat lower than those 
for a circular tube. The average Nusselt number 
was also found to be a function of wall temperature 
distribution and of Prandtl number. 


It was pointed out in the reference that, if 
secondary flows were present in the corner 
regions, the velocity distribution might be dif- 
ferent. Secondary flows, peripheral conduction, 
and turbulent transport in the fluid were neg- 
lected in the analysis and have not been ex- 
tensively investigated experimentally. 

Reference 7 presents an analytical treatment 
of a fairly common problem in reactor design, 
that of determining the temperature of a UO, 
fuel pellet encased in a cylindrical rod. The 
solution to Fourier’s equation was accomplished 
with some interesting modifications, however. 
The heat source was described by the following 
equation (cylindrical coordinates): 


S(r,0) = C,+C,r cos 6 


The cosine distribution was taken as an ap- 
proximation to the physical fact that part of the 
fuel rod “saw’’ the moderator and the re- 
mainder saw other fuel rods. The net neutron 
current would then be greater on the moderator 
side than on the other side. Because of this 
same geometry, the coolant-gas temperature 
was represented by the following function: 


Teoolant = C3 +C4 cos 8 


Circumferential conduction in the metal wall 
of the tube was included, and two cases were 
studied: (1) the pellet is centered in the tube, 
and (2) the pellet is located off-center in the 
tube. The detailed numerical results of the 
study are restricted to the fuel-element sizes 
and specific powers used in the calculations and 
will not be repeated here. In general, however, 
circumferential temperature differences of the 
clad wall and pellet surface (i.e., maximum 
wall temperature—minimum wall temperature) 
up to several hundred degrees were calculated. 
The equations given in the report should enable 


a hot-channel factor for this effect to be cal- 
culated with reasonably good assurance. 


Burnout 


The September 1959 issue of Power Reactor 
Technology contained a section presenting vari- 
ous correlations for burnout heat flux. Tothese 
correlations can be added several others re- 
cently proposed by Bell.’ A statistical analysis 
was performed on available burnout heat-flux 
data from various sources for water flowing 
vertically upward in uniformly heated rectan- 
gular channels at 2000 psia. The range of 
variables studied are: burnout enthalpy from 
540 to 1000 Btu/lb and mass velocity from 
0.2 x 10° to 5 x 10° Ib/(hr)(sq ft). Two correla- 
tions were determined, one for “best fit’’ of 
all the data and one for best fit in the quality 
region (672 = H;, = 1000 Btu/Ib) only. The form 
of the ccrrelation equations is 


; 1(G/10°)* 
$0 _ a(p +) (Bs —%) 
10° 10° 10° 


where B09 = burnout heat flux, Btu/(hr)(sq ft) 
a,b,c,d,e = correlation constants 
G = mass velocity, lb/(hr)(sq ft) 
Hg = enthalpy of saturated vapor, Btu/ 
lb 
H,, = bulk enthalpy at burnout location, 
Btu/lb 


The values for the correlation constants for the 
two correlations are presented in Table II-2, 
along with other pertinent data. 

A frequency distribution graph for the two 
correlations showed that nearly all the data 
fell within +20 from the mean, so a conserva- 


Table II-2 RESULTS OF STATISTICAL ANALYSIS 
0.2 < G/10® < 5 lb/(hr)(sq ft) 





Correlation 1* Correlation 2t 





Correlation 


constants: 

a 0.427 1.37 

b 2.59 1.19 

c 1.81 1.44 

d Pry 1.67 

e 0.406 0.474 
No. of data points 238 186 
Log mean standard 

deviation 0.100 0.0947 





*540 < H, <= 1000 Btu/lb. 
+672 < H, <= 1000 Btu/lb. 
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tive estimate of ¢B0 for design purposes would 
be 65 per cent of the value obtained from 
either calculation. A somewhat disquieting note 
was added, however, when a comparison of the 
correlations was made with 25 runs with a 
rectangular channel having an axial heat-flux 
distribution in the form of a chopped cosine 
(Bettis data). It was shown that the burnout 
data (presumably the burnout heat flux was 
taken as the heat flux at the point of burnout) 
for the cosine shape fell an average of 30 per 
cent below po calculated by the correlations, 
and four of the points fell below the correlation 
by factors ranging from about 2 to 2.9. The 
chopped cosine used hada maximum-to-average 
ratio of 1.38 and a maximum-to-minimum ratio 
of 4.0. In all cases burnout occurred between 
the halfway point, where the heat flux was 
highest, and the endofthe channel. This chopped 
cosine peak-to-average ratio is not particu- 
larly high; the Dresden reactor, for example, 
has an axial peak-to-average ratio of about 2. 
The detailed conclusions of Bell are quoted 
below: 


Based upon this study, the following conclusions 
can be made for vertical upflow of water in rectan- 
gular channels at 2000 psia: 

1, The burnout heat flux is primarily dependent 
on the fluid enthalpy at the burnout point with a 
varying dependency of the fluid mass velocity. 

2. The burnout flux decreases as the mass ve- 
locity increases, for a constant burnout enthalpy in 
the quality range. This inverse effect occurs for 
both uniform and cosine shaped axial heat-flux dis- 
tributions. The higher the quality, the stronger is 
this inverse effect. The opposite is true for sub- 
cooled burnout, with the crossover point occurring 
near the saturated liquid point where mass velocity 
has a minimum effect on the burnout heat flux. 

3. Channel length and length-to-thickness ratio 
(L/S) appear to have no effect on the burnout heat 
flux. Since the range of channel thickness was nar- 
row, no conclusions regarding the effect of this 
variable alone can be made. 

4. For the same conditions at the burnout point, 
the burnout heat flux for a channel having a cosine- 
shaped axial heat-flux distribution is significantly 
lower than for a uniformly heated channel. One 
must take this effect into consideration, then, in 
applying any correlations based upon uniform- 
heating burnout data to reactor design. 


Short Notes 


Present boiling-water reactor designs are 
primarily fuel temperature limited because of 


the economics of fuel-element design. As fuel- 
element technology improves, however, the 
limiting design condition could be burnout. The 
September 1959 issue of Power Reactor Tech- 
nology presents several empirical correlations 
of the critical heat flux in water at high pres- 
sure. Reference 9, recently published, presents 
an extensive analytical approach to the subject 
of pool boiling. The table of contents is quoted 
below: 


Chapter I. A Review of Nucleate Boiling 

II. The Problem of Bubble Growth 

III. Hydrodynamic Aspects of Nucleate 
Boiling 

IV. Hydrodynamic Aspects of Transitional 
Boiling 

V. The Minimum Heat Flux in Transi- 
tional Boiling from a Horizontal 
Surface 

VI. The Critical Heat Flux Density from 
a Horizontal Surface 


Since the condition of pool boiling from a hori- 
zontal surface occurs infrequently in reactor 
problems,-. the work is of general, rather than 
specific, interest. 

Compilations of physical properties of reac- 
tor materials are of continuinginterest. Refer- 
ences 10 and 11 are such publications. The 
former reference contains papers presented at 
the Symposium on Thermal Properties held in 
February 1959, under the sponsorship of the 
American Society of Mechanical Engineers. 
The latter reference investigates thermody- 
namic and transport properties of water and 
water vapor in the critical region. 
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CONTAINMENT 





-— 
In the design of secondary containment shells 


for water-cooled power reactors, current prac- 
tice* is to provide sufficient strength in the 
shell to withstand the pressure rise due to the 
escape of steam and water through the maximum 
credible rupture in the primary coolant system. 
Although credit may be taken for some conden- 
sation of steam on walls and other surfaces if 
the magnitude of the effect can be predicted, no 
deliberate attempt has been made to increase 
the condensation effect in previous containment- 
building designs. 

In recent months, however, a new concept of 
containment has been under investigation’ and 
has been considered for use on the Humboldt 
Bay reactor of the Pacific Gas and Electric 
Company. This concept can best be described 
by the schematic drawing shown in Fig. 2. The 
reactor is suspended in an air-filled cylindrical 
tank called the dry well. The dry well is con- 
nected to a pool of water by vent pipes. Over 
the pool of water is a vaportight container. 
Should a leak or rupture occur in the reactor 
vessel, the steam-water mixture would flow out 
of the vessel and into the dry well, displacing 
the air in the vent pipes. As pressure increased 
in the dry well, the air-steam-water mixture 
would flow into the pool of water wherein the 
steam would condense. If this condensation 
were rapid and complete, the pressure buildup 
within the vaportight enclosure would be small, 
and the water in the pool would act as a heat 
sink for the steam-water mixture. Ideally, the 
water would also scrub solid fission products 
from the vent pipe effluent, should any be pres- 
ent due to the nature of the accident. 

Reference 1 describes the development pro- 
gram conducted to investigate the containment 
system. After analytical investigation of the 
concept, two test facilities were constructed. 
The first facility, the condensing test facility, 
was designed to investigate the behavior of the 
vent units. The facility consisted of water en- 





*A general review of reactor containment was given 
in the June 1959 issue of this journal (Vol. 2, No. 3). 
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closed in a test tank; steam could be injected 
into the water through various experimental 
vents at rates up to about 100,000 lb/hr. The 
first series of experiments was done with atest 
tank 20 ft in diameter by 24 fthigh. The follow- 
ing range of parameters was tested: 


Vents 4-, 6-, 8-, and 14-in.- 
diameter single straight 
pipes and a set of three 
parallel 4-in.-diameter 
pipes 

Depth of 1 in. to 6 ft 

submergence 

Direction of Vertically downward and 

discharge horizontal 

Steam flow 10,000 to 93,000 lb/hr, 
100 psig steam at near 
saturation 

Tank water 50 to 150°F 

temperature 


In all but a few cases, the steam was completely 
condensed; the few exceptions were for the vent 
pipe discharge within a few inches of the sur- 
face of the water. Difficulty with tank vibration 
was experienced when the water pool tempera- 
ture was hotter than 120 to 130°F, but this 
should not be a serious restriction in an actual 
plant design. A second series of tests was done 
to investigate the vent discharge when a small 
volume of water was used. Into the large tank 
was placed a small compartment containing 
water; the compartment volume was varied but 
was under 100 cu ft. Ninety-one tests were run 
to investigate the effect of steam flow rate, vent 
diameter, and compartment geometry. Com- 
plete condensation was obtained, and mixing of 
the water in the compartment was said to be 
excellent. 

The second test facility utilized was the 
transient test facility. Basically, the facility is 
a model of the arrangement shown in Fig. 2. 
The volumetric size was '/9) that of the Hum- 
boldt Bay (California) Power Plant reactor. 
Discharge of the “reactor” water to the.dry well 
was accomplished by breaking a rupture disk by 
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Figure 2—Simplified schematic of containment sys- 
tem.° 


means of a blasting cap. The test parameters 
were varied as follows: 


Reactor pressure before 


firing, psi 600, 1000 
Fraction of reactor vessel 

filled with liquid before 

firing 0.66, 0.80 
Diameter of reactor dis- 

charge orifice, in. 0:6. 1.2,. 1:6 
Dry-well volume, cu ft 17.5, 20.3, 26.0 
Vent area, sq ft 0.5, 0.9, 2.1 
Vent depth of submer- 

gence, ft 0.75, 1.0, 1.5 


One of the purposes of the transient tests was 
to test the analytical calculations on peak dry- 
well pressures. Although some of the assump- 
tions used in the derivations of the analytical 
equations are given in reference 1, considerably 
more information is contained in reference 2. 
Since the equations are of particular interest to 
the specialist, they will not be repeated here; in 
general, the experimental results reported indi- 


cated that the analytical method of predicting 
peak dry-well pressure is satisfactory and con- 
servative. The following variables are listed in 
order of decreasing importance in determining 
the maximum dry-well pressure: orifice size, 
depth of submergence, dry-well volume, and 
vent area. Of these variables the most important 
one, orifice size, is, possibly, the variable over 
which the designer has the least control, since 
it is the size of the leak, hole, or rupture that 
occurs in the reactor primary system. Thus the 
design of the containment system hinges on the 
estimate of the maximum credible size of rup- 
ture. Figure 3 summarizes the calculations for 
the Humboldt Bay reactor. The maximum credi- 
ble operating accident postulated for that reac- 
tor assumes the following:’ 

1. Near-instantaneous main steam or feed- 
water pipe break. 

2. Broken parts or other failed parts become 
missiles. 

3. Blast effects and reactor-vessel reaction 
occur from the escaping steam and water. 

4. Imperfections in the containment permit 
some leakage. 
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Figure 3—The effect of peak dry-well pressure on 
vent area. 


On the other hand, a rupture of the reactor 
vessel, a large nuclear excursion, and an ex- 
plosive chemical reaction are deemed incred- 
ible. 

The feed-water line to the Humboldt Bay re- 
actor is an 8-in. line. The main steam line is a 
12-in. pipe; the calculations presented in Fig. 3 
assume that water, rather than a steam-water 
mixture, is flowing through the break. The upper 
left-hand portion of the curve, a horizontal line, 
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results when the vent area is so small that the 
peak dry-well pressure is attained before a 
significant amount of steam can escape through 
the vents into the pool water; for this condition 
the dry-well pressure depends only on the 
volume of the dry well. The horizontal line at 
the lower right-hand end of the curve represents 
the condition in which the vent area is large 
enough that it contributes little to the resistance 
to flow, and the dry-well pressure depends on 
the break size and the depth of submergence of 
the vent. The intervening portion of the curve 
is an area wherein the peak dry-well pressure 
depends primarily on the flow resistance of the 
vents, and the effect of dry-well volume and 
vent submergence is minimized. It is the region 
in which the Humboldt Bay containment will 
“operate.” If steam condensation in the dry 
well should tend to reduce the pressure below 
atmospheric, vacuum breakers are provided to 
equalize the pressure between the dry well and 
the pool air space. 

Simulated fission products were added to the 
model reactor vessel for some of the transient 
test facility operations. The specific objective 
of these tests was to determine the effectiveness 
of the water pool in retaining fission products 
introduced through the vents. The results were 
said to indicate that, in general, the water pool 
is an “effective barrier.” The results of the 
tests are shown in Table III-1. 


Table III-1 RESULTS OF TESTS ON ENTRAINMENT 
OF SIMULATED FISSION PRODUCTS IN WATER! 








Per cent 
Initial amount released to Concentration 
in pressure containment in containment 
Element vessel volume volume, mole % 
—50% —50% 
Krypton 250 cc (std.) 25 +100% 0.003 + 100% 
Xenon 250 ce (std.) Less than 45 Less than 0.005 
Iodine 1 lb 5.8 x 1075 1.83 x10 
7.2x10 at 2.3107 at 
5 hr 5 hr 
Sodium 100g ia x10 Lexi" 
iodide 
Zinc 110g 2i-xiio 3.7 xo” 
sulfide 





Note: containment volume, after firing, 78.5 cu ft (air); pool 
volume, 129 cu ft; pressure-vessel volume, 3.12 cu ft. 


This type of containment does have some 
complicating features. For a bottom-driven 
control-rod system, the rods are located within 


the dry well. Hence the Humboldt Bay plant pro- 
vides for access to the bottom of thedry well in 
order to maintain the rod drives. Further, all 
primary system pipes that could conceivably 
fail must be either contained within the dry well 
or provided with fast-acting valves to prevent 
discharge of the reactor water in the event of a 
rupture outside the dry well. 

Experiments of a rather similar nature have 
been performed by Sargent & Lundy. In tests 
recently reported,° a pressure drum, simulating 
a water-cooled reactor, was housed inside a 
vertical containment shell. The pressure drum 
was 42 in. in diameter by 23 ft in length, and 
the containment tank was 14 ft indiameter by 32 
ft high. The drum had a number of 12-in. open- 
ings along its side as well as one at the bottom, 
all of which could discharge downward toward a 
pool of cold water covering the bottom of the 
containment tank. A 600-psi rupture disk was 
placed over one of the openings, and the drum 
was filled with a given amount of water and sub- 
jected to heat; the pressure in the containment 
tank was measured after rupture of the disk. 
For most of the tests, 8000 lb of hot water was 
contained in the pressure drum, and 17,000 lb 
of cold water was held in the pool at the bottom 
of the containment tank. As the height of the dis- 
charge nozzle was varied, in successive tests, 
from a few inches above the cold-water surface 
to about 20 ft above the surface, the maximum 
pressure recorded in the containment tank in- 
creased from about 5 to about 50 psig. The 
theoretical pressure rise, if there had been no 
cooling of the discharged water other than by 
mixing with the gas in the containment tank, 
was 153 psig. When 4800 lb of water was dis- 
charged with no cold water in the containment 
tank, a maximum pressure of 95 psig was Oob- 
served, whereas the calculated “no cooling” 
maximum was 108 psig. In this test only, shock 
wave formation was reported. 
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FUEL CYCLES 





Reactivity Lifetimes with 
Plutonium and U** Enrichments 


Under the plutonium recycle program at Han- 
ford,* a study has been,made of the relative 
reactivity worths of U** and plutonium (ofvari- 
ous isotopic compositions) as feed for a uni- 
formly graded irradiation fuel cycle.! The uni- 
formly graded cycle is one in which the reactor 
fuel composition consists of a uniform mixture 


- of fuels having all degrees of radiation exposure, 


from zero (feed fuel) to the maximum (discharge 
fuel). Strictly speaking, the uniformly graded 
cycle can be achieved in an actual reactor only 
by discharging and reloading fuel ininfinitesimal 
amounts over the entire volume of the reactor. 
Practically, the cycle can be approximated by 
partial reloading schemes which reload perhaps 
20 per cent of the reactor fuel at atime. (A 
discussion of fuel reloading programs was given 
in the section on Refueling Methods for Solid 
Fuel Reactors in the September 1959 issue of 
Power Reactor Technology, Vol. 2, No. 4.) 

The study considered a single lattice but 
varied the neutron utilization by varying the 
fractional neutron leakage. The fractional leak- 
age of neutrons from a reactor may be specified 
by giving the value of k.. required for criticality: 
higher values of k,. represent higher fractional 
neutron leakage. In the study it was convenient 
to vary k./e rather than k.. The fast-fission 
factor, €, is not far from unity. In some re- 
spects the loss of neutrons by leakage may be 
considered as equivalent to loss by any other 
process. However, other losses will not vary in 
the same way with fuel composition. Other 
things being equal, only that portion of the leak- 
age which is made up of thermal neutrons will 
vary importantly with fuel composition. It will 
vary inversely with the effective thermal- 
neutron-absorption cross section of the reactor 
(2a). Neutron losses due to parasitic absorption 





*For other activities in the plutonium recycle pro- 
gram, see the discussion of the Plutonium Recycle 
Reactor (Sec. IX) in this issue. 


by such materials as fission products and struc- 
tural materials will vary inversely as the mac- 
roscopic absorption cross section of the fis- 
sionable isotope. However, since the study 
considers arbitrary values ofthe fractional neu- 
tron loss, these distinctions are not of great 
importance. 

The lattice used for the investigation is char- 
acterized by a constant value of the resonance 
slowing-down ratio 


Vi 


Spent ae 
Vm (EZs) 


where V;/V,, is the fuel-to-moderator volume 
ratio, >, is the slowing-down power of the 
moderator in the resonance region, and F isthe 
fuel disadvantage factor for resonance neutrons. 
This amounts to holding the U?** resonance- 
escape probability constant except for the small 
variation due to the slight change in u?*8 content 
as enrichment is varied. The lattice resonance- 
escape probability (p) is defined in the study as 
the ratio of the number of neutrons reaching 
thermal energy to the number of neutrons slow- 
ing down past the U?** fast-fission threshold in 
an infinite lattice. In this definition of p, the net 
neutron production by resonance fission is sub- 
tracted from the loss of neutrons by absorption 
in U’*® and other resonance absorbers. The 
values of p ranged from 0.927 to 0.949 for the 
Pu? cases considered and from 0.898 to 0.924 
for the U**> cases. The values for mixtures of 
plutonium isotopes were lower (down to 0.824) 
because of resonance absorption in the higher 
plutonium isotopes. The value of 3.00 for the 
resonance slowing-down ratio was stated to cor- 
respond to a moderator-to-fuel volume ratio of 
about 1.0 for H,O, 13 for D,O, and 55 for graph- 
ite, when the fuel is UO, of density 10 g/cm’. 
However, the value chosen for the thermal 
utilization eliminates H,O as a moderator to 
which the results might apply. The resonance 
slowing-down ratio also specifies, along with 
the fuel composition, the neutron spectrum in 
the thermal and near-thermal ranges. In the 
study the formulation of Westcott? was followed 
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Table IV-1 RELATIVE PERFORMANCES OF Pu’® AND U**> AS FISSIONABLE FEED IN REACTORS 
HAVING DIFFERENT EFFICIENCIES OF NEUTRON UTILIZATION! 








Percentage of 


Fissions per Mwd per net gram 








Reeaeneeet, Reactivity life, fissionable atom of fissionable total — 
cheese Mwd/ton in feed isotope destroyed ye AR 
able isotope - = U5 is the 

k../e iniead Py239« U235* — y235*/py239— py239 235235 /py239— py239 235235 ppy238 feed isotope 
1.05 0.71 6,896 4,927 0.715 1.145 0.805 0.703 2.18 2.67 1.225 63.4 
1.00 12,764 11,809 0.925 1.505 1.37 0.910 2.22 2.36 1.063 §2.2 
1.50 21,895 20,947 0.958 1.721 1.62 0.955 Ege 1.005 49.2 
1.15 1.00 5,251 4,180 0.796 0.872 0.683 0.734 1.05 1.30 1.238 71.6 
1.50 14,697 14,197 0.966 1.155 1.098 0.950 1.52 1.68 1.105 63.2 
2.00 21,496 20,958 0.976 1.267 1.235 0.975 1.51 1.57 1.040 61.5 
1.30 1.50 7,163 6,439 0.620 0.563 0.498 0.884 1.05 1.36 1.295 81.7 
2.00 11,640 11,832 1.015 0.686 0.739 1.077 1.05 1.21 1.152 76.9 

1 1 


3.00 21,222 24,153 1.138 0.834 0.934 1.120 


*The top column subheading Pu’ 


23 1.170 74.6 





signifies that the fissionable material in the feed is pure Pu’; the subheading u> 


signifies that the fissionable material in the feed is pure 235, the subheading u235 /py238 signifies the ratio of the numbers 


in the two preceding columns. 


in accounting for spectrum effects. The usual 
effective cross sections (6) were used 


6 = o (g+ rs) 


where 0, is the 2200 m/sec cross section; g and 
S are terms characteristic of the substance, 
which define the departure of the cross section 
from the 1/v law in the “thermal” and “reso- 
nance” regions, respectively; and yr is the frac- 
tion of epithermal neutrons in the spectrum. In 
the study the values of r ranged from 0.0650 to 
0.17 for the Pu’**-fed cases and from 0.0475 to 
0.1010 for the U**°-fed cases. The moderator 
temperature was in all cases 200°C. 

The value of the thermal utilization (/) was 
arbitrarily set at 0.9 for natural uranium, and 
the average macroscopic parasitic absorption 
cross section was held constant; thermal utiliza- 
tion was calculated for other fuel compositions 
by taking into account both the change in fuel 
absorption and a typical change for the disad- 
vantage factor. At 3 per cent enrichment the 
thermal utilization for the feed fuel had in- 
creased to 0.954 for the U’® case and to 0.976 
for the Pu’*® case. Xenon and samarium losses 
are included among the neutron losses repre- 
sented by k,,/€; other fission products are 
treated by a burnup equation. 

Three values of k./e were considered, and 
three enrichments of the feed fuel were con- 
sidered for each value of k../e. The results for 
the cases in which the feed is pure Pu’*® and 
pure U’*> are compared in Table IV-1. As would 
be expected, the reactivity life of the fuel in- 
creases rapidly with feed enrichment at a given 


value of neutron loss, and the reactivity lifetime 
at a given enrichment decreases rapidly as the 
neutron losses are increased. When the neutron 
loss is small or when the enrichment at a given 
neutron loss is sufficient only for a rather short 
reactivity life, the plutonium feed gives con- 
siderably longer reactivity life than does the 
u’*> feed. As the enrichment is increased to 
give longer reactivity life, the reactivity life 
given by the U’** feed becomes relatively better 
and, in the highest neutron loss case, surpasses 
that for the Pu’®® feed. The same trends are, of 
course, evident if the reactivity life of the fuel 
is expressed as the number of fissions occurring 
over the lifetime of the fuel per fissionable 
atom initially present in the fuel. 

The rather good performance of Pu*°* so far 
as reactivity lifetime is concerned is, at first 
glance, surprising, since the value of 7 (fission 
neutrons produced per neutron absorbed) for 
Pu?*® is less than that for U’*> The behavior of 
the Pu’® would be easier to discuss if values for 
the cross sections and neutron yields of the ura- 
nium and plutonium isotopes had been tabulated 
in the report.* However, the trend ofthe results 
does seem reasonable. 

Two types of effects must be recognized. 
First, the neutron balance for the U?* feed 
cases is not that which would result if all fis- 
sions were in U* since, over the life of the 
fuel, a substantial fraction of the fission occurs 
in plutonium isotopes. The last column of Table 
IV-1 shows, for the U’**-fed case, the percentage 


239 





*The values used were those of reference 2. 
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of total fissions occurring in U***. One would 
expect the neutron balance behavior to be a 
weighted average between that characteristic of 
the plutonium 7 and that characteristic of the 
u?*5 » and to approach the plutonium case as the 
fraction of fissions in U*’ decreases. 

The second effect is that Pu’®’, although it has 
a lower 7n value, has a higher effective fission 
cross section. The higher fission cross section 
reduces the neutron losses represented by (1 — 
f);it also yields a higher reactivity for the feed 
fuel. One would expect that these effects would 
be most pronounced at low enrichments that 
yield relatively short fuel lifetimes. Presuma- 
bly when one tries to attain long reactivity life- 
time, particularly in reactors with high neutron 
losses, the effect of the lower plutonium 7 on 
conversion ratio is overriding, and the fuel life- 
time relative to that obtainable with U’* de- 
creases. 

Further light is shed on the relation between 
the Pu?’ and the U?* feeds ifthe heat output per 
net gram of fissionable isotope destroyed is 
considered. This quantity can be determined 
from the isotopic compositions of the feed and 
discharge fuels as given in the reference re- 
port.! The quantity is closely related to the 
average conversion ratio of the reactor and is, 
of course, a quantity of primary importance for 
the effective utilization of nuclear fuel. It may 
be seen from Table IV-1 that this specific pro- 
duction is in all cases higher for the U** feed 
and is highest, relative to the value for Pu? 
feed, in those cases when most of the fissions 
occur in U*>. It is evident from these numbers 
that the U*°-fed reactors in all cases have 
higher conversion ratios than the Pu238 cases, 
and one would infer that the longer reactivity 
lifetimes exhibited by the plutonium-fed reac- 
torsare due to the higher reactivity of feed fuel 
at a given enrichment. Thus the U">_fed reac- 
tors, in producing a given quantity of energy, 
burn less net fissionable isotope (U?*> plus Pu2*® 
plus Pu”*!), but the fuel must be discharged, be- 
cause of reactivity deficiency, at a relatively 
higher fissionable-isotope content. 

The reference stresses the comparison of 
feed fuel composed of three additional mixtures 
of the plutonium isotopes. These mixturesare 
defined as mixtures which would be obtained if 
pure Pu”® (no U"** present) were exposed, in a 
reactor having aconstant epithermal ratio (rv) of 
0.08, until 50 percent, then 75 per cent, and then 
90 per cent of the original Pu?*® was burned. The 


Table IV-2 ISOTOPIC COMPOSITION OF 
STOCK PLUTONIUM! 
Per cent burnout Isotopic composition, wt 


of Pu?®® atoms Pu??? Pu?4? Pu’! Pu2# 
0 100.00 0 0 0 
50 75.20 19.50 4.92 0.38 
75 52.95 30.12 14.28 2.65 


Table IV-3 COMPARISON OF U*5 AND PLUTONIUM 
ISOTOPE MIXTURES! 


Maximum attainable 


iahite exposure 
ment, Fissions 
% fission- per fission- 
able Type of able atom 
Reo/€ isotope enrichment in feed Mwd/ton 
1.05 0.71 ans 0.81 4,930 
Pu239 1.15 6,900 
Pu, 50% burned 1.15 6,950 
Pu, 75% burned 122 7,300 
Pu, 90% burned 3.32 8,130 
1.30 3.0 u 0.93 24,150 
Pu239 0.83 21,200 
Pu, 50% burned 0.69 17,600 
Pu, 75% burned 0.63 15,900 
Pu, 90% burned 0.57 14,500 


compositions of the various plutonium mixtures 
are given in Table IV-2. Selected results of the 
reactivity lifetime calculations are listed in 
Table IV-3. In Table IV-3 the enrichment refers 
to the percentage of fissionable isotope (Pu’®® 
plus Pu’*!) in the fed fuel. It is to be noted that 
the “burned” plutonium mixtures show accentu- 
ated trends, relative to. U’**, of the same type 
shown by Pu?*®. This situation is a more com- 
plex one. Reference to Table IV-3 indicates that 
the number of fissions obtained per fissionable 
atom in the feed increases with the more highly 
burned plutonium in reactors with low neutron 
losses. Also, although not specifically calcu- 
lated, the comparison with U** in Table IV-1 
based on megawatt-days per net gram of fis- 
sionable material destroyed should improve for 
highly burned plutonium since the 7, or neutrons 
per absorption, for Pu! exceeds that for Pu*. 
Important factors contributing to the behavior 
are, no doubt, the high effective fission cross 
section of Pu™!, particularly in the low- 
enrichment case, and the important resonance 
absorption by Pe”, particularly in the higher 
enrichment cases. 
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Self-sustaining Plutonium Recycle 
in Sodium Graphite Reactors 


Reference 3 is a study of the useof plutonium 
recycle with natural-uranium feed in sodium 
graphite reactors. Existing designs of the so- 
dium graphite reactor require enrichments 
which are too high for this kind of operation, 
but it is considered a possibility provided a 
calandria type design is used and low-cross- 
section materials are employed for fuel jackets 
and structural parts. 


If the plutonium manufactured during fuel ir- 
radiation in a reactor is recycled as part of the 
feed for the reactor, it will supply some of the 
fissionable isotope needed for feed, and the 
proportion of U*®> in the feed may be reduced. 
In order for the U?** feed requirement to be as 
low as the U?* content of natural uranium, the 
reactor must make rather effective use of the 
fission neutrons. This required high neutron 
economy may be manifest as a relatively low 
enrichment requirement (U?** plus fissionable 
plutonium isotopes) for operation, or as a high 
conversion ratio. Both of these characteristics 
are favored if the parasitic loss of neutrons is 
kept low and if the number of neutrons produced 
by the fissionable isotope per neutron it absorbs 
(n) is high. The parasitic neutron losses depend 
primarily on the reactor composition, whereas 
n, at least for Pu’®®, depends rather strongly 
upon the neutron temperature in the reactor. 

The study reported in the reference, which is 
an approximate study, approaches the problem 
by calculating the rate of plutonium production 
and the isotopic composition of the plutonium in 
reactors having various conversion ratios, after 
the recycle program has reached an equilibrium 
state. These same reactors, if fueled exclu- 
sively with U** as the fissionable isotope, would 
exhibit somewhat different initial conversion 
ratios and would require, for operation, enrich- 
ments which may be specified. The relation be- 
tween the average conversion ratio of the re- 
cycled reactor and the initial conversion ratio 
of the U’*5-fed reactor involves only the 7 of 
Pu?*® if it is assumed that 7 for U*® and 7 for 
Pu"! are independent of neutron spectrum, as 
is done in the study. If one then choosesa value 
of n for Pu’*®, one can specify a set of U***-fed 
reactors, each characterized by a value of the 
initial conversion ratio and an enrichment re- 
quirement for operation, which would just op- 


erate with natural-uranium feed and plutonium 
recycle. Sets of reactors of this kindare speci- 
fied, for various values of 7 of Pu’*®, in Fig. 4, 
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Figure 4—Conditions for self-sustaining plutonium 
recycle (equilibrium plutonium concentration corre- 
sponds to zero irradiation per cycle). 


which is reproduced from the reference. For 
example, if the neutron temperature in the re- 
actor is consistent with a value of 1.90 for 7 of 
Pu?*® and if the reactor requires an enrichment 
of 1.1 per cent when operating on pure U’* and 
u?** fuel, and yields an initial conversion ratio 
under those conditions of about 0.7, then it will 
just operate with natural-uranium feed and 
equilibrium plutonium recycle. The attainable 
fuel exposure under these conditions would be 
infinitesimal. To obtain longer exposure at the 
7 value of 1.90, one must move to the left of the 
1.90 curve, choosing a reactor which, with U’* 
feed, either requires lower enrichment or yields 
a higher conversion ratio, or both. The farther 
to the left one moves, the higher will be the at- 
tainable fuel exposure. 

In Fig. 5, curves of Maxwell-average 7 values 
for Pu’*® as functions of neutron temperature 
are reproduced, as derived in the reference 
from three different sources. It is evident that 
the value of 7 is strongly affected by the 
thermal-neutron temperature; and if values 
from Fig. 5 are related to the curves of Fig. 4, 
it is also apparent that these changes represent 
quite important changes in the enrichment and 
conversion ratio requirements for self- 
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sustaining plutonium recycle. Since the neutron 
temperature will be at least as high as the 
moderator temperature, it appears probable that 
a sodium graphite reactor designed to operate 
on self-sustaining plutonium recycle must be 
one that can be made critical with uranium fuel 
of low enrichment, probably below 1 per cent 
y235, 

The basic design change, which is proposed 
for improving the neutron economy of the so- 
dium graphite reactor, is from a core using in- 
dividually canned hexagonal blocks of graphite to 
a calandria in which there is, in effect, a single 
block of graphite pierced by holes for the fuel 
elements. Each fuel element is located within a 
process tube through which the sodium flows, 
and the process tubes lie within the calandria 
tubes which constitute a barrier between the 
sodium and the graphite in the event of a 
process-tube leak. Uranium-metal fuel rods 
are used for maximum reactivity, in 7-rod 
clusters. The fuel jackets, process tubes, and 
calandria tubes are all of zirconium. The graph- 
ite is cooled by conduction to the sodium coolant, 
across the helium-filled gap between the calan- 
dria tube and the process tube. This results in 
a rather high graphite temperature, which is 
estimated as an average of 600°C. 

The general characteristics of a reactor of 
this type are given in Table IV-4. The initial 
conversion ratio of the reactor, when fueled 
with U**>, is 0.70. If the reactor is operated 
with plutonium-recycle and natural-uranium 
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Figure 5— Maxwell-average 7 values for Pu?9, 


Table IV-4 GENERAL SPECIFICATIONS OF AN SGR 
DESIGNED FOR SELF-SUSTAINING 
PLUTONIUM RECYCLE 





Power: 
Thermal power, Mw 630 
Electrical power, Mw 250 
Core: 
Diameter, ft 24 
Height, ft 18 
Lattice spacing, in. 14 
Average reflector thickness, ft 2 
Average moderator temp., °C 600 
Fuel: 
Material Uranium metal 
Enrichment 0.9 at.% v2 
Fuel-element configuration 7-rod cluster 
Rod diameter, in. 0.75 
Number of fuel elements 378 
Fuel specific power, Mw/ton 8.3 
Cladding material Zirconium alloy 
Cladding thickness, in. 0.010 
Process-tube inside diameter, in. 2.79 
Process-tube thickness, in. 0.050 
Primary coolant: 
Sodium inlet temp., °F 625 
Sodium outlet temp., °F 1000 
Sodium pressure drop in core, psi 15 


Heat flux (average), Btu/(hr)(sq ft) 236,000 


feed and with an approximation to uniformly 
graded irradiation in the radial direction, it is 
estimated that the average energy yield of the 
discharged fuel would be about 3750 Mwd/ton. 
At this exposure level, about 30 per cent of the 
initial U**> content of the feed would be burned. 
Thus the self-sustaining plutonium-recycle op- 
eration, in this reactor, although it would avoid 
the necessity for isotopic enrichment of the fuel, 
would be relatively wasteful of U?*. 

The reference considers also the burning of 
recycled plutonium in fuel elements which are 
separate from the uranium fuel elements. It is 
thought that this scheme might yield somewhat 
higher U**> burnup if whatever plutonium dis- 
card occurs during the recycle program could 
be made up of.the most exposed plutonium, thus 
preserving a somewhat higher ratio of fissiona- 
ble isotopes in the recycled material. 

The conclusion of the study is summarized in 
the following quotation: 


It appears technically feasible to design a SGR to 
operate on self-sustaining plutonium recycle. Such 
a reactor, however, would of necessity have a large 
core and low (for an SGR) volumetric power density 
and would require exclusive use of a low-cross- 
section material such as zirconium alloy for clad- 
ding and core structure. The fuel specific energy 
yield and U2** burnup attainable is uncertain but ap- 
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pears to be on the low side, perhaps 3000 to 5000 
Mwd/ton and 25 to 40 per cent U?* burnup. Such a 
system would almost certainly compare unfavorably 
with a slightly enriched SGR if the fuel inventory 
annual rate and the unit costs for U**® enrichment 
and for fuel fabrication and reprocessing, currently 
prevailing in the United States, are applied. Ex- 
pressed in a different way, it appears that an opti- 
mized plutonium-recycle case will require some 
enrichment of the uranium feed. These conclusions 
in large measure reflect the low reactivity of plu- 
tonium as a reactor fuel and, therefore, they apply 
in general to other thermal-reactor systems. 


Presumably the final statement is intended to 
apply only in a qualitative sense, for there are 
certain reactor types, notably some employing 
D,O moderation, which are expected to achieve 
relatively long fuel life with once-through 
natural-uranium feed. These reactors should 
show quite good economic performance with 
self-sustaining plutonium recycle, provided that 
the costs specifically attributable to the recy- 
cling process can be discounted. 


General Studies 


The results of an extensive fuel-cycle study, 
concerned primarily with the development of a 
computer code for calculating fuel cycles and 
fuel-cycle costs, has been reported in refer- 
ence 4. The code, named FUELCYC, has been 
written for the IBM-704 computer. It computes 
the distribution of neutron flux with respect to 
energy and position in the reactor, derives ef- 
fective cross sections for each nuclide at each 
point in the reactor, uses the effective cross 
sections to project the change of nuclide con- 
centration at each point with time, determines 
the conditions under which the reactor is just 
critical, and evaluates fuel-cycle costs. Geo- 
metrically, the code solves the two-dimensional 
(radial and axial) diffusion equations in an 
axially symmetric cylinder. Four major energy 
groups are treated: fast fission, fast, resonance, 
and thermal, with the resonance group further 
divided into four subgroups. Leakage is as- 
sumed to occur in the fast and thermal groups, 
and each fuel nuclide is assumed to absorb 
neutrons in one of the resonance subgroups as 
well as in the thermal group. The neutron- 
energy spectrum in the thermal region (below 
0.45 ev) is computed from the Wilkins equation, 
and effective thermal cross sections are derived 
by averaging over the spectrum. 


Although the code is adaptable to the computa- 
tion of fuel-cycle quantities for many fuel 
charging and discharging schemes, four such 
schemes have been specifically written into the 
code: 


1. Batch irradiation of fuel fixed in place in 
the reactor, with criticality maintained by a 
uniform distribution of control poison. 

2. In-out irradiation, in which there is acon- 
tinuous replacement of fuel elements, the new 
elements being added at the axis of the core, 
partially irradiated elements being shifted radi- 
ally outward, and fully exposed elements being 
removed from the periphery. 

3. Out-in irradiation, in which fuel elements 
are fed continuously, in the reverse direction to 
the in-out scheme. 

4. Graded irradiation, in which fuel elements, 
fixed in place in the reactor, are discharged in- 
dividually and replaced by fresh elements on 
such a schedule that the average composition of 
fuel in each local radial region of the reactor 
remains time independent. 


In all of the continuous reloading schemes (2, 3, 
and 4), the reactor remains critical without the 
addition of control elements or poison. In all 
cases studied, the fuel elements were assumed 
to be of full core length. 

Although the main contribution of the study is 
the development of the FUELCYC code, some of 
the fuel-cycle calculations that have been made 
with the code are of general interest. The re- 
sults of applying the four fuel reloading schemes 
to a pressurized-water reactor were studied in 
some detail. The reactor considered was simi- 
lar in geometry and composition to the Yankee 
Atomic Electric Company reactor. It was, of 
course, necessary to vary the enrichment of the 
fuel to meet the demands of the different fuel 
cycles. The actual reactor is fueled with 51,420 
lb of UO, at an enrichment of 3.4 wt.% U’*. It 
is cooled and moderated by light water at a 
pressure of 2000 psi and at a temperature of 
516°F. The reflector consists of a '/-in. steel 
baffle with an effectively infinite amount of 
water beyond. The fuel elements consist of UO, 
pellets, 0.29 in. in diameter, loaded into 20- 
mil-thick stainless-steel tubes. The tubes are 
placed on a 0.42-in.-square pitch to form a 
cylindrical core having a 3.1-ft radius and 
7.7-ft height. Other pertinent data are listed in 
Table IV-5. 
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Table IV-5 CHARACTERISTICS OF PRESSURIZED- 
WATER REACTOR USED FOR FUEL- 
CYCLE CALCULATIONS! 





Core dimensions, in.:* 


Equivalent diameter 75.68 

Total active length 92.257 
Temp., mean (Tmd), °F 516 
Power, Mw:7 

Thermal 480 

Net electric : 134 


Nuclear properties: 
Initial resonance-escape 
probability for U5 (3) 0.738* 
Initial resonance-escape 
probability for struc- 
tural materials and 


coolant (P-) 0,.942* 
Fast-fission factor (e) 1.0584* 
Fermi age (rT), cm? 51.5* 
Thermal disadvantage 

factor (w) 1.141* 
Reflector savings (6p = 

6x), cm 7.5* 


Microscopic 2200 m/sec 

absorption cross sec- 

tion for H,O, barns 0.575} 
Microscopic transport 

cross section for H,O, 

barns 70t 
Microscopic slowing- 

down power for H,O 

[(Eds) H,o]) barns 41.2} 
Microscopic 2200 m/sec 

absorption cross sec- 

tion for type 348 


stainless steel, barns 3.2068 
Atomic weight of type 
348 stainless steel 55.638 
Volume, Weight, 
Item cu in, lb 


Inventories: * 


UO, (wt.% U3? = 3.4) 141,091 51,420 
Stainless steel (type 348) 48,025 13,520 
H,O 208,330 5,881 
Zirconium 12,137 2,850 
Void 5,421 


Total 415,004 73,671 
*W. H. Arnold, Jr., private communications of Mar. 23 
and Feb. 26, 1959. 
+Yankee Atomic Electric Company Researchand Develop- 
ment Program, Reference Design, Report YAEC-1, 1957. 
tA. M. Weinberg and E. P. Wigner, ‘‘The Physical Theory 
of Neutron Chain Reactors,’’ The University of Chicago 
Press, 1958. 
§D. J. Hughes and R. B. Schwartz, Neutron Cross Sec- 
tions, Second Edition, Report BNL-325, 1958. 


Table IV-6 lists the more important results 
of the fuel-cycle calculations. The independent 
variable (first column) is the maximum local 
burnup of fuel, that is, the burnup atthe position 
of maximum exposure in the most highly ex- 


Table IV-6 PERFORMANCE CHARACTERISTICS OF 
PRESSURIZED LIGHT-WATER REACTOR FOR 
VARIOUS FUEL SCHEDULING METHODS 





Fuel scheduling 
method, max. 
local burnup, 
Mwd/ton Batch In-out Graded Out-in 





At.% U® in Feed 


10,000 3.14 2.93 2.95 3.02 
20,000 3.58 3.17 3.20 3.39 
30,000 4.10 3.42 3.53 3.79 
40,000 4.71 3.70 3.88 4.20 
50,000 5.34 4.06 4.24 

60,000 6.09 


Average Burnup, Mwd/ton 


10,000 4,200 7,900 7,900 7,900 
20,000 10,000 16,300 16,300 16,300 
30,000 16,800 25,100 25,100 25,100 
40,000 24,200 34,400 34,400 34,400 
50,000 32,200 44,000 44,000 

60,000 41,200 


Maximum -to-Average Power Density Ratio 


10,000 1.90* 3.28 2.54 1.80 
20,000 1.49* 4.20 2.39 1.47 
30,000 1.37* 5.16 2.30 1.38 
40,000 1.29* 6.17 2.23 1.35 
50,000 1.25* 7.20 2.18 

60,000 1.23* 





*At end of batch cycle. 


posed fuel element at the time of discharge. 
The three sections of the table show, as func- 
tions of this variable, the required enrichment 
of the feed fuel, the average fuel burnup, and 
the maximum-to-average power density ratioin 
the fuel. It is to be noted that the enrichment 
required for a given average burnup increases 
in the order: in-out, graded, out-in, batch; with 
the differences being relatively small except in 
the case of the batch reloading system. The 
average burnups for the continuous reloading 
schemes are less than the maximum local 
burnup because of the axial variation of burnup 
in a given element. The maximum-to-average 
burnup ratio decreases with increasing maxi- 
mum burnup because of the axial flux flattening 
caused by the preferential burnup near the cen- 
ters of the elements. The maximum-to-average 
burnup ratio for the batch case is considerably 
higher because it is determined by both the 
axial and the radial power distributions in the 
reactor. Again, the ratio decreases with in- 
creasing burnup. The ratio of maximum-to- 
average power density is given, for the batch 
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case, at the end of the batch lifetime. At the 
beginning of the lifetime, it has the value 2.70. 
Both the initial and final ratios are, of course, 
determined by the spatial distribution of control 
poison over the batch lifetime, and the initial 
ratio could be made somewhat lower, at the ex- 
pense of a somewhat higher final value than that 
tabulated, by an initial concentration of poison 
near the center of the core to flatten the power 
distribution. It is apparent that the maximum- 
to-average power-density ratio is quite poor 
for the in-out refueling system and becomes 
worse with longer fuel burnup, whereas the 
ratio for the out-in system is quite good and 
becomes progressively better with increasing 
burnup, at least up to about 35,000 Mwd/ton 
average burnup. It is to be emphasized, of 
course, that these maximum-to-average ratios 
are the smooth ratios for uniform lattices and 
do not include local peaking effects such as 
those due to water channels. 

The power-distribution patterns established 
in the reactor after long burnup of the fuel are 
particularly interesting. In Fig. 6 are repro- 
duced the contours of constant power density — 
in one quadrant of a plane which bisects the core 
and contains the core axis—for three different 
fuel-burnup conditions: the initial condition of 
the batch-loaded core, the final condition (just 
before discharge) of the batch-loaded core, and 
the equilibrium condition of the out-in loaded 
core. It is to be noted that in both of the latter 
cases the maximum power density occurs at a 
considerable distance from the core axis but 
that for the out-in case the maximum power 
density lies on the central plane of the core, 
whereas in the batch case it lies above and below 
the central plane. 

It is probable that, in choosing fuel reload- 
ing programs, considerations of maximum-to- 
average power density and maximum-to-average 
burnup will be more important in most cases 
than considerations of fuel enrichment. 

Reference 5 reports a study of reactivity life- 
times of fuels composed of mixtures of U?*, 
u?** and thorium. The objective of the study is 
to determine whether mixtures of U** and tho- 
rium as fertile material would yield reactivity- 
exposure curves which are flatter than those 
obtained with either material separately and 
would, therefore, minimize the problems of 
reactivity change with fuel burnup. Space- 
independent calculations for many different 
mixtures were made by means of an analog 
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Figure 6—Two-dimensional contour plot of power 
density in one quadrant of a vertical plant bisecting 
the reactor core. (a) Batch loaded, initial. (b) Batch 
loaded, final (23,000 Mwd/ton average). (c) Out-in 
loading cycle, equilibrium (23,000 Mwd/ton average 
discharge burnup). 
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computer, using a fixed set of nuclear constants, 
corresponding to a slightly hardened spectrum 
at 650°K. The results of the study are sum- 
marized in the quotation: 


1. By varying the atom ratio of a fuel consisting 
of Th?*?, u?55, and U?8, it is possible to attain, for a 
specific reactor, an efficient use of the original fuel. 

2. Removal of fission products during the opera- 
tion greatly increases the lifetime of a fuel. 

3. For systems using Th”? as fertile material, 
very large reactivity lifetimes are possible for an 
initial conversion ratio significantly less than the 
theoretical maximum (n — 1) ascompared with sys- 
tems using us only as fertile material. 

4. The characteristics of the 7 variations in a 
Th’*?-yj2%5 fuel are almost a mirror image of the 7 
variations in U”*8-y25 fuel. Therefore it is con- 
cluded that areactor which is loaded with a mixture 
of natural uranium and thorium can have a flat curve 
of reactivity versus irradiation time over its life- 
time. The lifetime would then be determined theo- 
retically by irradiation damage. 

5. It is desirable to have as low an enrichment 
and as high a value of [y = €(1 — p)L] as is consist- 
ent with good design, where L is the neutron-leakage 
probability. 


Since, in this type of general study, questions 
of criticality must be ignored, the results must 
be studied in reference to a particular reactor 
type to determine whether the values of y which 
give favorable reactivity lifetime curves are 


practically attainable. In considering the practi- 
cal aspects of uranium-thorium mixtures, it 
must be recognized that whatever U***'is present 
in the discharged fuel will probably have its 
value degraded because of the alpha and gamma 
activities due to U**? and U?*. [The buildup of 
radioactivity in isotopes produced by neutronir- 
radiation of thorium is described briefly in 
Power Reactor Technology, 1(4): 19 (Septem- 
ber 1958), and 2(2): 42—43 (March 1959). | 
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In conventional reactor fuel elements, the fuel 
material is canned in some kind of nonfissiona- 
ble jacket to prevent contamination of the reac- 
tor coolant system, and for other reasons. Ifthe 
outside of the jacket becomes contaminated with 
fissionable material before the element is placed 
in the reactor, this function of the jacket will be 
partially negated. The contamination may or 
may not be a hazard, depending upon its level, 
upon the degree to which it can escape the reac- 
tor coolant system, and upon the ventilation and 
disposal systems in the reactor room. For ex- 
ample, in submarine applications, where venti- 
lating air is recirculated, a relatively small 
amount of contamination may prove trouble- 
some. In any case the presence of fission prod- 
ucts originating in surface contamination can 
mask the effects of fuel-element failures, and 
for this reason surface contamination is quite 
undesirable. 

The avoidance of important surface contami- 
nation can in some cases be difficult because 
exceedingly minute amounts of uranium suffice 
to give appreciable effects: in contamination 
studies, the quantity of fissionable isotope pres- 
ent On a surface is usually measured in units of 
10~° g per square decimeter of surface. 

Reference 1 contains the results of a number 
of studies which were made at KAPL to trace 
down the sources of surface contamination that 
had been observed on some industrially fabri- 
cated zirconium-uranium alloy fuel specimens. 
The reference states that Zircaloy-2 will nor- 
mally contain a small amount of natural uranium 
and that a limit of 5 ppm has been proposed for 
the permissibie natural-uranium content of 
Zircaloy-2 used for fuel-cladding material. 
With this natural-uranium content, there would 
be approximately 18 x 10~°* g of U’*®* per square 
decimeter of surface within a layer of 0.3-mil 
thickness, the approximate thickness through 
which fission products could be released. This 
surface density of U*® may be considered to 
define the quantity of surface contamination 
which would be considered important. The ref- 
erence states that some of the specimens ex- 
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amined had surface quantities of u**5 in excess 
of 100 x 107° g/dm’. A systematic study of the 
fuel-element fabrication and assembly steps was 
made to determine the source of the contamina- 
tion. Initially, this was done by withdrawing 
sample elements from the fabrication process 
and analyzing them for contamination. Thetests 
indicated:! 

1. Minor contamination was present on the 
fuel elements after form rolling. 

2. Major contamination was produced in vac- 
uum annealing. 

3. Pickling in HNO,;-HF (45 HNO;, 5 HF, 50 
H,O) to a depth of 0.002 in. removed substan- 
tially all contamination. 

4. Finished elements were uncontaminated. 


A second series of tests was made by passing 
clean finished elements through single fabrica- 
tion operations. These operations included the 
autoclave corrosion testing which was normally 
done on all finished fuel elements. The tests 
showed: 

1. Major contamination was again produced 
by vacuum annealing. 

2. No contamination was produced by pickling, 
even in highly contaminated acid. 

3. Minor contamination was produced in weld- 
ing. 

4. Major contamination was produced in auto- 
clave corrosion testing. The level of contami- 
nation varied widely from autoclave to auto- 
clave. 


Tests of the furnace and accessory surfaces 
used in the annealing revealed substantial U’* 
contamination on retort walls and annealing 
fixtures. Examination of the autoclave pressure 
vessels revealed no loose contamination, but the 
fixed contamination level on the interior walls, 
measured by alpha-particle counting, was in 
some cases as high as 19,500 dis/min/dm’. 
These measurements, along with controlled 
laboratory tests on small specimens, showed 
that, during autoclave testing, contamination 
could be transferred from autoclave walls to 
Zircaloy surfaces and that, after its transfer, it 
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was apparently fixed in place by the oxide film 
formed on the Zircaloy during the autoclave 
testing. Such contamination could not be re- 
moved to any appreciable extent by successive 
washes with 50-50 nitric acid at room tempera- 
ture. 

The findings of the investigation emphasize 
the necessity for a high degree of cleanliness in 
those pieces of fabrication and testing equipment 
which contain the fuel elements at high tempera- 


tures and particularly in the autoclaves, whose 
function is to test the elements and to forma 
stabilized oxide coating on them. 
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The General Electric Company, in its Atomic 
Power Equipment Department, is conducting a 
research and development project for the AEC 
on the use of solid neutron-absorbing materials 
for reactor control. The first topical report on 
this project! has recently been issued. It in- 
cludes summary descriptions of what is known 
and what has been done previously in this field 
and of the nature of current programs and gives 
subjects and recommendations for future in- 
vestigation. The document is a good general 
review of the control-rod materials field and 
contains extensive bibliographies which will be 
useful to the reactor designer. 

Since the document itself is a review, the 
following discussion will attempt only to cover 
the points of most general interest to reactor 
designers and to add enough introductory ma- 
terial to point up the reasons for the various 
requirements imposed on control-rod materials. 


Reactivity Requirements 
and Control-rod Effects 


The basic purpose of control rods, beyondthe 
adjustment of reactor power level, is to com- 
pensate the excess reactivity which must be 
built into the reactor to make it operable over 
its projected lifetime. Although the control 
rods may perform additional functions such as 
flux flattening or control of spatial oscillations 
of xenon content, the total reactivity worth re- 
quired of all the control rods in the reactor is 
usually determined by the maximum excess 
reactivity which the reactor may have during 
its lifetime. Often there is a strong desire to 
provide this reactivity worth by using the smal- 
lest possible number of control rods; for the 
cost associated with each rod (the cost of the 
rod itself, its drive mechanism, and its instal- 
lation) is fairly high, and each additional control 
rod adds to the complexity of the plant. Inthese 
cases the designer desires a control-rod ma- 
terial of the greatest possible nuclear effective- 
ness. In other cases, other considerations, 
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such as the use of control rods for power 
flattening, may dictate that the total reactivity 
holddown capacity be distributed amonga larger 
number of rods. The desire for materials of 
high nuclear effectiveness may then be re- 
duced, although the designer may still wish to 
use effective materials in order to keep the 
size of the rods small. 

In any case the nuclear effectiveness of a 
control-rod material is its key property, and 
the number of control rods necessary to com- 
pensate all the excess reactivity available is a 
key question in the design of a reactor. 

The major reactivity changes in power re- 
actors usually result from temperature and 
density changes in the reactor materials, the 
effect of fission products, and the burnup of 
fissionable isotope. Characteristic magnitudes 
of these changes vary from reactor type to 
reactor type, and each reactor within a type 
has its own characteristic reactivity changes. 
The H,O-moderated reactors usually experience 
large reactivity changes over their lifetimes, 
partly because of the large volume coefficient 
of thermal expansion of the moderator, H,O, 
which may cause a reactivity change of several 
per cent between room temperature and reactor 
operating temperature. Some H,O reactors ex- 
perience total reactivity changes of 25 per cent 
or more over the lifetime of a fuel charge. 
Other reactors may experience considerably 
smaller variations, but it is doubtful that many 
enriched solid-fuel power reactors of present 
types will experience reactivity changes less 
than 6 or 7 per cent Ros. If it is required that 
the normal control-rod installation be capable 
of making the reactor subcritical by, for ex- 
ample, 5 per cent under all conditions, then 
most enriched reactors will have control-rod 
installations worth some amount in the range 
from about 12 to 30 per cent kere, whereas the 
natural-uranium reactors, with less available 
reactivity, may fall below this range. 

To translate the control-rod “holddown’’ re- 
quirement into a specification of the number of 
control rods is usually a difficult problem in 
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reactor physics. Reference 1 gives a com- 
prehensive bibliography of control-rod calcula- 
tion methods as well as a brief review of the 
present status of control-rod theory. (A brief 
review of control-rod physics was given in the 
June 1958 issue of Power Reactor Technology, 
Vol. 1, No. 3. Papers which have appeared on 
control-rod theory since the review in refer- 
ence 1 are listed as references 2 to 8.) Cal- 
culations are sometimes rather badly in error, 
not so much because of any inadequacy of re- 
actor theory but because many practical con- 
trol-rod installations present complex physical 
situations which are not described well by the 
approximations the calculator is temptedtouse. 
In commenting on the many pitfalls of control- 
rod calculation methods that employ currently 
popular approximations, reference 1 states: 
“However, using methods of this type, an ex- 
perienced reactor physicist can usually com- 
pute worth of a given control system with an 
error of less than 50 per cent in spite of the 
fact that using an inappropriate recipe can 
easily result in errors of a factor of 2.”’ 

This comment should suffice as a warning 
that the following discussion of control-rod 
worth, which is based on even more simplified 
concepts, is useful only for illustrating the 
relation between reactor physics and the prob- 
lems of control-rod materials. 

The absorber control rod produces a re- 
duction of reactivity by absorbing some of the 
neutrons which otherwise would cause fission. 
The reactivity change produced by the absorp- 
tion of a given fraction of neutrons in a reactor 
depends upon where in the reactor they are 
absorbed and upon the energy at which they are 
absorbed. However, in “thermal’’ reactors most 
of the neutron absorption (both in control rods 
and in other materials) occurs in the thermal 
energy range, and in most power reactors the 
complement of control rods is sprinkled fairly 
uniformly over the reactor volume. Conse- 
quently some insight can be gained by con- 
sidering the case of pure thermal-neutron- 
absorbing rods in a pure thermal-neutron 
reactor, with the rods distributed evenly over 
the reactor volume. In such a case, if the con- 
trol rods absorb a fraction C of all neutrons 
absorbed in the reactor, then the reactivity 
worth of all the rods is very nearly equal to 
this fraction C. In this simplified case the 
estimation of total control-rod worth reduces 
to the estimation of the fractional absorption of 


thermal neutrons by a regular array of ab- 
sorbers in a reactor core medium which can, 
to a reasonably good approximation, he con- 
sidered infinite in extent. 

Ordinarily, control rods are desired to be 
nearly “black’’ to thermal neutrons. A rod is 
biack to neutrons if all neutrons that cross its 
surface boundary are absorbed. Approximate 
blackness to thermal neutrons requires two 
characteristics: (1) the thermal-neutron- 
absorption cross section must be many times 
larger than the scattering cross section, and 
(2) the absorber material must have a thickness 
of several absorption mean free paths. 

Figure 7 illustrates the changes in neutron 
absorption when a black, thermally absorbing, 
infinite slab is inserted in a slot, previously 
empty, in an infinite thermally absorbing me- 
dium which is fed by a thermal-neutron source 
distributed uniformly over the volume of the 
medium. Prior to the insertion of the slab, the 
neutron flux has a uniform value, $9, over the 


medium, the value being determined by the 
source strength and by the macroscopic ab- 


sorption cross section in the medium (see Fig. 
7(2). When the black slab is present, it absorbs 
some of the neutrons, and its absorption is 
compensated by an equal reduction of absorp- 
tion in the medium adjacent to the slab. The 
neutron flux near the slab drops in such a 
manner as to produce just the proper reduction 
of absorption in the medium. The resulting 
thermal-flux distribution, ¢(x), is shown in 
Fig. 7(b). The fractional difference between 
¢y) and ¢(x) at any point x (shown as a shaded 
area in the figure) may be regarded as the 
fractional absorption of neutrons from that 
point by the slab. A little way from the slab 
surface the curve ¢)— (x) takes the form of an 
exponential, decaying as e7*/L where L is the 
diffusion length of thermal neutrons in the 
medium surrounding the slab. If this behavior 
persisted right up to the surface of the slab and 
if the thermal-neutron flux reached zero at the 
surface, the area of the shaded region in Fig. 
7(@) would be just ¢)L on each side of the slab. 
Thus each side of the slab would absorb a num- 
ber of neutrons equal to the number produced 
by the source in that volume of the medium 
which lies within one diffusion length of the rod 
surface. Actually the exponential behavior does 
not persist to the slab surface, and the actual 
number of neutrons absorbed is somewhat less 
than this. The actual number can be computed 
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Figure 7— Typical change of thermal-neutron flux distribution when a black thermal-neutron ab- 
sorber (b) is inserted in an empty slot (a) in an infinite medium having a uniform source of thermal 


neutrons. 


by taking into account the extrapolation iength 
at the rod surface, which is determined by the 
transport mean free path in the surrounding 
medium. If, for example, the surrounding me- 
dium had the same thermal-neutron absorption 
and transport properties as the average core 
composition of the Experimental Boiling Water 
Reactor (EBWR),’ the actual absorption in the 
slab would be equal to those neutrons produced 
within a distance of about 0.8L of the surface. 

If other absorber shapes are considered, the 
results are rather similar in that the black ab- 
sorber captures a quantity of neutrons equal to 
those produced within some characteristic frac- 
tion of a diffusion length. The exact value of the 
fraction will depend upon the shape of the ab- 
sorber and upon the transport mean free path 
in the surrounding medium, but the variation in 
size of the fraction is unlikely to be large over 
the range of practical situations. 

In the reactor core the situation is more 
complex than that considered above, because 
the thermal neutrons are not supplied by a 
uniformly distributed source but by the slowing 
down of fission neutrons which were generated 
at each point in the core at a rate proportional 
to the thermal-neutron flux at that point. Con- 
sequently the production of thermal neutrons 
near the control rod will be somewhat lower 


than in regions far from the rod. This is 
usually not an extremely large effect, however, 
because the diffusion of the neutrons during the 
slowing-down process makes the rate of pro- 
duction of thermal neutrons a good deal more 
uniform than the distribution of fission density. 

The net effect of these considerations is that, 
for most control installations involving a fairly 
uniform lattice of control rods, a rough idea of 
the fractional absorption of thermal neutrons 
can be obtained by assuming that it is equal to 
the fractional core volume which lies within a 
distance aL of the control-rod surfaces. The 
a@ is a fraction which depends on rod shape and 
on the transport mean free path in the core; it 
will seldom be less than 0.5 or greater than 
0.8; a value of 0.7 might be representative for 
most reactors. 

This approach might be used to examine, for 
example, the control-rod absorption in the 
EBWR. In the region of the core occupied by 
control rods, each rod can be considered to be 
located at the center of a “cell’’ having a square 
cross section with a side length of 12°/, in. 
(see Fig. 8). The rod itself is a 10-in.-span 
cross. If one draws along the rod a line which 
is everywhere about 0.7L from the rod surface 
(L=1.8 cm), the area enclosed between this 
line and the surface of the rod is 135 cm’, 





eo ee Pe Ew YY we eh OS CY 


“= — 


— = a F 





CONTROL-ROD MATERIALS 29 


whereas the cross-sectional area of the cell, 
exclusive of the rod, is 1030 cm’. This rough 
approximation indicates that the rod absorbs 
about 13 per cent of the neutrons, and the re- 
activity change corresponding to this fractional 
absorption should be 13 per cent Refs. The 
quoted value of all the control rods in the 
EBWR is 12 to 14 per cent Reff.* Little im- 
portance can be attached to the quantitative 
agreement or lack of it; for, actually, rods are 
not distributed uniformly over the volume of 
the EBWR, no allowance has been made for 
epithermal effects, the core composition is not 
really uniform up to the surface of the rod, and 


‘the rod is not inserted into a nonabsorbing slot, 


but actually displaces water. Nevertheless the 
estimate does illustrate some of the factors af- 
fecting control-rod worth. The method of esti- 
mating can be used to correlate roughly the rod 
requirements for specified reactivity worths in 
thermal reactors of widely varying properties. 
Table VI-1 compares the estimated worths with 
reported worths for four quite different reac- 
tors. The worth estimates are made simply by 
examining a cross section of the core with con- 
trol rods inserted and taking the ratio of the 
area within 0.7L of the rod surfaces to the area 
of core external to rod surfaces. 

The factors affecting rod reactivity worth 
can be further illustrated by considering the 
absorption by cylindrical rods in reactors of 
widely different compositions. Figure 9 illus- 
trates the absorption by rods ina reactor having 
a low macroscopic absorption cross section in 
the core and a correspondingly long diffusion 
length. The medium to which the figure applies 
has average properties similar to those in the 
core of the original Calder Hall reactor.'? The 
moderator is graphite, and the diffusion length 
in the core is 19.8 cm. If a black rod 2 in. in 
diameter is installed in such a medium, it ab- 
sorbs a number of thermal neutrons equal to 
those formed within a distance of about 13.4 
cm of the rod surface. If the rod diameter is 


increased to 3 in., the effective distance from 
which neutrons are drawn increases only 
slightly, to 14.4 cm. In considering these dia- 
grams, three points are obvious: 





* Reference 9. The worth listed in the reference 
(calculated) is given as the value of Refs which the 
rods will hold just critical. The rod worth used here 
is (kp — L/P ore: 


1. The volume of the core from which neu- 
trons are drawn by the rod is large compared 
with the volume of the rod. In this respect the 
rod is quite effective in influencing reactivity, 
and only a relatively small fraction of the core 
volume would have to be occupied by control 
rods in order to achieve a high degree of 
control-rod worth. 

2. The number of neutrons absorbed by the 
rod does not increase in proportion to the size 
of the rod. The reason for this is obvious. Be- 
cause of this effect, most low-cross-section 
reactors employ rods which are quite small 
relative to the size of the reactor. 

3. Evidently it would not be profitable to in- 
stall control rods very close together in a low- 
cross-section reactor, for then adjacent rods 
would have zones of effective absorption which 
would overlap and the total number of neutrons 
absorbed per rod would decrease. Actually the 
space necessary to avoid important mutual 
“shadowing’’ of the rods is somewhat greater 
than the distance indicated by the circles sur- 
rounding the rods in the figure. It would be 
unusual to find rods closer together than two 
or three diffusion lengths in any reactor. 
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Figure 8—Typical control-rod cell in EBWR. (a) 
Contours of constant thermal-neutron flux (adapted 
from reference 9). (b) Approximate area of effective 
absorption. 
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Table VI-i REACTIVITY WORTH OF CONTROL-ROD INSTALLATIONS IN SEVERAL REACTORS 








B 
(cross- 

A sectional Reported 

(cross- area within worth of 

Core Control- Control- No. of Core sectional distance Estimated control 

diam., rod rod control diffusion area of 0.7L of rod rod worth, rods 
ft shape size rods length, cm core), cm? surfaces), cm? B/A (Akers Ress) 
EBWR® (Experimental Boiling Water Reactor), Boiling H,O Type 

4 Cross 10-in. 9 1.8 11,700 1,240 0.11* 0.12-—0.14f 


span 


Adam" (Swedish Reactor for Heating), Natural Uranium Oxide -D,O Type 


a2 34 Cylindrical 90-mm 32 14.8 


diam. 


91,000 20,200 0.22 0.15 


SRE" (Sodium Reactor Experiment), Sodium Graphite Type 


6 Cylindrical 2.484-in. 8 12.9 


diam. 


26,300 3,480 0.13 0.135 


Calder Hall,!® Natural Uranium —-Graphite Type 


31 Cylindrical 1.75-in. 48 19.8f 


diam. in 
3.25-in. 
channel 


703,000 36,500 0.052 0.068 





* This value differs from the value estimated in the text because there is some core area which is not 
included in the nine ‘‘typical cells’’ surrounding the control rods. 

tIn reference 9 the (calculated) worths of the rodsare given as the value of keg; which the rods will hold 
just critical. To get the rod worth, as used here, take ker; — 1/keg, - 


t Intermediate core region. See reference 13. 
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Figure 9— Relative numbers of thermal neutrons ab- 
sorbed by 2- and 3-in. cylindrical black rods when 
placed in identical infinite media supplied by uniformly 
distributed sources of equal strength. Neutron trans- 
port and absorption properties are the same as in 
region B of Calder Hall, homogenized. 
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Figure 10—Relative numbers of thermal neutrons 
absorbed by 2- and 3-in. cylindrical black rods when 
placed inidentical infinite media supplied by uniformly 
distributed sources of equal strength. Neutron trans- 
port and absorption properties are the same as in a 
uniform, cold EBWR core.2 
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Figure 10, which is necessarily drawn to a 
quite different scale, illustrates the absorp- 
tions by the same two rods ifthey were installed 
in a highly absorbing reactor core. The actual 
medium considered has the average properties 
of the cold EBWR core.’ The main character- 
istics, which are as follows, are quite different 
than in the previous case: 


1. The volume of core from which neutrons 
are drawn by the rod is not large relative to 
the rod itself. 

2. The increase of absorption with rod size 
is more nearly proportional to the increase in 
surface area of the rod. 

3. Since the area of effective absorption of 
the rods is small, rods may be installed rela- 
tively close together without important mutual 
shadowing. 


Figure 10 illustrates some of the difficulties 
of devising efficient rod installations for water- 
moderated reactors. Quite evidently the ther- 
mally absorbing cylindrical rod is aninefficient 
means of compensating excess reactivity, for it 
occupies a rather large amount of core volume 
per neutron absorbed. Furthermore, when the 
rod is withdrawn from the core, the empty 
channel which it leaves poses a problem. If 
the channel is allowed to fill with water, the 
fast neutrons slowed down in the water will 
produce a peak in the thermal-neutron flux 
which will tend to overheat the adjacent fuel 
elements; and, in addition, the thermal-neutron 
absorption of the water, enhanced by the peak- 
ing of the thermal flux, may reduce the net re- 
activity gain resulting from withdrawal of the 
control rod. These are the reasons that rod 
shapes giving high surface-to-volume ratios, 
such as crosses, are used in water-moderated 
reactors. In the water reactor the diffusion 
length is so small that the rod surface must be 
brought into proximity to a substantial fraction 
of the uranium in the core in order to absorb a 
substantial fraction of the neutrons. Even when 
such extended surface rods are used, it is fre- 
quently difficult to arrive at a mechanical de- 
sign of the core which will place this surface 
adjacent to a sufficiently large fraction of the 
uranium. For that reason it is desirable in 
water reactors to increase the basic nuclear 
effectiveness of the control-rod material. Once 
the macroscopic absorption cross section of the 
material has been made sufficiently high that 
the rod is black to thermal neutrons, further 





increases do not cause additional absorption of 
thermal neutrons. The only remaining way that 
the effectiveness of the material can be in- 
creased is by causing it to absorb the neutrons 
whose energies are higher than thermal. 

It is difficult to absorb a large fraction of the 
epithermal neutrons, and there is no great in- 
centive to absorb epithermal neutrons unless 
the number of epithermal neutrons normally 
present in the reactor core is important rela- 
tive to the number of thermal neutrons. Ordi- 
narily those reactors which have highly absorb- 
ing cores, notably the water (H,O) reactors, 
will have high ratios of epithermal-to-thermal 
neutron flux. Consequently the use of epither- 
mal control-rod absorption in water reactors 
is favored for two reasons: it is needed to at- 
tain adequate rod worth; epithermal absorption 
in such reactors can be reasonably effective 
in reducing reactivity. 

There is a scheme which will make a control 
rod effectively an epithermal absorber even if 
the absorbing material of the rod has a high 
cross section only in the thermal energy range. 
This scheme is to use a relatively thick rod and 
to incorporate effective moderator in the rodso 


that some of the neutrons which enter the rod 
as epithermal neutrons will be slowed down to 


thermal energy within the rod and will be ab- 
sorbed. This scheme can be illustrated by Fig. 
10 if it is assumed that the rod consists of a 
cylindrical shell of absorbing material, just 
thick enough to be thermally black, and that the 
volume inside the shell is allowed to fill with 
water. Reasonably large control rods of this 
type can be very effective, for the pure water 
inside the shell is even more effective in ther- 
malizing neutrons than the mixture of water and 
other materials which makes up the core 
proper. The scheme poses the difficulty of the 
water-filled channel which exists when the rod 
is withdrawn. This difficulty can be circum- 
vented by providing a follower for the rod made 
of fuel elements similar to those in the rest of 
the core. These elements must, of course, be 
adequately cooled, and this may cause engineer- 
ing complications. However, the rod type has 
been used in the Materials Testing Reactor 
(MTR) and the Army Package Power Reactor 
(APPR).* It is sometimes referred to as the 
neutron rectifier control rod. 





* It is also used in the Argonne conceptual design of 
a prototype for an advanced boiling reactor, reviewed 
elsewhere in this issue. 
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If the scheme of incorporating moderator in 
the control rod is not used, absorption of epi- 
thermal neutrons can still be accomplished by 
the use of control-rod materials which have 
high epithermal absorption cross sections as 
well as high thermal absorption cross sections. 
Most materials considered for control rods 
have sufficiently high thermal cross sections 
that a reasonable thickness of the material is 
black or nearly black tothermal neutrons. Con- 
sequently the differences in effectiveness which 
different materials exhibit are principally due 
to differences in their epithermal cross sec- 
tions, and these differences will show up with 
greater or lesser importance accordingly as 
the ratio of epithermal-to-thermal neutron flux 
in the reactor is higher or lower. 


Effectiveness of 
Control-rod Materials 


The thermal absorption cross section of a 
material determines how thick a layer of the 
material must be used in order to be black to 
thermal neutrons. For practical purposes a 
layer of thickness that is two absorption mean 
free paths or greater can be considered nearly 
black. Stevens'4 has shown that an absorbing 
layer of thickness equal to two absorption mean 
free paths absorbs something like 97 per cent 
as many neutrons as a completely black sur- 
face. Table VI-2 lists a number of elements and 
isotopes with high thermal-neutron-absorption 
cross sections and gives for each the surface 
density, or number of grams of absorber per 


Table VI-2 AMOUNTS OF VARIOUS ABSORBERS REQUIRED FOR NOMINAL BLACKNESS 
TO THERMAL NEUTRONS" 


(Two Absorption Mean Free Paths) 








Average thermal 


Abundance of 


Surface density 





Element or cross section, important for blackness, Density, Thickness for 
nuclide barns isotopes, % g/cm? g/cm® blackness, cm 
Li 62 0.362 0.53 0.68 
Lié 818 7.52 0.024 0.53 0.045 
B 651 0.054 2.3 0.023 
B! 3,470 18.8 0.0093 2.3 0.004 
Rh 130 2.56 12.5 0.20 
Ag 54 6.47 10.5 0.62 
Agi"? 26 51.4 
Ag! 71 48.6 
Cd 2,210 0.165 8.6 0.0192 
calls 18,000 12:8 
In 165 2.26 7:3 0.310 
Sm 4,760 0.100 7.2% 0.0139 
Sm'48 57,200 13.8 
Eu 3,980 0.128 7.3% 0.0175 
Eu!t 7,800 47.8 
39,800 0.0128 7.6* 0.0017 
Gat 60,000 14.7 
Ga‘ 139,000 15.7 
Dy 950 0.55 8.1* 0.068 
Dy! 2,340 28.2 
Er 144 3.78 8.6* 0.44 
Lu 96 6.05 ? ~ 0.67 
Lu!’6 3,120 2.6 
Hf 91 6.38 13.4 0.471 
Hf!" 320 18.4 
Re 73 8.28 20: 0.414 
Ir 372 1.68 22.4 0.075 
Ir'*! 820 38.5 
Ir’ 104 61.5 
Au 85 7.52 19.3 0.392 
Hg 330 1.97 14.22 0.138 
Hg'*® 2,160 16.8 





* X-ray density for oxide. 
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square centimeter of surface, required to meet 
the blackness criterion of two absorption mean 
free paths. 

It must be remembered that a further crite- 
rion for blackness is that the absorption cross 
section of the material be much larger than the 
scattering cross section. Thus large thicknesses 
of moderately absorbing materials do not make 
very effective control rods, and the effective- 
ness of a highly absorbing material may be 
impaired if it is diluted too strongly with a 
material which is a good scatterer. Stevens’! 
cites the case of a B’® slab of thickness equal 
to two absorption mean free paths. This slab 
would fail to capture about 6 per cent of the 
neutrons crossing its surface because of its 
slight transparency. Such a slab would, of 
course, be quite thin. If it were diluted by 
steel sufficiently to form a slab 0.20 in. thick, 
the fraction of neutrons which avoided capture 
would increase to 8.3 per cent; some of these 
would pass through the slab, whereas others 
would be scattered out of the slab before ab- 
sorption. This degree of dilution increases the 
scattering cross section to 20 per cent of the 
macroscopic absorption cross section. It ap- 
pears that scattering cross sections of 10 to 
20 per cent of the absorption cross section are 
tolerable, provided the absorbing and scattering 
nuclei are intimately mixed so that the com- 
petition between scattering and absorption is in 
proportion to the macroscopic absorption and 
scattering cross sections. If the scattering and 
absorbing materials are segregated, the scat- 
tering loss of neutrons may be considerably 
greater. The worst form of segregation is, of 
course, that in which the scattering material 
covers the surface of the absorbing material, 
as when the absorber is encased in a jacket or 
cladding of low absorption cross section. Stevens 
shows that a 0.040-in. jacket of steel reflects 
about 8 per cent of the neutrons incident upon 
the jacket, whereas a 0.10-in. thickness of steel 
reflects about 17 per cent. Theactual reduction 
of rod effectiveness by the scattering effect of 
the jacket will be less than these percentages 
of reflection and will depend upon the surround- 
ing medium, for some of the reflected neutrons 
will diffuse back into the rod. Further, any 
material which lies between the control rod and 
the fissionable material in the reactor core can 
produce a somewhat similar effect. If the ma- 
terial between the rod surface and the fission- 
able material is a moderator, then the net effect 


on rod effectiveness is the difference between 
the effect of the diffusion barrier which the 
material represents and the effect of additional 
neutron slowing down in the vicinity of the rod 
surface. 

Since the contribution of epithermal-neutron 
absorption to control-rod worth depends upon 
the neutron energy spectrum of the reactor — 
principally upon the ratio of the epithermal- 
to-thermal neutron flux — it is difficult to gen- 
eralize upon the desirable characteristics of 
epithermal absorbers other than to say that the 
best materials will be those having high cross 
sections over large ranges of the epithermal 
energy region. The three most widely used 
control-rod absorbers (boron, cadmium, and 
hafnium) illustrate the possible range of epi- 
thermal absorption behavior. Figure 11, re- 
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Figure 11— Fraction of neutrons absorbed by various 


produced from the data of Stevens,‘ shows the 
fractional absorptions by slabs containing these 
materials as functions of neutron energy. The 
figure gives the fraction of neutrons absorbed, 
of those striking the slab surface, for four dif- 
ferent slabs: a '¢-in.-thick slab of cadmium, 
a 0.200-in.-thick slab of hafnium, a 0.200-in.- 
thick slab of steel containing 1.0 wt.% B", and 
a 0.200-in.-thick slab of steel containing 0.15 
wt.% B'®. The latter case would correspond to 
about 0.9 per cent natural boron in the steel. 
Cadmium is the classic example of the thermal 
absorber. The '/,,-in. thickness is black to 
neutrons up to an energy of about 0.4 ev. At 
this energy its cross section is decreasing very 
rapidly with neutron energy, and an increase in 
thickness of the slab is not very effective in 
increasing its neutron absorption: a doubling of 
the slab thickness would extend the range of 
blackness up to about 0.5 ev. 
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Figure 12—Neutron-flux spectrum in H,O containing 
2.9 barns of 1/v absorber per H atom.!® 


Since hafnium is not a very strong thermal 
absorber, the 0.2-in.-thick slab is not really 
black in the thermal energy region, although 
most of the thermal neutrons are absorbed. The 
slab is black in the vicinity of the absorption 
resonances which occur at 1.10, 2.38, and 
7.80 ev. An increase in thickness of the haf- 
nium would increase its absorption in the ther- 
mal energy range and in the epithermal ranges 
between the neutron resonances. 

Boron-10 is the classic example of the 1/v 
absorber. An increase in either the thickness 
of the slab or in the concentration of B*° ex- 
tends the range of blackness to higher and 
higher neutron energies. The 0.15 per cent 
B’® slab is about equivalent to the hafnium slab 
for neutrons of energies up to about 0.3 ev; at 
higher energies, it falls far below the hafnium 
curve. The 1 per cent B” slab, however, is 
black up to about 0.7 ev and is about as good as 
the hafnium slab at higher energies. 

Most epithermal absorbers are of the reso- 
nance type, like hafnium. Since these materials 
characteristically absorb in discrete energy 
regions, one of the methods of producing ma- 
terials with high epithermal absorption effi- 
ciency is to mix absorbers in such a way that 
the resonances of one fall into the energy re- 
gions between the resonances of another. The 
mixture of cadmium, silver, and indium is one 
such combination of this type. Mixtures of the 
rare-earth elements have also been considered. 

The importance of epithermal absorption in 
H,O reactors can be inferred from Fig. 12, 
which shows the energy distribution of neutron 
flux in a typical water-moderated medium, 
after the calculations of Amster.'® This flux 


spectrum, which was computed for water con- 
taining 2.9 barns of absorption cross section 
per hydrogen atom (boron absorption), is much 
like the spectrum that would be expected ina 
water- moderated reactor having an H,O/UO, vol- 
ume ratio of 2.0, with the uranium enriched to 
2 per cent U**. Since the importance of the 
epithermal energy range is distorted by the 
logarithmic energy scale, a better idea of the 
importance of the epithermal neutrons can be 
had from the integral, [¢(£) dE, which is also 
plotted in Fig. 12. [¢(£) and [¢(E) dE are plotted 
to different arbitrary scales.] It is evident that 
a substantial portion of the neutron total flux 
lies above the cadmium cutoff at 0.4 ev. 

The neutron energy spectrum and the ab- 
sorption-cross-section curve are not the only 
quantities involved in determining the epithermal 
effectiveness of an absorbing material. Just 
as, in the case of purely thermal absorption, 
the control rod must compete against the ab- 
sorption of the surrounding medium, so, in the 
epithermal case, the control rod must compete 
both against whatever epithermal absorption 
occurs in the surrounding medium and against 


Table VI-3 CONTROL-ROD EFFECTIVENESS* 
FOR VARIOUS MATERIALS! 








Effectiveness 
Absorptiveness, relative to 
Material z,tt hafnium 

3.0 wt.% B" in stainless 

steel (dispersion of 

minus 100-mesh par- 

ticles of 90 per cent 

enriched boron) 27.2 £32 
Hafnium 2.40 1.00 
0.97 wt.% B'° in stainless 

steel (alloy) 9.3 0.98 
15 wt.% Eu,O, in stainless 

steel (dispersion) 8.7 0.96 
Indium 3.68 0.93 
Silver 1.85 0.88 
Cadmium 77.5 0.80 
8.7 wt.% gadolinium in 

titanium 31.1 0.79 
Tantalum 0.60 0.71 
2.7 wt.% Sm,O; in stainless 

steel (dispersion) 3.20 0.71 
Haynes-25 alloy 0.96 0.68 
Titanium 0.16 0.24 
Zircaloy-2 0.05 
2S aluminum 0.02 





* Determined by R. R. Eggleston and a Critical Experi- 
ment Group of Bettis Laboratory on samples 2 by 2 by 0.2 in. 

+ Macroscopic thermal absorption cross section times 
thickness. 
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Table VI-4 MEASURED RELATIVE REACTIVITY WORTHS OF VARIOUS MATERIALS 
IN AN H,O-MODERATED ASSEMBLY 





Equivalent boron surface 








Surface Relative Relative ; ° 
density, Thickness, worth, * worth, * density, g/em 
Material g/cm? in. bare Cd covered Bare Cd covered 
Cadmium 0.022 0.001 0.494 0.931 0.0090 + 0.0005 
(Cd) 0.050 0.002 0.648 0.931 0.0192 + 0.0010 
0.101 : 0.005 0.728 0.928 0.0284 + 0.0015 
0.555 0.025 0.811 0.941 0.0440 + 0.0026 0.0060 + 0.0051 
| 1.104 0.051 0.838 0.948 0.0510 + 0.0032 0.0100 + 0.0053 
. 2.142 0.098 0.864 0.967 0.0593 + 0.0040 0.0200 + 0.0056 
| 4.389 0.201 1.000 0.0370 + 0.0063 
6.585 0.305 0.928 1.020 0.0890 + 0.0066 0.0500 + 0.0071 
Cobalt 5.601 0.252 0.855 1.046 0.0563 + 0.0038 0.0810 + 0.0105 
(Co) 3.418 0.154 0.759 1.011 0.0345 + 0.0018 0.0510 + 0.0072 
2.187 0.098 0.663 0.988 0.0206 + 0.0010 0.0323 + 0.0061 
Dysprosium 3.047 0.300 1.039 1.158 0.2050 + 0.0190 0.2400 + 0.0250 
oxide 1.753 0.174 " 0.989 1.119 0.1380 + 0.0120 0.1660 + 0.0185 
(Dy20s) 0.505 0.050 0.801 1.034 0.0420 + 0.0024 0.0710 + 0.0087 
Erbium 3.405 0.299 0.920 1.077 0.0840 + 0.0061 0.1060 + 0.0130 
oxide 2.009 0.175 0.839 1.026 0.0520 + 0.0033 0.0650 + 0.0083 
(Er,03) 0.558 0.049 0.570 0.979 0.0137 + 0.0007 0.0290 + 0.0060 
Europium 3.375 0.289 1.253 1.332 1.0900 + 0.0710 1.0000 + 0.1150 
oxide 1.972 0.174 1.197 1.265 0.6850 + 0.0570 0.6100 + 0.0600 
(Eu,03) 0.561 0.050 1.032 1.130 0.1940 + 0.0180 0.1850 + 0.0200 
Gadolinium 3.388 0.301 1.042 1.135 0.2100 + 0.0190 0.1950 + 0.0210 
oxide 1.989 0.177 0.984 1.082 0.1330 4 0.0110 0.1180 + 0.0140 
(Gd,03) 0.567 0.050 0.885 1.011 0.0680 + 0.0047 0.0510 + 0.0072 
Gold 14.531 0.297 1.038 1.164 0.2040 + 0.0190 0.2540 + 0.0260 
(Au) 9.606 0.196 0.987 1.127 0.1360 + 0.0110 0.1800 + 0.0200 
4.926 0.101 0.868 1.076 0.0610 + 0.0041 0.1100 + 0.0135 
2.518 0.051 0.719 1.032 0.0273 + 0.0014 0.0700 + 0.0088 
1.318 0.027 0.564 1.004 0.0127 + 0.0007 0.0460 + 0.0068 
Hafnium 0.175 0.276 1.106 1.217 0.3530 + 0.0360 0.4050 + 0.0390 
(Hf) 6.833 0.206 1.051 1.189 0.2220 + 0.0200 0.3200 + 0.0310 
3.422 0.103 0.920 1.130 0.0840 + 0.0062 0.1890 + 0.0205 
1.476 0.045 0.726 1.072 0.0284 + 0.0015 0.1050 + 0.0130 
0.856 0.025 0.572 1.039 0.0133 + 0.0007 0.0760 + 0.0096 
Holmium 3.365 0.302 0.745 1.087 0.0314 + 0.0016 0.1280 + 0.0150 
oxide 1.955 0.176 0.602 1.038 0.0152 + 0.0008 0.0740 + 0.0094 
(Ho203) 0.564 0.051 0.302 0.983 0.0036 + 0.0003 0.0310 + 0.0059 
Indium 5.554 0.301 1.036 1.137 0.2000 + 0.0190 0.2000 + 0.0215 
(In) 3.697 0.201 0.998 1.105 0.1470 + 0.0120 0.1470 + 0.0167 
1.858 0.100 0.908 1.058 0.0785 + 0.0057 0.0940 + 0.0120 
| 0.928 0.050 0.785 1.029 0.0386 + 0.0021 0.0660 + 0.0086 
0.468 0.025 0.632 1.003 0.0176 + 0.0009 0.0450 + 0.0068 
| Lutetium 2.244 0.175 0.772 1.059 0.0361 + 0.0019 0.0930 + 0.0120 
oxide 1.605 0.125 0.693 1.041 0.0241 + 0.0012 0.0780 + 0.0099 
(Lu,O3) 0.638 0.050 0.465 0.992 0.0078 + 0.0004 0.0360 + 0.0062 
Samarium 3.357 0.302 1.063 1.155 0.2490 + 0.0240 0.2340 + 0.0240 
oxide 1.966 0.175 1.018 1.113 0.1730 + 0.0150 0.1580 + 0.0173 
(SmgOs3) 0.558 0.050 0.896 1.008 0.0730 + 0.0052 0.0500 + 0.0070 
Silver 7.970 0.299 1.021 1.172 0.1760 + 0.0155 0.2710 + 0.0270 
(Ag) 5.226 0.196 0.944 1.123 0.0995 + 0.0076 0.1700 + 0.0190 
2.746 0.103 0.798 1.065 0.0414 + 0.0023 0.0990 + 0.0125 
1.376 0.052 0.617 1.022 0.0167 + 0.0009 0.0600 + 0.0080 
Tantalum 12.733 0.302 0.855 1.166 0.0680 + 0.0047 0.2610 + 0.0265 
(Ta) 8.511 0.202 0.775 1.124 0.0370 + 0.0020 0.1740 + 0.0190 
4.214 0.100 0.572 1.059 0.0133 + 0.0007 0.0980 + 0.0120 








(Table Continues) 
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Table VI-4 (Continued) 
Equival f 
Surface Relative Relative greg — 
density, Thickness, worth,* worth,* hd 
Material g/cm? in. bare Cd covered Bare Cd covered 

Terbium 2.287 0.175 0.829 1.067 0.0482 + 0.0030 0.1010 + 0.0135 
oxide 1.624 0.125 0.765 1.040 0.0347 + 0.0019 0.0730 + 0.0090 
(Tb,O7) 

Thulium 2.089 0.175 0.735 1.056 0.0292 + 0.0015 0.0890 + 0.0110 
oxide 1.490 0.125 0.649 1.037 0.0193 + 0.0010 0.0730 + 0.0090 
(Tm,03) 0.580 0.049 0.418 0.992 0.0064 + 0.0004 0.0370 + 0.0063 

Tungsten 9.900 0.204 0.693 1.077 0.0257 + 0.0013 0.1120 + 0.0136 
(W) 9.581 0.197 0.675 1.053 0.0237 + 0.0012 0.0890 + 0:0110 

Ytterbium 3.584 0.300 0.579 1.028 0.0138 + 0.0007 0.0670 + 0.0084 
oxide 2.100 0.175 0.440 0.994 0.0070 + 0.0004 0.0345 + 0.0061 
(Yb,0Os) 0.581 0.050 0.194 0.957 0.0017 + 0.0002 0.0150 + 0.0054 

Boron glass 0.04787 0.101 0.8144 1.0051 

0.1575 0.104 1.0062 1.1161 
0.0335 0.152 0.7503 0.9903 
0.0880 0.200 0.9282 1.0560 
0.3026 0.206 1.0867 1.1806 
0.1206 0.275 0.9728 1.0882 
0.4241 0.290 1.1326 1.2232 
0.0664 0.302 0.8756 1.0393 





* Edge corrected. 
+ Surface densities of boron. 


the moderating neutron collisions which tend to 
slow the neutrons down past the control-rod 
absorption resonances. These effects can, of 
course, be evaluated theoretically, but experi- 
mental results are more convenient for survey- 
ing the approximate ranges of improvement of 
effectiveness that can be expected from the use 
of good epithermal absorbers. 

Table VI-3 lists measured effectiveness for 
various control-rod materials relative to a 
0.20-in.-thick hafnium sheet, as reported in 
reference 1. The measurements were made 
with small samples, 2 by 2 by 0.2 in. The 
spectrum in which the measurements were 
made is not specified, but it can probably be 
assumed that it is a typical water-moderated 
spectrum and that the relative values of ef- 
fectiveness would be roughly the same for many 
water-moderated reactors. 

Table VI-4 lists the results of relative worth 
measurements made in the Critical Experiment 
Facility at Vallecitos Atomic Laboratory.'® The 
materials were used in the form of 3- by 6-in. 
plates having various thicknesses as recorded 
in the table. The critical assembly was H,O 
moderated; it was fueled with UO, pellets, 
0.491 in. in diameter, in 0.562-in.-O.D. alu- 
minum tubes of 30-mil wall thickness. The fuel 
rods were located on a uniform-square lattice, 


with 0.8-in. spacing between centers; the en- 
richment was 1.6 wt.% U**. Measurements were 
made of the reactivity worth of each absorbing 
plate, bare, and of each plate plus a constant 
cadmium box enclosing it. The relative worth 
is expressed, arbitrarily, as the ratio of the 
reactivity worth of the arrangement in question 
to the reactivity worth of a 0.201-in.-thick 
cadmium plate in the cadmium box. The results 
have been edge corrected to reduce the edge 
effect to that which would characterize a plate 
of zero thickness. The equivalent boron surface 
density is that surface density of natural boron 
which, if substituted for the sample (bare or in 
the cadmium box as the case might be), would 
yield the same reactivity worth. The epither- 
mal effectivenesses of the materials are illus- 
trated by their cadmium-covered relative worths 
and by the cadmium-covered equivalent boron 
surface densities. 


Other measurements on rare-earth oxides 
and mixtures of rare-earth oxides may be found 
in reference 17. 


Burnup of Control-rod Materials 


Some control rods spend little time in the 
operating reactor; they may be used simply to 
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provide extra shutdown capacity or to compen- 
sate such reactivity excesses as the inverse 
temperature defect or the increase due to the 
decay of equilibrium xenon, which do not re- 
main long in the operating reactor. Other rods 
are exposed to the reactor flux over long pe- 
riods during their lifetime. They may be rods 
which compensate the excess reactivity pro- 
vided for fuel burnup or for transient xenon 
override. These rods may suffer from radia- 
tion damage or from burnup of the absorbing 
isotope. Radiation damage is usually important 
only in those materials (boron and lithium) 
which undergo (m,q@) reactions. The effects are 
discussed for boron in a later section. When 
appreciable burnup of the absorbing isotope is 
anticipated, extra material, beyond that needed 
for blackness, must be incorporated in the rod. 

Whether deterioration of the rod is due to 
radiation damage or to burnup, its extent is 
determined by the number of neutron absorp- 
tions in the rod. A rough estimate of the num- 
ber of absorptions can be made by recalling 
that the reactivity worth of the rod, at a posi- 
tion of average importance in the reactor, is 
just equal to the fraction of total neutron ab- 
sorptions in the reactor which take place in the 
rod. For criticality reasons most of the ab- 
sorptions in the reactor must take place in 
u*> or other fissionable isotopes. In almost 
all cases the fractional absorption in the fis- 
sionable isotopes will lie in the range 0.6 to 
0.9. The destruction of 1 gram-atom of U’*® 
yields a thermal energy output of about 186 
Mwd. Consequently an estimate of the average 
number of gram-atoms of control-rod absorber 
destroyed is given by 


Gram-atoms of absorbing isotope 


destroyed = — x 


where p is the worth of the control rod or con- 
trol rods in question, £ is the energy, in ther- 
mal megawatt-days, which the reactor produces 
while the rod is inserted, and the most pes- 
simistic value of the fractional absorption in 
fissionable isotope has been used. Thus, if a 
B” rod installed in a position of average im- 
portance in a reactor of 100 Mw(t) output is 
worth 0.25 per cent Ror, the burnup of boron 
in the rod over a period of 10 years of equiva- 
lent full power is approximately 


(0.0025) (100) (10) (365) 
112 





Burnup = 


= 8.2 gram-atoms 


or about 82 g of B’. It must, of course, be 
remembered that a rod in a position of lesser 
importance will absorb more neutrons for a 
given rod worth. Further, the burnup in the rod 
will be distributed along the length of the rod 
in proportion to the power distribution in the 
reactor along the rod axis. If the rod is only 
partially inserted, there will be a further con- 
centration of absorption near its tip. The lat- 
ter effect may be very large if a number of 
strong rods are moved in unison as a bank. 

When appreciable burnup in the rods must be 
anticipated, there is little benefit in the use of 
absorbing materials having extremely high 
microscopic thermal absorption cross sections; 
for, although quite small quantities of such 
materials may suffice for blackness, the burn- 
up requirement may set the specification for 
the amount of absorber incorporated in the rod. 
If neutron absorption in the absorber yields a 
stable or long-half-life isotope which itself has 
a high absorption cross section, then, of course, 
the burnup requirement is reduced. Europium, 
for example, has two isotopes (151 and 153), 
both of which yield high-cross-section isotopes 
upon neutron absorption. 


Other Requirements 
of Control-rod Materials 


Of the remaining requirements for control- 
rod materials, the two most important, other 
than the obvious mechanical requirements, are 
compatibility with the environment within the 
reactor and susceptibility to adequate cooling. 
Little can be said in a general way on the first 
of these requirements; for a given reactor, the 
characteristics of the materials being con- 
sidered must be checked in detail against the 
conditions imposed by the environment. The 
two materials that have been used fairly widely 
without protective claddings are hafnium and 
borated stainless steel. Both are compatible 
with water at the temperatures characteristic 
of boiling- and pressurized-water reactors. It 
is, of course, quite feasible to use materials 
which must be jacketed against the environ- 
ment. However, jacketed rods sometimes in- 
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troduce problems in addition to that of main- 
taining the integrity of the jacket. They may 
increase the mechanical engineering problems, 
decrease somewhat the absorption effective- 
ness, and increase the problems of cooling. In 
H,O reactors, where there is a strong desire 
to keep control-rod sections as thin as possible, 
a jacket may add an undesirable extra thickness. 


There are three possible sources of control- 
rod heating: the energy produced by the neutron- 
absorption reaction in the rod, the energy de- 
posited by fission gamma radiation from the 
surrounding fuel, and the energy deposited by 
gamma radiation generated in the inelastic col- 
lisions of neutrons. Of these sources, the 
latter is of the least consequence and the first 
is usually the most important, both because it 
usually produces the highest average energy 
deposition in the rod and because it may be 
more highly concentrated, particularly at the 
ends of partially inserted control rods. In 
(n,@) absorbers (boron) practically all of the 
energy of the absorption reaction is deposited 
in the control rod. In (n,y) absorbers an ap- 
preciable fraction of the gamma-radiation 
energy escapes from the control rod. In either 
case the heat-generation rate in control rods in 
high-power-density H,O-moderated reactors 
can be quite high. Adequate coolant flow must 
be provided for the rod, and, in the case of 
jacketed absorbers, adequate provision must 
be made for heat transfer from the absorber 
to the jacket. 


Although it is important to assess the prob- 
lem of control-rod heating adequately in high- 
power-density reactors and to provide for 
adequate cooling, it is also important to recog- 
nize that in some lower power-density applica- 
tions the problem of control-rod cooling may 
be a trivial one. Whenever the problem of 
control-rod heating is not great, the incentive 
for control materials compatible with the re- 
actor-core environment is considerably re- 
duced, and it may be possible to produce 
adequate control rods quite inexpensively by 
simply canning an absorbing material in an 
appropriate jacket without great regard for the 
problems of heat removal. 


Boron Steel! 


The boron-stainless steel alloys are rela- 
tively cheap control materials which have been 


used for control rods in a number of applica- 
tions, including applications in high-temperature 
water. The boron content used usually lies in 
the range between 1 and 2 wt.%; concentrations 
much above 2 per cent do not presently appear 
to be attractive for metallurgical reasons. The 
absorption mean free path of a 2 per cent alloy 
of natural boron is about 0.15 cm. Since the 
natural element contains about 18.8 per cent 
B" (the highly absorbing isotope), this absorp- 
tion mean free path can be reduced by about a 
factor of 5, at considerable expense, by the use 
of highly enriched B'®. The main disadvantages 
of boron—stainless steel as a control-rod ma- 
terial are the relatively low percentage of boron 
which can be incorporated in the alloy and the 
alloy’s susceptibility to radiation damage which 
manifests itself as embrittlement of the alloy 
and some growth. 

The alloys may be fabricated by conventional 
means. Austenitic alloys containing up to 3.2 
wt.% boron have been cast by normal induction- 
melting methods in vacuum or inert atmos- 
pheres; alloy ingots containing up to 2.4 wt.% 
boron have been fabricated into strip by forg- 
ing and rolling at temperatures near 1160°C; 
alloys containing up to 2.2 wt.% boron have been 
extruded after heating to 1150°C. Both the 
sigma and the Heliarc welding techniques are 
applicable to the alloys, the sigma welding being 
preferred. 

The addition of boron to a stainless-steel 
alloy causes an increase in tensile strength and 
decreases in the ductility and impact strength. 
Neutron irradiation of the alloy causes further 
changes of the same kind. (Some irradiation 
tests on boron-—stainless steel were reviewed 
in the September 1959 issue of Power Reactor 
Technology, Vol. 2, No. 4.) Boron-10, when it 
absorbs a neutron, disintegrates into Li’ and 
an alpha particle. These products occupy more 
space than the original boron atom, and growth 
of the alloy occurs after sufficient irradiation. 
It is stated, however, that no growth appears to 
take place until nearly 0.7 per cent of the total 
atoms in the alloy have been transmuted by the 
(n,a@) reaction and that, after the irradiation 
has reached this point, growth is rapid and 
proceeds at an increasing rate. Metallographic 
examination of boron-—stainless steel alloys 
after irradiation has shown voids in the vicinity 
of boride particles; these voids were attributed 
to the helium gas (alpha particles) generated in 
the nuclear reaction. Other experiments which 
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Table VI-5 TENSILE-TEST RESULTS ON STAINLESS STEEL ALLOYED WITH 
NATURAL BORON (B%) AND ENRICHED BORON (B"")! 
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Nonirradiated Irradiated 
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, 0.2% Ultimate 0.2% Ultimate 
Boron Test yield tensile yield tensile 
e content, temp., strength, strength, Elongation, strength, strength, Elongation, Thermal 
y wt.% “¢ psi psi % psi psi Io nvt 
2 R 
t 0.0 (304 S.S.) 23° 37,500 90,000 50.2 82,000 104,800 2.7 x 10°° 
0.5 Bi? 23 29,500 81,500 29.6 98,800 121,600 2.4 Et <10" 
0.5 BY 315 23,200 61,900 19.8 82,600 92,500 0.5 6.7 x 107° 
2 0.85 BY. 23 32,950 85,000 21.1 111,800 122,800 0.3 8.0 x 107° 
> 0.85 B’” 315 23,800 65,700 10.3 90,900 0.1 8.0 x 107° 
. 1.18 BY 23 40,200 96,100 18.4 132,000 154,000 3.4 x 107° 
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1 1.64 BN 315 71,300 106,000 0.97 « 107° 
> 1.8 BN 23 41,300 84,200 8.0 136,000 151,000 5.8 x 107° 
‘ 1.8 BN 315 85,500 119,000 5.8 x 10°° 
2.04 BN 23 47,000 89,000 5.4 97,800 120,000 13x10" 
y 2.04 BN 315 71,300 106,000 7.3 x 107° 
2 were carried to exposures of 1 x 107! thermal Table VI-6 RESULTS OF IMPACT TESTS ON 
: : "a T N J NT + 
- nvt were reported to cause no gross dimensional era -ignpcgsereceenilitnypinetbent gt 1 
: BORON (B!°) ALLOYED WITH STAINLESS STEEL 
- changes or warping of the samples. weil . a 
A Tables VI-5 and VI-6 summarize some of the Nonives- has 
- data on tensile strength and impact strength Boron Test diated, diated, 
: contained in the reference.' Other experiments a | E.,* Izod Izod Thermal 
, P a rays 4 wt.% °C ft-lb ft-lb nvt 
n are cited in which samples containing various 
P ! amounts of natural boron up to 2 wt.% were ir- 347 S.S. 23 15.9 15.6 3.0 x 107° 
e | radiated to an exposure of 8 x 10° nvt. It is 0.5 a 23 9.5 1.48 3.0 x 10% 
: - 5 0 ¢ } 3.0> 20 
7 | stated that the yield strength increased by penne oe pes ‘= wep Steps 
; t 0.85 B 23 5.67 0.24 9.0 x 107 
250 per cent and the ultimate tensile strength 0.85 BY 315 5.70 0.18 9.0 x 102° 
1 increased by 80 per cent, whereas the elonga- 1.18 BN 23 2.3 0.61 4.35 x 10” 
: N ‘ 20 
d tion decreased to as low as 1 per cent of the 1.18 BN 315 1.8 0.66 = 4.35 x 10" 
original value. Since there appeared to be no eee ro ~" eda taticss* 
— et ee ee : 1.8 BY 315 1.0 0.35 1.94 x 10%” 
r difference in the irradiation behavior of the 2.04 BN 23 0.64 0.22 6.8 x 10° 
n alloys which could be attributed to the different 2.04 BN 315 0.64 0.05 6.8 x 107° 
d boron concentrations, it was postulated that the 
y tensile strength reached a saturation value at a : 
t reasonably low level of neutron exposure. 
d Marked embrittlement of alloys by irradiation Table VI-7 CORROSION BEHAVIOR OF SOME 
e was observed in other experiments at expo- STAINLESS STEEL— BORON ALLOYS IN 
19 600°F WATER 
h sures as low as 8 x 10°” nvt. 
" The use of boron-—stainless steel control e Cumulative weight gain, mg/cm? 
- oron 
D rods in water-cooled bad reactors has dem- Melt content, 156 514 1203 2013 3000 4008 De- 
1 onstrated that the corrosion resistance in high- No. wt.% hr hr hr’ hr- hr_ hr scaled 
e temperature water is acceptable. Experimental 
n measurements of corrosion rates are not en- 2 O38 BE 2D) 1D) SH BA 08) 88 
d i : SS-2 O75. 17 23. 13, 32 4S. 33..:=49 
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Cc ditions of boron in the 1 to 2 wt.% range to ss-4 $66 2S 212 32-38 38 -292 88 
S stainless-steel alloys do not affect greatly N-272 2.32 2.6 2.8 21 44 3.4 2.0 —-9 
y the corrosion resistance in dynamic tests and Water resistiv- 
d that irradiation does not have an important ef- ity at end of 
7 . test period, 
n n. Som tic tests, however 
fect oh corrosto nervinian wl om mee " megohm-cm 0.8 0.3 0.3 0.2 0.2 0.2 
h have shown much greater weight gains for a 
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boron—304 stainless steel alloy than for an 
ordinary 304 stainless steel alloy and a much 
greater rate of weight gain after irradiation 
than before irradiation for the boron alloy. The 
latter difference was attributed to stress cor- 
rosion induced by the helium formed in the 
(n,@) reaction. The results of corrosion tests 
on unirradiated boron—stainless steel alloys at 
KAPL are given in Table VI-7. Presumably 
these were static tests. 


Other Boron-containing Materials! 


Alloys of boron with zirconium and titanium 
have also been made. The addition of boron to 
zirconium and titanium has effects qualitatively 
similar to those resulting from its addition to 
stainless steel: the tensile strength is increased 
and ductility is decreased; neutron irradiation 
produces a further very marked decrease in 
ductility. Boron contents between 1 and 2 wt.% 
appear to be practical. The unirradiated alloys 
with stainless steel, titanium, and zirconium 
show corrosion resistance to high-temperature 
water which decreases in the order named. 
The postirradiation corrosion-resistance trend 
is similar but more pronounced. Zirconium- 
base alloys irradiated to relatively high boron 
burnup appear to have poor corrosion resist- 
ance. 

Of the boron compounds, boron carbide is the 
most widely used. It is a refractory material 
with a theoretical density of 2.51 g/cm’; actual 
densities up to 2.4 g/cm® can be obtained in 
massive bodies. The material, although brittle, 
has relatively good thermal-shock resistance 
because of its high compressive strength 
(414,000 psi at room temperature), its rela- 
tively small value of thermal-expansion coef- 
ficient (45 x 107"/°C), and its reasonably high 
thermal conductivity (0.24 and 0.15 watts/ 
(cm)(°C) for theoretically dense material at 
200 and 800°C, respectively). The corrosion 
resistance of the material for high-temperature 
water is not considered good but is sufficiently 
good to prevent catastrophic corrosion rates if 
jacketed material is exposed to water through 
failure of the jacket. The corrosion resistance 
is much better for high-density material than 
for low-density material. 

Neutron irradiation of boron carbide produces 
anisotropic dimensional changes in the boron 
carbide crystals. These effects appear to an- 


neal out at 700 to 900°C. However, the deposi- 
tion of helium in the lattice is a further cause 
of radiation damage. Samples irradiated to a 
burnup of 36 per cent of the B’* atoms were 
reduced to powder. The. effects of elevated 
temperatures on radiation damage are not well 
known, and it is thought that high temperatures 
may be beneficial. At low temperatures, only 
a small fraction of the helium is released from 
the material even after the irradiation has 
proceeded to the point where the material 
breaks up into granules. Subsequent heating of 
such a sample to 800°C caused the release of 
about 19 per cent of the contained helium. 

The borides of a number of metals, including 
titanium, hafnium, zirconium, niobium, and 
tantalum, have received some study, particu- 
larly with regard to aqueous corrosion, radia- 
tion damage, and helium retention. As yet, 
none of these materials has seen service as a 
control-rod material. * 

Dispersions of boron or boron compounds in 
ductile metals may be used as alternatives to 
the boron alloys or to canned boron compounds. 
The helium problem is still present in the dis- 
persion; it is usually the metal of the dispersion 
which is counted on to contain the helium. The 
most common material of this type is Boral, a 
dispersion of boron carbide in aluminum. This 
material is, of course, limited to low-tempera- 
ture applications. The absorbing material for 
the first set of APPR control rods consisted of 
a dispersion of 3.23 wt.% boron (90 per cent 
enriched) in iron, clad with type 304 stainless 
steel. After approximately two years of serv- 
ice, this material has been replaced by a Eu,0;— 
stainless steel dispersion that is stainless- 
steel-clad. The B'-Fe dispersion rods are 
being given postirradiation metallurgical ex- 
amination. 

Boron dispersions in stainless steel, tita- 
nium, f and Zircaloy have been investigated. 
The titanium-base dispersions exhibit excellent 
corrosion resistance in hot water when the 
boron content is relatively low, up to perhaps 
3 per cent. Titanium-clad specimens of the 
titanium-boron matrix, purposely defected, ap- 
peared to be good up to boron contents of at 
least 12.9 wt.%. The corrosion does not appear 





* Reference 18 discusses briefly some of the borides 
of interest. 

{Two recent papers on boron-titanium dispersions 
are listed as references 19 and 20. 
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to be greatly increased by irradiation. The 
stainless steel—boron dispersions appear to be 
reasonably good when the boron content is low 
and rather poor when the boron content is high. 
A sample containing 1.6 wt.% boron lost 38 mg/ 
dm’ in 70 days at 345°C, whereas a 6 wt.% B" 
dispersion lost 1250 mg/dm? under the same 
conditions. : 

The addition of boron in dispersion lowers 
the strength of the base material at room tem- 
perature, although at higher temperatures its 
effect may be small. Irradiation embrittles the 
dispersion, but the degree of embrittlement 
seems to vary with the base material. Titanium 
appears to be the least affected, whereas steel 
is quite sensitive in thisrespect. Longirradia- 
tion produces dimensional changes of a few 
per cent and may produce blistering, apparently 
as a result of the helium accumulation. Blister- 
ing has been observed in a titanium-B’ dis- 
persion containing 1 per cent B” after the ex- 
posure had reached 40 per cent burnup. 

Dispersions of boron carbide in titanium and 
zirconium have shown the same pattern of 
aqueous corrosion behavior as have the boron 
dispersions. In general, the susceptibility to 
corrosion increases with the amount of boron 
carbide; the titanium dispersions can, on the 
whole, be described as having good corrosion 
resistance, but the zirconium and Zircaloy dis- 
persions are rather poor. 

The addition of boron carbide as a dispersion 
increases the brittleness. In the case of tita- 
nium-base dispersions, this increase was found 
to be greater than that produced by an equal 
quantity of pure boron in dispersion. The pri- 
mary effects of neutron irradiation, again, are 
embrittlement and some increase of volume. 

Some work has been done on dispersions of 
the diborides of zirconium, titanium, and haf- 
nium, in zirconium and titanium matrices. The 
combination zirconium diboride in a titanium 
matrix presently appears to be most promising. 
Again the titanium base exhibits the best resist- 
ance to aqueous corrosion. 

The effects of irradiation on a number of 
dispersions containing boron or boron com- 
pounds are summarized in Table VI-8. 


Cadmium and Cadmium Alloys! 


Metallic cadmium has beenused rather widely 
as a control-rod material in low-temperature 





reactors. Its low melting point (321°C) pre- 
cludes its use for high-temperature applica- 
tions. It is a soft metal, of low tensile strength 
(2500 to 9000 psi ultimate tensile strength for 
as-cast specimens); it is ordinarily used with 
a jacket or cladding to provide mechanical 
strength. It does not have good corrosion re- 
sistance to water but does not react rapidly 
with water at temperatures up to the boiling 
point. 


A ternary alloy of silver, indium, and cad- 
mium has been developed, primarily at the 
Bettis Plant of the Westinghouse Electric Cor- 
poration, for use in pressurized-water reactors. 
Previous work had been done on alloys of silver 
and cadmium, but these binary alloys are in- 
ferior to the ternary alloy in both corrosion 
resistance and neutron absorption. 


The alloy composition which appears to be 
nearly optimum for corrosion resistance is 
81.5 Ag—13.6 In—4.9 Cd, where the composi- 
tion is given in weight per cent. This alloy in 
water showed a corrosion rate of about 0.18 
mil/year. The rate did not vary much with 
temperature over the range tested, up to 288°C, 
in degassed and hydrogenated water, nor was 
the rate affected much by a change from neutral 
to high PH water or by a change from static to 
dynamic test conditions. However, in oxygen- 
ated water (8 to 12 cm*/kg) at 288°C and high 
pH, the corrosion rate was high, amounting to 
6 mg/cm’ at the end of 28 days. 


The absorption of neutrons converts silver 
into cadmium, converts a little of the cadmium 
into indium, and converts indium into tin. Con- 
sequently, after long irradiation the composition 
of the alloy will change. An alloy 76.5 Ag— 
12.1 In—9.4 Cd-—2.0 Sn, which simulated the 
one obtained by irradiation of the ternary alloy 
for 3000 hr in a flux of 1 x 10" neutrons/(cm’) 
(sec), was also corrosion tested and was found 
to have better corrosion resistance than the 
ternary alloy. 

The neutron-absorption effectiveness of the 
ternary alloy approaches that of hafnium. The 
quaternary alloy corresponding tothe irradiated 
material (alloy number 2115) is somewhat less 
effective. Table VI-9 shows the relative ef- 
fectivenesses of various alloys as determined 
experimentally in a critical facility which was 
said to be designed to simulate the neutron-flux 
spectrum expected in full-scale pressurized- 
water reactors. 
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occur in approximately equal abundance, and 
each has an appreciable thermal-neutron cross 
section. The daughters produced by neutron 
capture also have high cross sections; conse- 
quently, europium is attractive for rods which 
must absorb many neutrons over their life- 
time. Europium is one of the less abundant and 
more costly of the rare earths. Two of the 
gadolinium isotopes, comprising about 30 per 
cent of the natural material, have high cross 
sections. However, there is no continuous se- 
ries of high-cross-section daughter decay prod- 
ucts. Samarium has only one isotope of high 
cross section, of 13.84 per cent abundance. 
Consequently, samarium will burn out rapidly 
unless used in relatively high concentrations. 
Dysprosium has a lower cross section than the 
other three rare earths of interest and has been 
considered primarily for use as a burnable 
poison. 

The rare earths, in the metallic form, are 
reactive chemically and must be protectedfrom 
most reactor coolants. Work is under way to 
produce alloys with corrosion-resistant ma- 
terials, but it is probable that such alloys will 
need protective cladding in most applications. 
Alloys with titanium and stainless steel have 
been made. In general, increasing amounts of 
the rare earths cause an increase in brittle- 
ness of the rolled alloy material. The corro- 
sion behavior of these alloys appears very 
promising for rare-earth contents of a few 
per cent, but corrosion resistance deteriorates 
as the percentage of rare earth is increased. 
Little or no work has yet been done on the ef- 
fects of irradiation on these alloys. 

The rare-earth oxides hydrate rapidly when 
exposed to hot water; and, consequently, if the 
rare earths are used in oxide form, these too 
must be protected. Most of the present work is 
directed toward the use of dispersions of the 
rare-earth oxides in a metal matrix; this cermet 
is then clad with an appropriate material to 
eliminate corrosion. There is an economic 
incentive for using the rare earths in the oxide 
form rather than the metallic form since the 
cost of producing rare-earth metals is rela- 
tively high. Among the metals which have been 
used in rare-earth oxide dispersions are stain- 
less steel, titanium, zirconium, and nickel. 
Much of the work has been done with samarium 
oxide rather than with the more expensive 
europium oxide, but the behavior is expected 
to be nearly the same for the two materials. 


Stainless steel—Sm,O; dispersions, which are 
stainless-steel-clad and which contain from 15 
to 40 vol.% of the oxide, have been successfully 
fabricated. Since the rare earths capture neu- 
trons by an (n,y) reaction, no very strong ir- 
radiation-damage effects would be expected. 
However, bend tests on irradiated Eu,O;—stain- 
less steel cermets showed some loss of ductility 
during irradiation. It is stated that examina- 
tion of micrographs showed that the damage 
was adjacent to the Eu,O; particles and that the 
cause of this damage was later shown to result 
from the reaction between the silicon in the 
stainless steel and the Eu,O3. Impact tests on 
9.9 wt.% Eu,O, dispersions in stainless steel 
after two cycles of MTR irradiation showed a 
decrease in impact strength by only 5 to 10 per 
cent in the temperature range 24 to 315°C. 
Bend tests on irradiated Gd,O, cermets, ir- 
radiated for two cycles in the MTR, showed 
that the material retained 90 per cent of its 
preirradiation ductility. However, growth was 
noted of the order 0.3 per cent in all three di- 
mensions. 


Table VI-12_ CONTROL MATERIAL Costs? 





Cost per pound, 





Material dollars Basis of cost 

Boron (natural) 450 (99.6% Supplier 

purity) 
15 (commercial 
grade) 

B!° (93% enriched) 1500 Private 
communi- 
cation 

B,C 5-10 Supplier 

Boral 8-18 (per Supplier 

square foot of 
,- to %-in.- 
thick sheet) 

Boron- stainless steels: 

2 wt.% BN 6.35 Supplier 
1 wt.% B® (90% 52.40 (500-Ib Supplier 
enriched) lot) 

Cadmium 2-3 (metal) Supplier 

Hafnium 50-100 (crystal AEC esti- 

bar) mate 

Dysprosium oxide 50* (99% purity) Supplier 

Europium oxide 790 (99% purity) Supplier 

Europium metal 1500 Supplier 

Gadolinium oxide 35* (90% purity) Supplier 

Samarium oxide 25* (90% purity) Supplier 


Titanium sponge powder 
Stainless-steel powder 
Low-cobalt stainless 
steel equiv. to AISI 
304 E/C (ingots) 
Cd-In-Ag billets 


6 
20 


5 
20 


Reference 23 
Reference 23 


Reference 23 
Reference 23 





* For lots of over 500 lb. 
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Material Costs and Applications} 


Although the cost of a complete control rod 
involves many items besides the cost of the 
absorbing material in the rod and may be, in 
some cases, many times the cost of that ma- 
terial, some of the materials are rather ex- 
pensive; their cost will, therefore, play a part 
in determining whether or not they will be at- 
tractive for a given application. The costs of 
some of the absorbing materials which have 
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been discussed above are listed in Table VI-12, 
from reference 1. 

Manufacturing experience sufficient to per- 
mit estimates of the relative costs of complete 
control rods of the various types has not yet 
been accumulated. In the EBWR, where some 
of the rods were of hafnium and some were of 
boron—stainless steel, a relatively clean com- 
parison of manufacturing costs is possible. The 
total fabrication cost for a single hafnium 
cruciform rod, '-in. thick by 10-in. span by 


Table VI-13_ CONTROL-ROD TYPES IN VARIOUS POWER REACTORS! 





Reactor Control-rod environment 


Control-rod type 





Calder Hall CO, 


outlet temperature 637°F 


EBWR (Experimental Boiling- 
water Reactor) 


Boiling H,O 
outlet temperature 488°F 


VBWR (Vallecitos Boiling- 
water Reactor) 


Boiling H,O 
outlet temperature 545°F 


ALPR (Argonne Low Power Boiling H,O 


Reactor) outlet temperature 420°F 
Dresden Nuclear Power Boiling H,O 
Station outlet temperature 542°F 


Sodium 
outlet temperature 800°F 


Enrico Fermi Atomic 
Power Plant 
(fast-neutron reactor) 


OMRE (Organic Moderated 
Reactor Experiment) 


Santowax OM 
cutlet temperature 


APPR (Army Package Power Pressurized H,0 


Reactor) outlet temperature 450°i" 
Shippingport Atomic Pressurized H,O 

Power Station outlet temperature 539°F 
Consolidated Edison Pressurized H,O 

(Indian Point) reactor outlet temperature 510°F 
Yankee Atomic Electric Pressurized H,O 

Plant outlet temperature 514°F 
STR (Submarine Thermal Pressurized water 

Reactor) 
ARE (Aircraft Reactor Helium rod coolant, 

Experiment) separated from the 1500°F 


molten-salt fuel mixture 
Helium rod coolant, 
separated from the 960°F, 
outlet temperature, 
sodium reactor coolant 


SRE (Sodium Reactor 
Experiment) 


530-710°F 


Number of rods adjustable, approximately 40; boron 
steel cylinders encased in low-cobalt stainless- 
steel sheaths 

Nine cruciform rods, 10-in. span by 46 in. long 
(five rods are '4-in. hafnium with Zircaloy fol- 
lowers; four rods are /,-in., 2 per cent natural 
boron-—stainless steel with Zircaloy followers) 

Three central blades, 18 in. wide by '/, in. thick, of 
1 wt.% B'° stainless steel; four outer blades con- 
tain '/,-in.-OD cylindrical steel tubes filled with 
B,C 

Five cruciform rods, 14-in. span, 32 in. long, 

0.220 in. thick, composed of 0.060-in.-thick 
cadmium canned in 0.080-in.-thick aluminum- 
nickel alloy; followers are of aluminum-nickel; 
four T-shaped rods of the same composition are 
provided for larger fuel loadings 

Eighty cruciform rods; 2 wt.% natural boron-— 
stainless steel, 0.375 in. by 6'/,-in. span by 
102 in. 

Ten cylindrical rods, enriched B,C hollow cylinders 
encased in stainless-steel tubes; the stainless- 
steel tubes act as pressure vessels to contain the 
helium generated 

Twelve eylindrical rods, composed of stainless-steel 
tubes, 1/, in. OD by 36 in. long, containing com- 
pacted B,C 

Seven hollow square tubes, 2, in. square, with fueled 
followers; absorbing material is 3.23 per cent B!” 
in iron, clad with stainless steel; these original 
rods were replaced by others using Eu,O;—stain- 
less steel dispersion in the spring of 1959 

Thirty-two cruciform hafnium rods, 0.226 in. by 
3%%-in. span by 71 in. long 

Twenty-one cruciform hafnium rods, ¢ in. by 7.5-in. 
span by 135% in. long 

Thirty-two extruded cruciform rods, of 80 wt.% 
silver, 15 wt.% indium, and 5 wt.% cadmium 

Hafnium control rods 


One regulating rod of stainless steel; three shim rods 
of segmented cylindrical B,C canned in stainless 
steel 

Eight cylindrical rods of 2.5 wt.% boron—nickel alloy 
rings encased in stainless-steel tubes; the shim 
rods use 18 rings, 2.5 in. OD by 1.75 in. ID by 4 in. 
long; the safety rods are slightly smaller 














46 GENERAL RESEARCH AND DEVELOPMENT 


46 in. long, with a 39-in.-long Zircaloy-2 fol- 
lower, has been given as $9500. The total 
fabrication cost for the boron—stainless steel 
rod, containing 2 wt.% natural boron, and simi- 
lar to the hafnium rod in dimensions except for 
an increased thickness to '/, in., has been stated 
as $2750. Of this cost, it is estimated that ap- 
proximately $1900 is the cost of the Zircaloy-2 
follower. The fabrication cost of control rods 
does not, of course, tell the complete story; 
for some rods, because of better burnout char- 
acteristics, better resistance to radiation dam- 
age, and better corrosion resistance, may last 
longer in the reactor than other rods. This 
consideration is reflected in the choice of the 
EBWR rods. The hafnium rods, which are 
superior with respect to burnout and radiation- 
damage resistance, are used to compensate the 
excess reactivity provided for fuel burnup; 
whereas the boron-steel rods are used pri- 
marily for shutdown and, on the average, re- 
ceive considerably less neutron exposure. Thus 
the boron-steel rods, although much cheaper 
than the hafnium rods, do not do an equivalent 
job. 

Current practice in the choice of absorbing 
materials for control rods is summarized in 
Table VI-13, which lists the materials for a 
number of power reactors which are either in 
operation or under construction. Reference 1, 
the source of the information in Table VI-13, 
contains a much more extensive list of the 
same kind. 
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DIRECT CONVERSION OF NUCLEAR ENERGY 





The operating principles and status of develop- 
ment of direct energy-conversion schemes were 
reviewed in a series of papers presented at the 
American Nuclear Society Symposium held in 
Pittsburgh on Mar. 19, 1960. The basic and 
practical limitations on obtaining increases in 
efficiency and power density were discussed. 
Conceptual methods of integrating direct-con- 
version devices with reactor heat sources, and 
some of the foreseeable major problem areas, 
were also presented. 

The direct-conversion schemes described 
are: (1) thermionic emitter, (2) thermoelectric, 
(3) magnetohydrodynamic (MHD) generator, and 
(4) electrochemical energy converter or fuel 
cell. 

To attain high efficiencies in the direct- 
conversion heat engines, it is necessary to re- 
duce the irreversible heat losses. The heat 
losses are due to straight thermal conduction 
and radiation. The heat leak paths are gen- 
erally the same as the electron flow paths. 
Hence it is difficult to reduce heat leaks without 
increasing the /’R drop internally and/or in the 
lead wires. 

Dr. C. E. Zener of the Westinghouse Electric 
Corporation gave a measure Of the irreversible 
heat losses occurring in presently developed 
direct-conversion devices. He expressed this 
measure as follows: 


n = 6, 


where 7 = over-all efficiency 


_ electrical power output from device 
heat input into device 





N- is the Carnot efficiency, and 5 is the degra- 
dation coefficient. 


5 
Regenerative H, fuel cell 0.10 
Thermoelectric (bismuth telluride) 0.17 
Thermionic emitter (2 mm spacing, 
cesium-coated anode) 0.30 
MHD >0.85 
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For the MHD generator, the lead losses result 
in only a small degradation of energy. 

By proper design, it should be possible to 
improve the degradation coefficient 5. The goal 
is to uncouple the electron flow path from the 
heat leak path. 

With the thermionic emitter, high power 
densities (about 30 watts per square centimeter 
of cathode surface) have already been produced. 
Laboratory models have been built and oper- 
ated at 13 to 15 per cent efficiency. The process 
does require high-temperature operation (1500 
to 2800°K) which presents materials problems. 
J. R. Reitz of Case Institute of Technology 
described the several regimes of operation of 
a thermionic emitter, namely, (1) vacuum, (2) 
low-density plasma, and (3) high-density plasma. 
For high current densities (of the order 100 
amp/cm?’), it is necessary to operate in the 
high-density plasma regime. The Los Alamos 
plasma thermocouple operated in this regime 
to deliver about 60 amp/cm? at 2500°K cathode 
temperature. 

S. J. Angello of the Westinghouse Electric 
Corporation presented a brief review of the 
thermoelectric effects. From an elementary 
derivation, he obtained 


7 2 
= 2G 
4 \pk 
where T) = hot-junction temperature 
a@ = thermoelectric power 


p = electrical resistivity 
k = thermal conductivity 


Since a, p, and k are basic material properties, 
it is obvious that material development is 
crucial in increasing the efficiency of thermo- 
electric devices. The thermoelectric materials 
must withstand high temperatures and, when 
used in a reactor environment, must also with- 
stand radiation damage. Present type and p 
type materials are available which permit 
operation up to 1200°C to give an over-all 
theoretical efficiency of 16 per cent. Attention 


bE 
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must be paid to extraneous losses in the design 
of the thermoelectric device and in integration 
of the device with the heat source. Heat leak 
paths not directly through the elements them- 
selves must be minimized. Care must be taken 
to minimize contact resistance to current flow 
between the elements and the heat-collector 
and heat-sink skins. 

A great deal of material research work is 
required to develop thermoelectrics capable of 
higher temperature operation and having ma- 
terial properties leading to higher efficiencies. 
Between 1956 and 1960 the efficiency has been 
increased from 6 to 16 per cent. It is projected 
that with sustained effort on material develop- 
ment an efficiency of 30 per cent canbe attained 
by the year 1970. 

R. J. Rosa of Avco-Everett Research Labo- 
ratory described the MHD generator. The MHD 
generator extracts power directly from the 
motion of ionized gas through a magnetic field. 
It consists of a duct through which an ionized 
gas flows, coils which produce a magnetic field 
across the duct, and electrodes at the top and 
bottom of the duct. The ionized gas, by virtue 
of its motion through the magnetic field, has 
an electromotive force generated in it which 
drives a current through it, the electrodes, and 
the external load. The generator has no moving 
parts or close tolerances, and hence its cost is 
expected to be low compared to that of conven- 
tional conversion machinery for the same power 
output. Its efficiency is comparable to that of 
a gas turbine, i.e., it will have approximately 
the same enthalpy extraction from the gas for a 
given pressure drop through the device. 

The most important quantity determining the 
feasibility and characteristics of MHD gen- 
erators is the gas conductivity. For air seeded 
with 1 per cent potassium vapor, sufficient 
thermal ionization (~1 mho/meter gas con- 
ductivity) is obtained at temperatures in excess 
of 2000°K for a reasonably practical design. 
With this value of gas conductivity, a 100-Mw- 
output MHD generator would be about 30 ft in 
length and would require a 10,000-gauss mag- 
netic field strength. It is evident that MHD 
generators tend to be quite large. 

R. J. Rosa described an MHD device having 
a 10-kw power output which was successfully 
operated as a demonstration of the basic work- 
ing principles of the device. The MHD gen- 
erator extracted about 1'/, to 2 per cent of the 
enthalpy of the gas. 


Schemes for irtegrating the direct-conversion 
devices with a nuclear heat source were dis- 
cussed conceptually, principally with a view of 
attempting to define the problems involved. 
The application of immediate interest is for 
space-satellite auxiliary power supplies wherein 
an important requirement is maintenance-free 
operation at full power for time periods of a 
year or more. The power requirements are 
quite modest, ranging from about 1 to 20 kw of 
net electrical power output. For this power 
range the reactor power density is low enough 
that it may prove practical to conduct the re- 
actor heat to the surface of the reactor. Obvi- 
ously, to keep the temperature at the center of 
the reactor within reasonable bounds and to 
minimize thermal stresses in core and re- 
flector, structural material of high conductivity 
is required. The heat-collector skin of the 
thermoelectric device can be wrapped around 
the reactor to receive its heat by conduction 
and/or radiation from the reactor surface. The 
heat-sink skin radiates the waste heat to the 
space environment. 

For the thermionic emitter application, the 
cathode surface is wrapped against the cylin- 
drical surface of the reactor. The anode metal 
surface is also cylindrical and is concentric 
with the cathode surface; the intervening space 
is filled with cesium vapor. Again, the waste 
heat from the anode is rejected by radiation. 

At the higher reactor power densities re- 
quired for central power stations, the reactor 
heat could be removed in the conventional 
manner by forced circulation of a coolant (high- 
pressure gas or liquid metal) through the re- 
actor. The high-temperature coolant would then 
be pumped past the heat-collector surface of 
the direct-conversion device. In this concept, 
the reactor heat source and the direct- 
conversion device are physically separated. 
The high temperatures required pose serious 
material problems and present a need for new 
reactor concepts. 

Considerations are given to placing the 
thermoelectric material in contact with the fuel 
rods so as to receive the heat directly. A 
similar scheme for using the thermionic emitter 
as a topping device for nuclear reactors has 
been proposed and has been previously de- 
scribed. (See the section on Thermionic Energy 
Conversion in the December 1959 issue of 
Power Reactor Technology, Vol. 3, No. 1.) An- 
other scheme is to place the thermionic emitter 
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and the thermoelectric device in series within 
the reactor. In this scheme the heat from the 
fuel rod first passes through the thermionic 
emitter. The cathode surface is the inner wall 
of a double-walled jacket for the fuel rod. The 
cathode is heated directly by fission in the fuel. 
The annular space within the double-wailed 
jacket is filled with cesium vapor. The outer 
wall is the anode surface of the thermionic 
emitter and is also the heat collector for the 
thermoelectric elements. The heat rejected by 
the thermionic emitter passes through the 
thermoelectric elements to the heat-sink skin. 
The heat-sink skin is cooled in the conventional 
manner by the reactor coolant. With this scheme 
the thermoelectric device works as a topper in 
series with the conventional nuclear energy 
conversion system and, in turn, the thermionic 
emitter is a topper onto the thermoelectric 
device. 

The problem of integrating the reactor heat 
source with the MHD generator appears to re- 
quire development of unusual reactor concepts, 
such as the nuclear cavity reactor, in order to 


achieve the required gas temperatures without 
melting down the structures of the reactor. 

J. C. Danko of Westinghouse Electric Corpo- 
ration presented some results obtained on ir- 
radiation effects on thermoelectric materials. 
The prime subject for investigation was 
Li,Ni;_,O. In-pile measurements of the Seebeck 
coefficient and electrical resistivity were made 
at reactor ambient temperatures and elevated 
temperatures. In addition, pre- and postirradia- 
tion measurements of the Seebeck coefficient, 
the electrical resistivity, and the thermal con- 
ductivity were performed. 

Postirradiation measurements on samples of 
Lip, o143Nip.g70 after exposure toa total integrated 
thermal-neutron flux of 1.4 x 10'* revealed that 
the electrical resistivity had returned to the 
preirradiation value; the Seebeck coefficient 
had increased slightly; and the thermal con- 
ductivity had remained essentially unchanged. 
In-pile measurements of the Seebeck coefficient 
and electrical resistivity of Lig 9;Nig,,;0 at 
464°C and up to 8x 10 nvt (thermal) did not 
show much variation from the out-of-pile values. 





PROGRESS ON SPECIFIC REACTOR TYPES 


DESIGN STUDIES AND EVALUATIONS 





A considerable volume of literature covering 
design and evaluation studies of specific reac- 
tor types, prepared by or sponsored by the 
Atomic Energy Commission, has been issued 
within the last several months. It is the pur- 
pose of this note to list these reports and to 
clarify their relations to one another. 

In 1958 the AEC was authorized, under Public 
Law 85-590 (Sec. 101), to perform a series of 
design and engineering studies of four large 
power reactors. The results of those studies 
have been published as references 1to5. These 
reports were reviewed in the September 1959 
issue of Power Reactor Technology, Vol. 2, No. 
4, shortly after their publication. Since then, 
four more reports (references 6 to 9) have 
been issued which contain additional studies 
and information supplementing the original re- 
ports. 

In June 1959 the Division of Reactor De- 
velopment of the AEC undertook the task of 
defining a long-range program for civilian 
power-reactor development, generally referred 
to as the “Ten-year Program,”’ and a further 
group Of evaluation studies was made as part 
of that project. The large-power-reactor studies 
in references 1 to 9, although they undoubtedly 
contributed to the evaluations, did not con- 
stitute a formal part of the “Ten-year Pro- 
gram’’ study. 

At the November meeting of the Atomic In- 
dustrial Forum, Frank K. Pittman, Director of 
the Division of Reactor Development, sum- 
marized the results of the studies to evaluate 
the potentials of the several reactor types. His 
paper was reviewed in the March 1960 issue of 
Power Reactor Technology, Vol. 3, No. 2. Since 
then, three reports covering the studies and the 
proposed development program (references 10 
to 12) have been issued. A fourth report (ref- 
erence 13) in several parts, covering the tech- 


nical status studies for the individual reactor 
types, is yet to be published. The three re- 
ports issued to date have been summarized in 
reference 14. 
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THE PLUTONIUM RECYCLE TEST REACTOR 





The Plutonium Recycle Test Reactor (PRTR) is 
an AEC facility, under construction at Hanford, 
which will be used to develop the technology re- 
quired for the utilization of plutonium fuels in 
thermal heterogeneous power reactors. Specifi- 
cally, the reactor provides facilities for the 
following:' 


1. Irradiation testing of plutonium-bearing 
fuel elements. 

2. Direct investigations of reactivity and ex- 
posure effects from isotope buildup on the 
uranium-plutonium fuel cycle. 

3. Production of irradiated fuels for reproc- 
essing and refabrication studies. 

4. Investigation of reactor operating prob- 
lems for plutonium-recycle operation. 

5. Demonstration of the economics and prac- 
ticability of various fuel cycles. 


Operationally the facility will be utilized for 
research and development to test the feasibility 
of various plutonium-recycle concepts. Empha- 
sis will be placed on obtaining physics data 
necessary to evaluate plutonium recycle, e.g., 
data pertaining to isotope cross sections and 
data on control characteristics. 

The reactor is of the vertical pressure-tube 
type, moderated, cooled, and reflected with 
heavy water, and is designed to operate at a 
maximum of 70 Mw(t). The moderator is cor- 
tained in an unpressurized aluminum calandria 
tank that has an 84-in. ID and is 115 in. in 
height. Eighty-five fuel channels, 18 shim 
control-rod channels, and 13 instrumentation 
channels are mounted vertically within the 
calandria with the fuel channels arranged on an 
8-in. triangular pitch. To insulate the coolant 
from the moderator, a double-wall construction 
is utilized for the fuel channels. The outer 
aluminum tube is an integral part of the calan- 
dria; the inner process tube is fabricated of 
Zircaloy-2 and contains the fuel elements. Each 
process tube is connected to coolant manifolds 
by means of jumpers to allow maximum flexi- 
bility in monitoring and cooling individual as- 
semblies. Heavy water flows vertically upward 


within the process tubes at a pressure of about 
1000 psi. Heat generated in the reactor is used 
to make 425-psia steam in a primary heat ex- 
changer, and this steam is condensed and dis- 
posed to the Columbia River. The heavy-water 
reflector is contained in an annular cylinder 
coaxial with the calandria. 

The primary reactivity control system of the 
reactor is based on controlling the level of the 
heavy-water moderator within the calandria. A 
schematic of the primary control system is 
shown in Fig. 13. The moderator level in the 
calandria is maintained by means of a helium- 
gas-pressure balance system. Heavy water is 
continuously withdrawn from the storage tank 
and pumped into the calandria via a moderator 
heat exchanger. When the moderator level is in 
the normal operating range, heavy water is re- 
turned to the moderator storage tank from the 
calandria by flow through discharge ports atthe 
top of the calandria and by flow over a weir at 
the bottom of the calandria. The pressure head 
required to maintain a given level is developed 
by means of a helium compressor; at steady- 
state operation most of the gas pumped through 
the helium compressor is returned to the com- 
pressor inlet through the control valves, but a 
small amount is sent to the helium purification 
system and is used to sweep the upper part of 
the calandria. The position of the control valves 
is determined by an automatic controller, a 
manual valve operator, or a safety interlock. 
This system is stated to be capable of main- 
taining the moderator within +0.05 in. over a 
range of moderator levels from 36 to 111 in. 
from the bottom of the calandria. Scram of the 
reactor is accomplished by rapid drainage of 
the moderator from the calandria. Uponreceipt 
of the scram signal, the four 8-in. dump valves 
open, equalizing gas pressure between the top of 
calandria and the weir. The moderator then 
flows over the weir and into the dump chamber 
immediately surrounding the weir. Drainage 
from the dump chamber into the storage tank is 
through five 8-in. pipes; the dump chamber thus 
acts as a “receiver” for a significant part of the 
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Figure 13—Primary-control-system flow diagram. 


moderator during a scram, and this flow takes 
place rapidly. It is stated that, within 0.85 sec 
after receipt of a scram signal, the moderator 
level has fallen 2 ft, decreasing reactivity by no 
less than 0.018 kos (for the moderator level 
initially at its maximum). 

Gross adjustments of reactivity are provided 
by the shim controls. Eighteen multiple-shim 
control units drive a total of 54 Inconel shim 
rods, worth a total of 0.115 kers. The rods are 
cooled by the heavy-water moderator and are 
top driven; when withdrawn from the core, they 
are stored in tubes extending down from the 
calandria into the bottom shield. The shim rods 
are operated manually by switches located 
within the control room. The shim rods are not 
intended as safety devices and are not actuated 
by scram signals. 

The PRTR initial loading will contain two 
different types of fuel elements. The Mark-I 
element consists of a 19-rod cluster assembly 
employing a plutonium-aluminum fuel alloy or 
sintered UO, pellets. The fuel material is clad 
with Zircaloy-2 tubing, and 12 ofthe 19rods are 
wrapped with a 72-mil Zircaloy-2 wire to main- 
tain alignment of the tubes and thus provide for 


adequate flow of coolant. The Mark-II elementis 
composed of a rod of UO, surrounded by two 
concentric annular rings of UO,; Zircaloy-2 is 
also used as a cladding. The annular coolant- 
flow passages are maintained by means of ribs 
on the surface of the clads. Other plutonium fuel 
elements are also being considered, such as 
uranium dioxide—plutonium dioxide and ele- 
ments wherein the active material is carried in 
a ceramic matrix. 

The fuel-handling system provides for the 
charging and discharging of fuel elements when 
the reactor is shut down. The basic component 
of the system is the charge-discharge machine. 
This machine can be indexed over any process 
tube by means of rotating shield plugs, and a 
connection can be made to the process tube, 
after removing a shielding plug in the floor of 
the rotating shield. A hook is then lowered into 
the process tube to engage the fuel element, 
which is removed and storedina shielded region 
of the charge-discharge machine. The machine 
is designed to handle the process tubes, should 
this ever be required. Underwater storage of 
irradiated elements is also provided in a pool 
exterior to the reactor containment building. 


AV 


a VD ew 


Table IX-1 PLUTONIUM RECYCLE TEST REACTOR (PRTR) #NGINEERING DATA! 
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General Reactor Data 


Power level (nominal) 
Heat to primary coolant loop 
Heat to moderator coolant loop 
Heat to reflector 
Heat to shields 
Primary coolant flow 


Primary coolant pressure 


Reactor inlet temp. 
Reactor outlet temp. (bulk) 
No. of process tubes 
Lattice 


Fuel loading: 

No. of UO, tubes 

No. of Pu tubes 
No. of effective power tubes 
No. of effective flow tubes 
Maximum tube power 
Maximum tube flow 
Maximum tube outlet temp. 
Av. neutron flux (thermal) 


Goal fuel exposure, U 

Xenon override time 

Time required for: 
Shutdown (normal) 

Cool and depressurize primary 
coolant system (normal) 
Pressurize primary coolant 

system startup (from 
critical to full power) 
Cold 
Hot 
Power and neutron-flux distribu- 
tion (unflattened): 
Radial flux distribution 


Vertical flux distribution 
(7.4-ft fuel) 


70 Mw 

66.4 Mw 

2.2 Mw 

0.6 Mw 

0.8 Mw 

8400 gal/min D,O at 
478°F 

1025 psig at pressur- 
izer 

478°F 

530°F 

85 

8-in. equilateral 
triangle 


42-75 

10-43 

60 

68 

1200 kw 

123 gal/min 

542°F 

8.3 x 10'3 neutrons/ 
(cm?)(sec) radial 

7.7 x 10!3 neutrons/ 
(cm*)(sec) axial 

5—8000 Mwd/ton 

Up to 2 hr 


10 Mw/min 


2 hr 


Less than 5 min 
2 hr (min.) 
13 Mw/min (max.) © 


0.72 av./max, (1.39 
max. /av.) 

0.833 av./max. (1.20 
max./av.) 


Reactor Piping: Process Tubes 


Material 

Design stress 

Dimensions ID (active zone) 
Wall thickness (active zone) 
Length 
Wall thickness (inlet end) 
OD (inlet end) 

Weight 

D,O flow area for empty tube 
(active zone) 

D,O flow area with cluster ele- 
ment, Mk-I, in tube 

D,O flow area with concentric 
element, Mk-II, in tube 

D,O flow/tube, maximum 

D,O velocity, empty tube (123 
gal/min) 


Zircaloy-2 

14,400 psi 

3.250 in. + 0.010 in. 

0.154 in. + 0.008 in. 

17 ft 5 in. 

0.235 in. 

2'f, in. 

85 lb 

8.3 sq in. (0.0576 
sq ft) 

3.56 sq in. (0.0246 
sq ft) 

2.27 sq in. (0.0158 
sq ft) 

123 gal/min 


4.75 ft/sec 





D,O velocity, 19-rod cluster ele- 
ment (123 gal/min) 

D,O velocity, concentric element 
(123 gal/min) 

Flow split in concentric element: 
Inner channel 
Middle channel 
Outer channel 

D,O inventory (85 tubes and noz- 
zles) 

Pressure drop across cluster 
element 

Pressure drop across concentric 
element 


Fuel Elements 


19-rod cluster elements: 
Length of active core in ele- 
ment 
Length of element 
Length of hanger 
Length of fuel assembly 
(hanger and fuel-element 
assembly) 
Bare rod diameter 
Sheath material 
Sheath thickness 
Core cross-sectional area (19 
rods) 
Total cross-sectional area (19 
rods) without wire wrap 
Core weight: 
Uo, 
Pu-Al 
Fuel-element assembly weight 
(19 elements) 
uo, 
Pu-Al 
Hanger weight 
Fuel-assembly weight: 
uO, 
Pu-Al 
Concentric elements: 
Length of UO, in element 
Length of element 
Length of hanger 
Length of fuel assembly 
(hanger and fuel-element 
assembly) 
UO, diameters: 
Rod OD (Mark-II C) 
Inner tube ID 
Inner tube OD 
Outer tube ID 
Outer tube OD 
Sheath material 
Sheath thickness 
UO, cross-sectional area 
(total) 
Total cross-sectional area 
UO, weight (10.2 g/cm’) 
Fuel-element assembly weight 
Hanger weight 
Fuel-assembly weight 
Maximum heat flux 


11.1 ft/sec 
17.7 ft/sec (av.) 


17% 

49.4% 

33.6% 

450 gal with 19-rod 
Mk-I element 


~T7 psi 


~15 psi 


7 ft 4 in. 
8 ft 3 in. 
6 ft 2 in. 


14 ft 5 in. 
0.504 in. 
Zircaloy-2 
0.030 in. 


3.76 sq in. 
4.76 sq in. 


122 lb 
32.8 lb 


160 lb 
70 lb 
40 lb 


200 Ib 
110 lb 


7 ft 4 in. 
8 ft 7 in. 
5 ft 10 in. 


14 ft 5 in. 


0.548 in. 
1.082 in. 
1.782 in. 
2.328 in. 
2.948 in. 
Zircaloy-2 
0.060 in. 


4.34 sq in. 

6.03 sq in. 

141 Ib 

195 lb 

40 lb 

235 lb 

400,000 Btu/(hr)(sq ft) 





(Table Continues) 
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Table IX-1 


(Continued) 





Primary Coolant System 


Heating load 
D,O flow rate 
D,O inventory (exterior of ring 
headers) 
Steam generator: 
Type and number 
Total surface area 
Tube material 
Total evaporation rate (full 
power) 
D,O velocity 
D,O pressure drop 
Pumps: 
Number and type 
Spares in circuit 
Operating head 
Capacity (each) 
Electrical power (each) 
Material 
Piping: 
Size 
Material 
D,O velocity 
Pressurizer: 
Pressurizing method 


Pressurizer size 
Control pressure 
Steam system: 
H,0 flow 
Steam generator pressure 


66.4 Mw 
8400 gal/min 


2200 gal 


1, U tube and shell 
6658 sq ft 
304 stainless steel 


214,000 lb/hr 
12 ft/sec (max.) 
14 psi 


3, mech, seal 

i 

110 psig 

4200 gal/min 

350 hp 

316 stainless steel 


14 in. sch. 100 
304 stainless steel 
24 ft/sec 


Helium bleed from hp 
tanks 

109 cu ft D,O 

1025 psig 


535 gal/min 
425 psia 


Moderator Coolant System 


Heating load 

D,O flow rate 

D,O inventory 

Temp. range 

Heat exchangers: 
Type and number 


Total surface area 
Pumps: 
Number and type 
Electrical power (each) 
Spares 
Operating head 
Capacity 
Secondary cooling water: 
H,O flow rate 


2.2 Mw 

1086 gal/min 
2544 gal 
137—149°F 


1, countercurrent 
shell and tube 
1000 sq ft 


3, mech. seal 
60 hp 

1 

110 psig 

600 gal/min 


750 gal/min at 68°F 
inlet 


Heavy Water (General) 


D,O inventory, total 
Purity specification 


8210 gal 
99.8% 





D,O density at 505°F 0.858 g/cm? (53.4 


lb/cu ft) 

D,O heat capacity at 505°F 1.16 Btu/(lb)(F) 
Resin quantity (IX-1) . 22 cu ft 
Resin operating temp. 130°F 
Total capacity, corrosion prod- 

ucts 19.5 g equivalent/cu ft 
Storage capacity (moderator 

storage tank and primary D,O 

storage tank) 7000 gal 


Reactor Core Components 


Calandria: 
Over-all length (top of vessel 
to bottom of vessel) Outside 9 ft 10'4 in. 
ID of moderator shell 7 ft 0 in. 
Wall thickness of shell 1% in. 
Over-all diameter of calan- 
dria, including dump cham- 
ber 11 ft 0 in. 


Material of construction Aluminum 
Design pressure 5 psig minimum 
Test pressure 6.25 psig 
No. of process-channel shroud 

tubes 85 
OD of process channel shroud 

tubes 4.25 in.* 
Wall thickness of shroud tubes 0.065 in.* 


Gas gap between shroud and 
process tubes 0.258 in. minimum 
No. of flux-monitor channels 13 


OD of monitor-channel tubes 1.5 in. 
Wall thickness of monitor- 

channel tubes 0.049 in. 
Material of monitor-channel 
tubes Aluminum 
No. of shim control channels 18 
OD of shim control channels 2.5625 in. 
Wall thickness of shim 

control-channel tubes 0.065 in. 
Volume of heavy water in 

moderator 284 cu ft 
Bottom drain line 8 in. 
Top overflow line 8 to 4 in. 

Reflector: 


Type of reflector used 
Length of radial reflector (top 
to bottom) 


Heavy-water tank type 


6 ft 74 in. 


Thickness of radial reflector 23% in. 
Over-all diameter of outer re- 
flector tank 11 ft 0 in. 
Total metal thickness between 
inner reflector and modera- 
tor 14 in. of aluminum 





*Center shroud tube: 6 in. OD; 0.085 in. thick. 


The containment vessel houses all the PRTR 





process areas with the exception of the storage 
pool mentioned previously. The vessel is an 
all-welded steel cylinder with a hemispherical 
dome and ellipsoidal bottom head; the diameter 


is 80 ft, and the over-all height is about 120 ft. 
The top of the rotating shield plug is approxi- 
mately at grade level, and the space inside the 
containment vessel below grade is divided into 
three cells by 5-ft-thick concrete shielding 
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walls. One of the cells contains the components 
necessary for reactor operation with the excep- 
tion of the ion-exchange resins, which are 
housed in a vault exterior to the containment 
vessel. The second cell contains experimental 
equipment, and the third cell is an instrument 
and hot shop cell. Adjacent to the experimental 
equipment cell is an air-cooled pit used for the 
examination of irradiated ‘fuel tubes and process 
tubes. 





Additional engineering data on the PRTR can 
be found in Table IxX-1. 


Reference 


1. N. G. Wittenbrock, P. C. Walkup, and J. K. Ander- 
son, eds., Plutonium Recycle Test Reactor Final 
Safeguards Analysis, USAEC Report HW-61236, 
Hanford Atomic Products Operation, Oct. 1, 1959. 








THE HOMOGENEOUS REACTOR EXPERIMENT NO. 2 





The Homogeneous Reactor Experiment No. 2 
(or HRT) was discussed in the December 1959 
issue of this journal, Vol. 3, No. 1. A detailed 
analysis of the power transients mentioned in 
that review has recently been published.’ Table 
X-1 lists characteristics of anumber of positive 


tive excursions in which the power increased 
more than 7 per cent, and there were over 70 
negative excursions in which the power fell by 
at least 7 per cent of normal power. The fre- 
quency of occurrence of the positive excursions 
over 7 per cent is shown in Fig. 14 as a function 


Table X-1 CHARACTERISTICS OF A NUMBER OF POWER EXCURSIONS! 











rags en Peak Rate of reactivity 
Run Pos AP/Py,; peak, period, reactivity addition, %%/sec 

Time Date No. Mw % sec sec addition, % Maximum Average 
0149 8-15-58 17-A 3.3 -—5.0 0.6 -7 —0.030 —0.070 —0.050 
0450 8-15-58 17-A 3.3 5.4 10 40 0.060 0.011 0.006 
0713 8-15-58 17-A 3.4 6.4 5 30 0.066 0.014 0.013 
0627 8-15-58 17-A 3.4 6.6 3 16 0.058 0.024 0.019 
0445 8-15-58 17-A 3.3 8.7 8 26 0.062 0.020 0.008 
0600 8-15-58 17-A 3.4 re | 4 16 0.066 0.026 0.016 
0708 8-15-58 17-A 3.2 15 4 15 0.082 0.064 0.020 
1700 9-4-58 17-B 2.2 31 3 6 0.13 0.079 0.043 
1103 9-7-58 17-B 2.2 32 4 11 0.16 0.046 0.040 
0135 9-6-58 17-B py 47 5 9 0.19 0.094 0.038 
0259 9-8-58 17-B 2.1 61 8 5 0.35 0.10 0.044 
1838 9-4-58 17-B 2.4 150 7 f 0.58 0.15 0.083 
1330 8-30-58 17-B 3.7 580 8 2 2.2 0.7 0.17 





Table X-2 CONDITION OF REACTOR! DURING 
RUNS 13, 14, 16, AND 17 








Run No. Condition 
13 D,O reflector; powers up to 1.5 Mw 
14 D,O reflector; powers up to 6.3 Mw; termi- 
nated by occurrence of hole in core vessel 
16 Blanket fuel concentration about 36 per cent 
of that in core; powers up to 5.8 Mw 
17 Corrosion-sample assembly present in core; 
blanket fuel concentration about 32 per cent 
of that in core; powers up to 5.2 Mw 
7-A First 20 days of run 17 
7-B Twenty days following addition of 3.44 moles 
of acid on August 20 
17-C Two days following addition of 5.02 moles of 


acid on September 10 





and one negative power excursion. Table X-2 
summarizes the condition of the reactor during 
the runs. Additional information may be obtained 
from the December 1959 Review. During runs 
14, 16, and 17, there were more than 800 posi- 
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of power level. A plot of the frequency of oc- 
currence of negative excursions is also shown 
in the reference, but there appeared to be little 
correlation to reactor power level. After the 
addition of acid for runs 17-B and C, the fre- 
quency of positive excursions decreased. Inad- 
dition, for 2 hr during run 17-B, the core 
solution-circulating pump was operated in re- 
verse, giving 40 per cent of normal forward 
fluid flow in the core. This point is also shown 
in Fig. 14. During run 17-B, the reactor was 
operated at average core temperatures from 
243 to 272°C at two power levels, 1.9 and 3.6 
Mw. At the lower power level, the excursion 
frequency was insensitive to temperature; 
whereas, at higher power level, the excursion 
rate increased with increasing temperature. 
The frequency numbers shown on the figure are 
averages taken over a relatively long period of 
time; the time between excursions was found to 
be quite irregular, ranging, for one run, from 
less than 1 min to as much as 141 min. 
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TOTAL POWER, MW 
Figure 14— The effect of total reactor power on fre- 
quency of positive excursions. O, run 17-A (before 
acid addition). 0, run 17-B (after acid addition). A, 
run 17-C (after second acid addition. Data represent 
434 excursions with the average core temperature 
near 270°C. 


A statistical analysis was made of the rela- 
tion between the relative frequency of excursions 
and their size. It was found that Eq. 1 expressed 
this relation for all but five of the excursions. 


f(x) = 6 e~8 (x x9) (1) 


where /(x) = relative frequency 
x =power increase/initial power = 
AP/P, 
64x) = constants for a given run* 


The five excursions that were eliminated from 
the correlation by means of statistical “rules” 
are, however, the four largest excursions in 
run 17 and the largest in run 16; these are not 





*For run 16, @ = 0.194; x9 = 4.5. 





statistically similar to the remainder of the 
excursions. 

To analyze the excursion data, a program es- 
timating the time variation of reactivity associ- 
ated with a power excursion was coded for the 
IBM-704. It was stated that the reactivity ob- 
tained with the code is that reactivity associated 
with effects other than those produced by changes 
in delayed-neutron concentration and reactor 
temperature. Excursion histories were then 
analyzed using the computer program. The re- 
activities obtained were equivalent to moving 
several grams of uranium from outside the 
core into its center, or about 10 g from the 
upper screen to an average position in the core. 
The weights, of course, depend on which detailed 
excursion was analyzed; the numbers cited are 
for a peak reactivity addition of about 0.08 per 
cent. In general, the reactivity changes during 
the excursions could be satisfactorily explained 
by the movement of a relatively few grams of 
fuel of high density (solid) suddenly into the 
core where it would circulate and dissolve or 
pass out of the core. Repeated peaks in the ex- 
cursion traces could be associated with the re- 
turn of uranium tothe core after passage through 
the external piping, since the time between peaks 
coincided closely with the transit time in the 
external flow loop. 

A peculiar power trace that occurred during 
run 17 was analyzed. This trace exhibited 
nearly a step increase in reactor power followed 
by a rapid decrease to below the original power 
setting and a gradual ramp back to the set- 
point; at least five such traces were observed 
in the power trace records. This behavior could 
be explained by uranium falling from the top of 
the core to the bottom of the core, and itis 
postulated that some, if not all, of these traces 
were caused by uranium-bearing corrosion 
specimens falling from their original position 
at the top of the core to the screen area. (See 
the December 1959 issue of Power Reactor 
Technology, Vol. 3, No. 1, for additional details 
on the corrosion specimens present during run 
17.) The recovery of reactivity might then be 
caused by dissolution of the uranium from the 
specimens due to better flow patterns around 
the screens or some other mechanism. This 
might explain the four excursions in run 17 
which would not fit the statistically determined 
frequency equation, Eq. 1; but the largest ex- 
cursion in run 16, of course, could not be ex- 
plained by a movement of corrosion specimens, 
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since presumably none were present in the re- 
actor during run 16. 

In addition to the excursions previously dis- 
cussed, the HRT also exhibited a power varia- 
tion termed “power fluctuations.” These were 
small changes in power, usually less than 1 per 
cent dP/P,. Shown in Fig. 15 is the parameter 
dP/P, versus reactor power for both the power 
fluctuations and power excursions. Figure 16 
shows the frequency of positive peaks for the 
fluctuations and also a plot of Eq. 1. From the 
two curves it is evident that the fluctuations, 
although they occur more frequently (for “low” 
powers), have considerably lower amplitude. 
Data presented in the reference show that the 
fluctuations are apparently insensitive to core 
temperature. A hypothesis was advanced that 
the power fluctuations were due to slugs of 
“cold” fuel solution moving in convective flow 
within the core and/or blanket. The PDQ code 
was utilized to examine the reactivity changes 
associated with the movement of various amounts 
of cold fuel solution about the core, and, in 
general, the results appeared consistent with the 
behavior of the reactor. The operating experi- 
ence with the HRT tends, however, to disprove 
this hypothesis. If convective flow were respon- 
sible, then the changed convective conditions, 
during the 2 hr in run 17-B when low flow was 
experienced, should have resulted in a change 
in frequency of the fluctuations. No such change 
was observed during the low-flow portion of 
run 17-B. A run was made on an analog ma- 
chine to determine whether ‘the small power 
fluctuations are a natural property of the circu- 
lating system of the HRT; however, the results 
showed none of the oscillatory character of the 
fluctuations, and stable behavior was observed. 
At the time of writing, the source of the fluc- 
tuations had not been identified. 

The power-fluctuation traces shown in refer- 
ence 1 bear a remarkable similarity to the 
power records of boiling-water reactors, shown 
in reference 2, page 89. Void formation in the 
HRT was studied using a one-dimensional, two- 
region diffusion code, but the application to the 
power fluctuations is not discussed in refer- 
ence 1. 

Since the discovery of the hole in the HRT in 
April 1958, the reactor has been operated as a 
one-region machine. The latest run, lasting for 
105 days of continuous operation, was termi- 
nated when a second hole developed in the zir- 
conium core tank.* The new Opening measures 
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Figure 15— The effect of core power level on (a) av- 
erage size of positive excursions and on (b) power 
fluctuations. 
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Figure 16— The effect of power increase on frequency 
of excursions and fluctuations. 





THE HOMOGENEOUS REACTOR EXPERIMENT NO. 2 


about '/, by 1% in. and is located just below the 
equator of the spherical core tank. It is now 
planned to repair both holes in the core tank of 
the HRT and resume experimental operation 
with the flow in the core reversed, since it is 
thought that this will eliminate local stagnant 
areas which are believed responsible for both 
penetrations. ’ 
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BOILING-WATER REACTORS 





Modifications of the EBWR 


The EBWR, in its original form, has been ade- 
quately described in the literature.'? A more 
recent publication deals with the modifications 
of the reactor to permit safe operation at 100 
Mw(t) and evaluates the possible hazards asso- 
ciated with such operation.° 

The EBWR has so far been successfully op- 
erated at a power of 61.7 Mw(t) with a core 4 ft 
in diameter by 4 ft in height. A combined ana- 
lytical and experimental investigation of the 
stability of boiling-water reactors (see the 
December 1959 issue of Power Reactor Tech- 
nology for additional details) has indicated that 
stable operation of the EBWR at 100 Mwi(t) is 
probably feasible if the core diameter is in- 
creased to 5 ft; this size of core can be installed 
in the existing EBWR pressure vessel. The 5- 
ft-diameter core will be composed of the follow- 
ing set of fuel-element assemblies: 


112 partially enriched (1.44 per cent) uranium 
(existing) 
7 natural uranium (existing) 
1 source element 
28 highly enriched uranium (new) 


148 total 


The elements marked “existing” represent ele- 
ments already fabricated for EBWR operation 
and are of two different types; some are of 
“thin plate” design and some are of “thick 
plate” design. These existing elements are 
composed of six fuel plates with an active length 
of 4 ft. The new fuel assemblies are of the rod 
type, wherein a 0.325-in.-OD fueled rod is clad 
with a *%-in.-OD Zircaloy-2 tube. The fueled 
rod is a dispersion of highly enriched U,O, in 
an X-8001 aluminum matrix; 49 such rods are 
contained in the fuel-element assembly. The 
active length of the new elements is 5 ft. Since 
the remainder of the core is only 4 ft in height, 
the new elements are of the reversible type, 
being symmetrical with respect to length. Re- 
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movable upper and lower end fittings are pro- 
vided to facilitate this reversible feature. 

The core loading for 100-Mw operation em- 
ploys the highly enriched elements as “spikes.” 
The 28 spikes will be located in a square, sur- 
rounding 36 of the shorter, existing fuel ele- 
ments. The location of the spikes in this posi- 
tion results in a considerable degree of radial 
power flattening. Control of the reactor is ef- 
fected by nine cross type control rods with a 
10-in. blade width. The rods are followed by a 
fueled plate follower of similar span; this is a 
departure, of course, from the old EBWR 
control-rod follower. The absorber section, 2 
per cent boron in stainless steel, is 58 in. in 
length, whereas the fueled section, highly en- 
riched uranium in a 304 matrix, clad with 304 
stainless steel, has a length of about 36 in. 
Addition of 1 per cent boron-stainless steel 
strips to the sides of the spike elements is 
being considered should additional control be 
necessary. Since the old EBWR rod had an ab- 
sorber section length of 46 in., redesign of the 
control-rod-drive mechanisms to accommodate 
this longer rod has been necessary. 

Some changes to the pressure-vessel inter- 
nals are necessary. To achieve a higher ratio 
of steam production rate to average steam void 
fraction in the core, increased recirculation of 
the water through the core is desired. To 
achieve this increased flow, it is proposed to 
put a core riser or chimney over the new core. 
This riser produces a greater head difference 
than was previously available between the mix- 
ture leaving the core and the water inthe down- 
comer space. In order to accommodate the in- 
creased flow of steam from the pressure vessel, 
the outlet steam-collecting system is to be 
modified. This modification includes a larger 
steam outlet nozzle (increased from 6 to 10 in.) 
and a baffle system designed to permit operation 
with a higher water level in the reactor and still 
provide for steam-water separation. A larger 
feed-water inlet nozzle and distribution ring are 
also needed. The present turbogenerator will 
remain unchanged and is capable of absorbing 
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Table XI-1 


PLANT PARAMETERS, MODIFIED EBWR® 








Genera 


Type 
Heat output, kw 


Gross electrical output, kw 
Operating pressure, psig 
Operating temp., °F 
Coolant 

Moderator 

Fuel 


Core 


Active diameter, ft 

Active height, ft 

Total uranium content, tons 

Total u2% content, kg 

Average water-to-uranium 
volume ratio 

Structural material 


l 


Direct-cycle boiling 
reactor 

20,000 kw for turbine; 
80,000 kw for heating 

5000 

600 


5 
6.4 
65-90 


3.64 
Zircaloy-2 and 
aluminum 


Fuel Assemblies 


Total number in core 


Cross-sectional dimensions, in. 
No. of plates per assembly 

No. of rods per spiked assembly 
Side plates 


Cladding 


End fittings 
Composition of fuel ‘‘meat’’ 


Composition of spikes 


Thickness of meat in thin 
plates, in. 

Thickness of meat in thick 
plates, in. 

Thickness of meat in rods, in. 

Width of water channels 
between thin plates, in. 

Width of water channels 
between thick plates, in. 

Power generation in average 
assembly, Mw 


Nuclear Data 


Average thermal flux 


Neutron lifetime, sec 

Reactivity to control, % k 

9 control rods with 8 fuel 
followers, %k 


147 (~20% of assem- 
blies are spiked) 

3.75 x 3.75 

6 

49 

0.060-in.-thick 
Zircaloy-2 

0.020-in.-thick 
Zircaloy-2 

Stainless steei 

93.5% U, 5% Zr, 
1.5% Nb (by weight) 

Zr, Al structure with 
U30, (>90% enriched) 


0.174 


0.239 
0.325 diameter 


0.428 
0.360 


=0.68 


3 x 10!3 neutrons/ 
(cm?)(sec) 

6 x 107° 

~<=13 


~15 


Heat Transfer and Fluid Flow 


Average power density in core 
coolant: 
Spiked zone, kw/liter 
Unspiked zone, kw/liter 


135 
115 








Maximum power density in core 
coolant: 

Spiked zone, kw/liter 
Unspiked zone, kw/liter 

Steam flow (at 100 Mw), lb/hr 

Average steam voids in heated 
channels, % 

Average steam voids in total 
moderator, % 

Average exit-steam quality, % 

Pounds of water circulated per 
pound of steam 

Subcooling at core inlet, °F 

Nonboiling height of core, % 

Average coolant-inlet velocity, 
ft/sec 

Average heat flux: 
Spiked zone, Btu/(hr)(sq ft) 
Unspiked zone, Btu/(hr)(sq ft) 

Maximum heat flux, 
Btu/(hr)(sq ft), based on 
maximum -to-average heat 
flux ratio of 3.1 

Estimated burnout heat flux, 
Btu/(hr)(sq ft) 

Maximum surface temp., °F 
Average fuel-centerline temp.: 
Thick plate (enriched), °F 

Thin plate (enriched), °F 
Rod (spike), °F 

Maximum fuel-centerline temp.: 
Thick plate (enriched), °F 
Thin plate (enriched), °F 
Rod (spike), °F 


418 
356 
300,000 


18.9 


16.5 
3.2 


32 
12 
40 (av.) 


6.3 


136,000 
157,000 


485,000 


>1,000,000 
514 


579 
567 
558 


767 
724 
694 


Control Rods 


Total number 

Spacing, in. 

Shape 

Material 

Thickness, in. 

Boron content, wt.% 

Penetration of absorber into 
core, in. 

56-in. travel time, sec 

Maximum withdrawal rate: 
in./sec 
% k/sec 

Strength of 9 rods, % k 


9 

12.75 x 12.75 
Cruciform: 10 by 10 in. 
Boron-stainless steel 
0.250 

2 


Pressure Vessel 


Diameter (inside), ft 

Height, ft 

Working pressure, psig 

Design pressure, psig 

Thickness of cylindrical portion, 
in. 

Material 


Cladding 


Thermal shield 


2° 

SA-212 Grade B boiler 
plate 

0.1-in. type 304 
stainless steel 

1-in. 18-8 stainless 
steel containing 1% 
boron 





(Table Continues) 
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Table XI-1 


(Continued) 





6, 50, 700, 725, 750, 
and 775 

Total weight, tons 60 

Weight of contained water, tons 14 


Relief valve settings, psig 


Power Plant 
7 


Condenser pressure, in. Hg abs. 2h 
Flow rate to cooling tower, 
gal/min 13,650 (max.) 
Feed-water flow rate, gal/min 600 
Feed-water temp., °F 120 
Generator output voltage, volts 4160 
Steel Building 
Diameter, ft 80 
Height, ft 119 


Net volume, cu ft ~400,000 








Design pressure (internal), psig 15 
Design pressure (external), psig 0.5 
Maximum leakage rate at 15 
psig, cu ft/day . 1000 
Type steel ASTM A-201 Grade B 
Firebox Quality 
% (cylindrical) 
% (elliptical) 
% (hemispherical) 
2 ft (min.) concrete 
1-ft-thick concrete 


Sides of tank, inches thick 
Bottom of tank, inches thick 
Top of tank, inches thick 
Lining of tank below grade 
Lining of tank above grade 


Capacity of water reservoir, gal 15,000 
Flow rate of sprinkler system, 

gal/min 1000 
Ventilation air flow rate, 

cu ft/min 3000 


Air-lock doors 
Bolted and gasketed 
door 


Personnel access 
Freight access 





20 Mw(t); the remaining power will be either 
sent into the steam distribution system at Ar- 
gonne or dumped to the air via air-cooled 
steam-condensing units. To prevent any possi- 
bility of contamination of outside steam lines, a 
secondary steam generator is postulated. Other 
design detail on the modified reactor is givenin 
Table XI-1. 

The capability of the modified reactor to sup- 
ply 100 Mw of heat isnotassured. Before modi- 
fication, the reactor was operated at amaximum 
power of about 62 Mw, and extrapolation of 
oscillation tests indicated that the reactor might 
reach the limit of stability at about 68 Mw. If 
the naive assumption were made that average 
power density is the quantity which determines 
stability, then the limiting power output for the 
reactor would be nearly proportional to the 
number of fuel elements (neglecting the extra 
foot of length of the new elements). This as- 
sumption would indicate a maximum power of 
68 x 147/114 = 86 Mw. The most important fac- 
tors which will modify this estimate are (1) the 
extra circulating head provided by the chimney 
and (2) the changes in radial power distribution 
and in effective void coefficient of reactivity due 
to the changed fuel loading. 

Aside from the question of stability, it is 
possible that a deterioration of the effectiveness 
of steam-water separation at high steaming 
rates might limit the maximum output. 

If steam bubbles were entrained for a signifi- 
cant distance in the downcomer, the net driving 
head would be reduced. The resulting reduction 
in the recirculation flow rate would increase the 


core steam volume fraction and possibly hasten 
the onset of unstable operation. Equally as 
severe is the entrainment of water in the satu- 
rated steam flowing from the reactor. The use 
of the core riser decreases the volume availa- 
ble for disengagement, and the increased steam 
production rate increases the velocity in the 
steam dome to a superficial steam velocity of 
1.89 ft/sec. The degree of steam separation that 
will take place under these rather adverse con- 
ditions is not amenable to calculation. Itis true, 
however, that, even if the reactor were to de- 
liver steam of too poor a quality to admit di- 
rectly to the turbine, it is, inprinciple, possible 
to utilize external steam separation. 

Since publication of reference 3, the modifica- 
tions have been completed, and the EBWR has 
achieved criticality again. An extensive period 
of experimentation is planned to assess the ef- 
fects of .each of the structural modifications 
before very high power operation is attempted. 


Advanced Boiling-water Design 


A recent report’ presents the results of a 
relatively brief conceptual design study by Ar- 
gonne National Laboratory for the purpose of 
arriving at a design of a 50-Mw(e) prototype for 
a large boiling-water reactor of advanced de- 
sign. An engineering study of essentially the 
same reactor has also been made by Sargent & 
Lundy.° 

There has been rather general agreement that 
a profitable direction of development for the 
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Table XI-2 CORE DESIGN AND PERFORMANCE DATA, PROTOTYPE BCILING-WATER REACTOR‘ 








Fuel Description 


Fuel material UO, 

Pellet diameter, in. 0.410 

Gap material He (2 atm. pressure) 
Gap thickness, nominal (cold), ir. 0.0015 
Clad material Zircaloy-2 
Clad thickness, in. 0.025 
Fuel-rod OD, in. 0.463 
Fuel-section length, in. 40.5 

Axial gap, in. 1.75 
Active-zone length, in. 81 
Triangular lattice pitch, in. 0.6 


Core Description 


Active-zone length, in. 82.75 
Active diameter, ft 5.0 
Fuel-assembly shape Hexagon 
Lattice pitch, in. 5.125 
Rods per assembly 61 
Assemblies per core (including 

control rods) 121 
Rods per core 7381 
No. of control rods 19 
Power density, kw/liter 47 


Core Materials 


Weight of UO), lb 29,459 
Volume Zircaloy-2/volume inside 

clad 0.312 
Volume H,O + steam/volume 

inside clad 1.52 
Average density UO,/unit volume 

inside clad, g/cm® 10.0 

Physics Data 

Fuel enrichment (feed), at.% 2.65 
Fuel enrichment (discharge), at.% 1.60 


Pu?*® and Pu‘! content (discharge), 
at.% 

€ 

p 

tT, cm? 

B’, cm 

Rog (initial, hot, voided Xe + Sm 
poisoned) 

dk (Xe, Sm), % 

6k (initial to equilibrium), % 

5k (voids), % 

6k (between discharge points), % 

dk (required, hot), % 

6k in rods (hot), % 

Burnup at discharge, Mwd/ton U 

Total maximum-to-average flux 


2 


0.64 
1.05 
0.71 

90 
0.00104 


1.046 
3.7 
2.7 
8.9 
4.8 
18.5 
38 
12,000 
4.90 


Hydraulic and Thermal Data 


Total reactor power (design), Mw 
Average power/assembly (design), 
Mw 
Total heat-transfer area, sq ft 
Total core flow area, sq ft 
Average heat flux, Btu/(hr)(sq ft) 
Maximum heat flux, Btu/(hr)(sq ft) 
Fraction nonboiling length 
Void fraction in boiling zone 
Void fraction (uniform flux): 
Total 
Exit 
Steam exit quality 
Recirculation ratio 
Slip ratio 
Total flow rate (design), gal/min 
Inlet subcooling, °F 
Steam rate (design), lb/hr 
Average core inlet velocity, ft/sec 
Core pressure drop (water), ft 
Average fuel-centerline temp., °F 
Maximum fuel-centerline temp., °F 


165 


1.36 
6,030 
10.75 
93,500 
369,000 
0.25 
0.267 


0.20 

0.445 

0.0575 

16.4:1 

1.57 

30,000 

9.5 (12.5 Btu/lb) 
644,800 











boiling-water reactor would be toward the 
achievement of higher power density inthe core, 
provided this could be done without the addition 
of very expensive design features. It is, of 
course, probable that higher power density could 
be obtained in a straightforward way simply by 
increasing the subdivision of the fuel (smaller 
fuel-rod diameter) and by providing sufficient 
coolant flow velocity by forced circulation to 
keep the amount of reactivity invested in steam 
voids to something like a conventional level. 
Although this approach should require little or 
no further research and development, it has the 
disadvantage that the greater fuel subdivision 
increases the fuel fabrication costs per unit 
mass of fuel, and the provision of high circula- 
tion rates will be reflected in pump and piping 


costs. An alternate approach is to attemptto 
pack more fuel into a given core volume (i.e., 
increase the fuel-to-water ratio}, to attempt to 
improve the spatial distribution of power, and 
to accept, if necessary, a high reactivity invest- 
ment in steam void in order to avoid excessive 
coolant-circulation rates. The Argonne design 
lies in the direction of the second approach. It 
utilizes a high (for boiling reactors) UO,-to- 
water ratio and accepts the accompanying high 
steam coefficient of reactivity, which, along 
with the coolant flow rate, results in a calcu- 
lated investment of 8.9 per cent Ref in steam 
voids at rated power. The power distribution 
achieved is not outstanding, partly because of the 
skewing of axial power distribution by the axial 
distribution of steam voids and partly because 
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Table XI-3 PRESSURE-VESSEL AND CONTROL-ROD- 
DRIVE DATA, PROTOTYPE BOILING-WATER REACTOR‘ 





Pressure Vessel 


Inside diameter, ft sl, 
Inside length, ft 31'/, 
Main side-wall thickness, in. 3.25 
Material: 
Shell plate SA-212-B 
Flanges SA-105-II 
Nozzles SA-106-B 
Shipping weight, lb 210,000 
Operating pressure, psia 1000 
Design pressure, psia 1100 
Operating temp., °F 545 
Design temp., °F 650 
Vessel liner: 
Thickness, in. 0.187 


Type 304 stainless 
steel 


Material 


Nozzles 

Coolant inlet: 

Inside diameter, in. 13.938 

No. required 4 

Flow velocity, ft/sec at 30,000 

gal/min 16 

Coolant outlet: 

Inside diameter, in. 20.938 

No. required 4 

Flow velocity, ft/sec 16 

Control-rod Drive 
No. required 19 
Type of operation Gang-hydraulic- 
mechanical] 

No. of gangs 6 
No. of drives per gang 3 
Center rod drive 1 
Normal driven speed, in./min 12 
Fast shutdown time, sec 3 
Control-rod weight (including fuel), 

lb 375 





the fuel management program selected gives a 
relatively poor radial distribution. 

Certain arbitrary goals affected some of the 
design decisions. In particular, a nominal value 
of 50 Mw/liter was selected as the average core 
power density, and the design was fitted to this 
specification. A length-to-diameter ratio of 1.30 
was selected for the core in order touse rather 
long fuel elements which would be more nearly 
prototypes for the elements of a large reactor. 
The decision for a high fuel-to-water ratio led, 
directly or indirectly, to a number of other de- 
sign decisions. A triangular fuel lattice was 
chosen because fuel-rod spacings become un- 
comfortably close in a “dry” square lattice. 
The triangular spacing led naturally to the se- 


lection of a hexagonal shape for the fuel-element 
assemblies. In a dry core the reactivity ef- 
fectiveness of thermal-neutron absorption by 
blade type control rods is reduced, and the 
power peaks due to water channels are in- 
creased. These effects led to the selection of 
control rods in which the absorbing section 
consists of a hexagonal shell of absorbing mate- 
rial, filled with water. A fueled follower is 
provided to fill the water gap as the absorbing 
section is withdrawn. This construction is cal- 
culated to provide a negative reactivity capacity 
of 38 per cent Resp in 19 control rods. 

Other features of the design are: the upward 
direction of control-rod motion for reduction of 
reactivity (similar to Dresden) to improve the 
axial power distribution; the proposed fuel 
management program, which utilizes small batch 
reloading into the center of the reactor, with 
outward radial shift of the fuel elements as their 
exposure increases. Inversion of fuel elements 
is also employed to equalize the axial distribu- 
tion of burnup. To accomplish the latter objec- 
tive, fuel is handled in the form of hexagonal 
cartridges, each cartridge consisting of two 
half-core-length fuel assemblies placed end to 
end in a hexagonal holder. The holder consists 
of a hexagonal cage with end fittings to engage 
the top and bottom core grids. To invert fuel, 
the cartridge is removed from the reactor tothe 
fuel storage pit. The two assemblies are re- 
moved from the holder and are replaced with 
their positions interchanged, the assembly ini- 
tially in the top position being changed to the 
bottom position. This interchange, in effect, 
turns the fuel “inside out,” placing the fuel 
which was originally near the ends of the core 
near the center of the core. After the inversion, 
the cartridges are placed in the reactor in a 
new radial position. The loading of new fuel into 
the center of the reactor minimizes the enrich- 
ment requirement for a given fuel lifetime but 
yields a relatively poor radial power distribu- 
tion. The scheme does, however, fit in well with 
the axial inversion program, since the fuel car- 
tridges, after the inversion which places the 
most reactive fuel near the center plane of the 
reactor, are occupying radial positions of re- 
duced power density. 

The characteristics of the reactor core are 
summarized in Table XI-2, and the pressure- 
vessel and control-rod-drive data are sum- 
marized in Table XI-3. The rated coolant flow 
at full power is 30,000gal/min. This is supplied 
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by two pumps, each of 25,000 gal/min capacity, 
for experimental flexibility. Since the fuel- 
assembly holders are simply open cages, there 
is no partitioning of the coolant flow between 
assemblies, and the possibility of some lateral 
flow of coolant exists. It is stated that the re- 
sistance to lateral flow is about seven times the 
resistance to axial flow. The coolant pressure 
is 1000 psia, and the corresponding steam satu- 
ration temperature is 544°F. The reactor is 
designed for direct-cycle operation; steam- 
water separation is accomplished in external 
drums. The feed-water return temperature is 
350° F. It was stated that a higher design pres- 
sure might be optimum for the prototype reac- 
tor but that the pressure of 1000 psi was se- 
lected as nearer the optimum for a very large 
plant. 


Boiling Reactor Bibliography 


A bibliography of selected literature on 
boiling-water reactors has been issued by the 
Technical Information Service Extension of the 
Atomic Energy Commission. It is listed as 
reference 6. 
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Design studies of various reactor types were 
made in connection with the formulation of the 
AEC Ten-year Program for Civilian Power Re- 
actor Development. Reference 1 is such astudy 
by Atomics International, covering an investiga- 
tion of the effects of core design parameters on 
power costs for a300-Mw(e) organic-moderated 
reactor. The approach used in the study is, 
first, to define a reference design for the reac- 
tor, which, although not optimized, is expected 
to be somewhere near the optimum design. The 
effects of variation of specific design parame- 
ters are then investigated by varying those 
parameters individually (along with any de- 
pendent parameters) around the reference de- 


Table XII-1 


sign point. The effects of the variations are 
evaluated in terms of the change which they 
produce in the total power-generation cost. 

An earlier design study of a large organic- 
moderated reactor had been made as one of the 
studies carried out under the authorization of 
Public Law 85-590 (Sec. 101) and was reviewed 
in the September 1959 issue of Power Reactor 
Technology, Vol. 2, No. 4. The later study’ 
differs primarily in that it does not employ 
nucleate boiling of the organic coolant and in 
that it employs uranium—3.5 per cent molybde- 
num fuel clad with finned aluminum, rather than 
uranium dioxide jacketed in SAP (sintered alu- 
minum powder). It is stated that the design is 


ORGANIC-MODERATED REACTOR! 





Reference Design Data 


Organic moderated 
and cooled 


Reactor type 





Gross electrical power, Mw 314.4 
Gross cycle efficiency, % 31.1 
Thermal power to primary coolant, 

Mw 1010 
Net electrical power, Mw 291.5 
Net cycle efficiency, % 28.9 
Thermal power leakage from core, 

Mw 10 
Total thermal power, Mw 1020 
Coolant Santowax R 
Coolant inlet temp., °F 524 
Coolant outlet temp., °F 600 
Total fuel elements 372 


Type of fuel element 
Fuel (aluminum clad) 


Initial enrichment, at.% 


Finned plates 
Uranium —3.5% 

molybdenum alloy 
2.1 


Tot2i lattice spaces in core 408 
Core diameter, in. 164 
Core height, in. 156 
Core operating pressure, psi 90 
Total weight of fuel in core, kg 120,400 
Total weight of uranium in core, kg 115,600 
Average burnup, Mwd/metric ton of 

uranium 7,000 
Peak burnup, Mwd/metric ton of 

uranium 11,000 
Initial reactivity (first core), % 9.5 
Depleted enrichment 1.38 
Plutonium in spent element, g/kg of 

uranium 3.67 
Initial conversion ratio 0.6 








Moderator-to-fuel ratio 3.42 
Control-rod type Cruciform 
No. of control elements 73 
Total worth of control rods, % 15 
Peaking factors: 

Over-all 4.08 

Axial 1.57 

Radial 1.5 

Element (adjacent to rod) 73 
Hot-channel factors: 

F. 1.08 

Fy, 1.30 

Fo 1.43 
Maximum coolant velocity, ft/sec 15 
Average coolant velocity, ft/sec 10 
Maximum permissible fuel temp., 

“Fe 1200 
Maximum cladding surface temp., 

by 750 
Maximum fuel-element power, Mw 4.1 
Core pressure drop, psi 80 
Maximum heat flux, Btu/(hr)(sq ft) 94,000 
Burnout heat flux, Btu/(hr)(sq ft) 500,000 
Maximum specific power, kw(t)/kg 

of uranium 35.7 
Average specific power, kw(t)/kg 

of uranium 8.8 
Core power density, kw(t)/liter 18.8 
Mode of heat transfer Forced convection 
Total coolant flow, lb/hr 95,000,000 
No. of primary loops 4 
Steam pressure, psig 550 
Steam temp., °F 580 


Containment 


Totally enclosed 
primary system 
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based upon current technology and that it is 
capable of construction today with no additional 
research anddevelopment. The estimated power 
cost for a 300-Mw(e) plant is 9.02 mills/kw-hr, 
which is higher than that estimated in the origi- 
nal large-reactor design but lower than that 
cited for a 300-Mw(e) organic-moderated plant 
in the AEC summary’ of technical and economic 
status as of 1959 (11.45 mills/kw-hr). Most of 
the latter difference appears to be due to the 
use of an average fuel exposure at discharge of 
7000 Mwd/metric ton of uranium in reference 1 
rather than a value of 4500 Mwd/metric ton of 
uranium as in the AEC report. 

The main characteristics of the reference de- 
sign are summarized in Table XII-1. The reac- 
tor core is 13.7 ft in diameter and 13 ft high. It 
is made up of 372 fuel elements arranged in a 
close-packed lattice. Each element, in cross 
section, consists of 10 metal fuel plates, of 
uranium--3.5 per cent molybdenum alloy, clad 
with a finned aluminum sheath. The 10 fuel 
plates are enclosed within a square stainless- 
steel box having a wall thickness of 0.060 in. 
and an outside dimension of 6.2 in. The ele- 
ments are 180 in. long and have anactive length 
of 156 in., made up of 20 subassemblies, each 
15.6 in. long and consisting of 10 pairs of half- 
width plates. The fuel meat thickness is 0.117in. 
The cladding has a minimum thickness of 0.02 
in. and provides 51 fins, 0.20 in. high, on each 
plate. The inside of the cladding is coated with 
a 0.001-in. nickel diffusion barrier. Thecoreis 
surrounded by a 6-in. organic reflector and two 
annular steel thermal shields, 6 and 1.5 in. in 
thickness. The core vessel is 17.2 ftindiameter 
and 60 ft high and has a wall thickness of 2% in. 

The parameter perturbation study was limited 
to parameters affecting the reactor cost, pri- 
mary heat-transfer equipment cost, net fuel 
cost, and organic-coolant makeup cost. These 
comprise about 65 per cent of the cost of nuclear 
power from the reference plant. The detailed 
parameters varied are as follows: core and fuel 
geometry; fuel-plate thickness; resonance es- 
cape probability; stable fission-product poisons; 
fuel exposure; fuel-element box thickness; gross 
axial and radial flux-peaking factors; element 
flux-peaking factor; hot-channel factors; fuel- 
element surface temperature; high boiler con- 
tent in the coolant; core length-to-diameter 
ratio; coolant velocity; and average specific 
power. The results of the study are presented 
in graphical form, an example of which is given 
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Figure 17— Effect of maximum fuel-element surface 
temperature on differential power cost (adapted from 
reference 1). 


in Fig. 17, for the effect of variation of maxi- 
mum fuel-element surface temperature. 


The authors present the following conclusions:! 


The most important cost determining factors 
among the parameters investigated are the thermal 
parameters of fuel-cladding surface temperature 
and_hot-channel factors. Each of these variables 
increases in importance as the fuel exposure in- 
creases. 

If the moderator-to-fuel ratio (M/F) is changed 
by varying the fuel-plate thickness while keeping 
the number of plates constant, there is an economic 
optimum M/F which is dependent on fuel exposure. 
For a fuel exposure of 5000 Mwd/metric ton of ura- 
nium, the economic optimum M/F is 4.5. For an 
exposure of 7000 Mwd/metric ton of uranium, the 
optimum M/F is around 5.0, while for 9000 Mwd/ 
metric ton of uranium, the optimum is in excess of 
5.0. 

Considerable savings can be realized by flattening 
the radial and fuel-element neutron-flux distribu- 
tions, with more savings accruing for higher fuel 
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exposures. Less savings can be realized by flatten- 
ing the axial neutron-flux distribution, with little 
difference for various fuel exposures. 

The resonance escape probability has little effect 
on power costs at low fuel exposures. This parame- 
ter becomes more important as the fuel exposure 
increases. 

The average fuel exposure need not be extended 
beyond 8500 Mwd/metric ton of uranium for the type 
of metal-fueled reactor considered in this study. 

The quantity of structural steel in the fuel ele- 
ments has a pronounced effect upon fuel costs and 
should be kept to a minimum. 

The optimum coolant velocity for a reactor of 
this type is 11.2 ft/sec for the hot channel. 

Lattice spacing, core length-to-diameter ratio, 
high boiler content between 20 and 40 wt.%, and 
stable fission-product poisoning between 40 and 60 
barns/fission have little effect on power costs. 


In general, the effect of varying the individual 
parameters changed the power cost by some- 
thing less than 1 mill/kw-hr, e.g., Fig. 17. 
Several of the results deserve further com- 
ment. The power cost was surprisingly insensi- 
tive to moderator-to-fuel ratio. The moderator- 
to-fuel ratio in many reactor types has an 
important effect on power density, and high 
power density is often rather effective in de- 
creasing capital cost. Presumably, this insen- 
sitivity in the case of the organic reactor is 
partly due to the relatively low pressure of the 
coolant and to the absence of stored mechanical 
energy in the coolant system. Thus the major 
effects of large reactor size, in increasing the 
costs of the reactor vessel and of the contain- 
ment structure, are considerably tempered. 
Although the possibility of decreasing power 
cost by extending fuel exposure may be an 


academic one in a reactor employing metallic 
fuel elements, the conclusion, that increases in 
fuel life beyond 8500 Mwd/metric ton of ura- 
nium do not decrease power cost, is at first 
glance surprising. Reference 1 points out that 
the minimum occurs at this relatively low value 
because the ground rules for computing fuel cost 
arbitrarily specify that inventory charges are to 
be based on 1.6 times the fuel in the reactor. 

The rather low value for the optimum coolant 
velocity results partly from the circumstance 
that, as coolant velocity is increased, the thick- 
ness of the stainless-steel fuel-element box 
must be increased to withstand the flow pres- 
sure drop. The effect of the stainless-steel ab- 
sorption on neutron physics produces an adverse 
effect on the power cost. 

Although the effect of varying the high boiler 
residue content is investigated in the report, 
there is no discussion of the effect of variations 
in the basic rate of organic decomposition on 
power cost. 

Reference 1 points out that the variations 
studied are primarily variations in design pa- 
rameters and predicts that larger improve- 
ments can be achieved through extensions ofthe 
existing technology. 
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